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Seven bands of ammonia between 1620 and 1450A have been photographed in absorption 
with a dispersion of 1.3A/mm. These are the 0,2 —0,8 bands of the second electronic state above 
the normal. Each band has one P and one R branch, and a broad intense line as a Q branch. 
Within the limit of observational errors, which were large in this region, the moment of inertia 
A’ is only slightly larger than A’’, and C’ is about equal to C”, so that the dimensions are not 


greatly changed in the excited state. 


REPORT has been made recently of the 

absorption spectrum of ammonia in the 
ultraviolet between 2300 and 850A.' The spec- 
trum was examined with a vacuum grating 
spectrograph with a dispersion of about 8.5A/ 
mm. Some bands in the region 1700—1400A were 
found to possess considerable fine structure, 
which appeared partly resolved. Since the 
moments of inertia appeared to remain small in 
excited electronic states, almost complete resolu- 
tion seemed possible at the highest dispersion 
available. It was hoped that an analysis of the 
rotational structure could be made, in spite of the 
large observational errors which must be present 
in a region where 1A is equivalent to 40-50 wave 
numbers. 


EXPERIMENTAL 


These bands were photographed with a 21-foot 
focus normal incidence vacuum spectrograph, 


1A. B. F. Duncan, Phys. Rev. 47, 822 (1935). 
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a description of which has been published by one 
of us.? The grating used had 160,000 lines ruled 
30,000 to the inch on speculum. The dispersion in 
the first order of the region used was approxi- 
mately linear and averaged 1.328A/mm. The use 
of higher orders proved impossible, because of 
the fact that filters transparent to this region but 
opaque to the near ultraviolet and visible were 
unavailable. 

The Lyman continuum was produced by the 
discharge, once a second, of an 8-mf condenser 
charged to 3200 volts by a transformer and UV 
866 rectifier. The life of the 1-mm glass capillary 
(30 mm long) was about 90 minutes (5400 dis- 
charges). Exposures were 90-180 minutes‘on very 
sensitive Hilger-Schumann plates. The densities 
of the images were very far from being propor- 
tional to the exposure times, 270 minutes (3 
capillaries) giving little if any more blackening 


2G. R. Harrison, Rev. Sci. Inst. 2, 600 (1931); 4, 651 
(1933). 
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than 180 minutes. For this reason the contrast 
was very poor, especially at high pressures of am- 
monia necessary for the weak lines. Also for this 
reason visually estimated intensities have little 
meaning and are not given for these bands. 
Hydrogen, purified by diffusion through palla- 
dium, was pumped through the discharge tube at 
about 0.5 mm pressure. 

Alignment of the discharge tube with the 
spectrograph presented unusual difficulties in 
this case, and the use of a long absorption column 
between the source and slit greatly increased 
these difficulties. For this reason the spectro- 
graph itself was used as the absorbing column and 
the discharge was placed as close to the slit as 
possible. In some cases a window of lithium 
fluoride, transparent to 1000A was placed in 
front of the slit, to remove possible higher orders. 
This proved unnecessary, but the window served 
to prevent the spectrograph filling with hydrogen. 

The ammonia was purified by distillation over 
sodium. The liquid ammonia was stored over 
solid ammonium nitrate which reduced its vapor 
pressure to about 0.5 atmosphere. Due to the 
high speed evacuation of the spectrograph, a 
continuous flow of ammonia was undesirable, 
since a uniform pressure could not be obtained in 
this way. Consequently the valve to the pump 
was closed during exposures, and ammonia was 
admitted to the desired pressure. The ammonia 
disappeared slowly in the spectrograph, either 
by absorption or reaction, or more probably 
through the valve into the pump, which was de- 
signed to stand pressure differences in the oppo- 
site direction. The pressure was kept constant by 
admitting small quantities of ammonia at inter- 
vals of 2 or 3 minutes. The pressure was meas- 
ured by a McLeod gauge close to the slit, and for 
any one exposure probably represented a constant 
fraction of the true pressure in the main body of 
the spectrograph. The pressures of ammonia 
varied from 0.017 to 0.09 mm on different plates. 

Seven spectra were measured and every line 
was measured on at least two spectra. The wave 
numbers reported represent the unweighted 
mean of two or more observations. 

The plates were measured in the usual way. 
Emission lines in the source used as standards 
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were C IV 1548.214, 1550.774; N I 1745.246, 
1742.734.8 


RESULTS 


Lines in seven bands were measured, which 
are the 0,2—0,8 bands of the second electronic 
state.! The 0,0 and 0,1 bands did not appear at 
the highest pressure, and the 0,2 and 0,3 bands 
are very incomplete, lines near the head only 
appearing. The other five bands are much better 
developed and show two branches, with a broad 
central line. The heads in the R branches are not 
prominent at high dispersion, and are indeed 
difficult to locate at all. 

The absorption lines were rather broad at the 
lowest pressures, often as much as 0.1 mm. For 
this reason an actual resolving power better than 
3-5 cm is not claimed. This is about one-tenth 
the theoretical at 1500A. 

The lines observed appear in Table I. The lines 
of 41612 and 41587 could not be numbered but 
they appear to form the R branches of these 
bands. The lines of 41562, 1537, 1514, 1490, and 
1467 (wave-lengths of heads) are numbered. This 
numbering will be discussed later. 


EXPECTED STRUCTURE OF THE BANDS 


In its excited states ammonia is either a regu- 
lar pyramid or a plane. In the latter case we have 
the special relation 2A’=C’, while the general 
relation for a pyramid is A’—C’/2=yh?. Here C 
is the moment of inertia about the symmetry 
axis and A the moment around an axis in a plane 
through the center of mass perpendicular to the 
axis of symmetry, / the height of the pyramid 
and yw the reduced mass. The rotational energy 
levels are those of the symmetrical top in both 
cases. Our measurements are not sufficiently 
accurate to detect any but large deviations from 
the symmetrical top case, and these do not appear 
to be present. 

The rotational levels of the symmetrical top 
have been discussed in detail by Dennison, and 
by Dieke and Kistiakowsky.*:* The latter 


3 J.C. Boyce and C. A. Rieke, Phys. Rev. 47, 654 (1935). 

4D. M. Dennison, Rev. Mod. Phys. 3, 280 (1931). 

5G. H. Dieke and G. B. Kistiakowsky, Phys. Rev. 45, 4 
(1934). 
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authors have discussed particularly selection 
rules and intensities in perpendicular bands 
deviations from the symmetrical top, and have 
applied the theory to certain bands of formalde- 
hyde. The theory for parallel bands (to which 
type these ammonia bands belong) is much more 
simple and will be briefly stated here. 
The general form of the rotational energy is 





(I) K=0; A’=A”. 


(Il) K=0; A’#A”. a single P, Q and R branch. 


h? ¢J(J+1) 1 1 
rae MEY A 
Srl A c a 


The selection rules are AJ=0, +1; AK=0. The 
branches to be expected are gP, pQ, and gR.5 
The appearance of the bands will depend on the 
relative magnitudes of A’ and A”, C’ and C”, 
and on the absolute value of K. Several cases 
may be considered. 


a single P and R branch; a single line for the Q branch. 


(111) K#0; A’=A”; C’=C”. same as (I), since the separate P branches, Q branches and R branches coincide. 


(IV) K#0; A’#A"; C’#C”. 


The separations of the branches may have any 
values in case IV. The separation of the Q 
branches for the same value of J but two different 
values of K will be given by 


Av(cm~) 








41612 41587 A1514 41490 








R branch ? R branch ? J PW) RJ) PU) RV) 
62037 62964 0 65709, 66665 
61991 933 , 57 ae 
61962 908 }2 — — — 66732 
61932 898 | 3 641 782 610 
905 882 | 4 634 818 582 762 
878 861 5 609 829 557 794 
845 840 6 590 —_— - 800 
820 818 zz 871 522 829 
761 | 3 a . 837 
741 9 525 906 485 862 
ieirseecmcancegten ae ee 935 —- 895 
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J PW) RJ) PW) RJ)/|12 — 961 - 930 
0 63763 6477 }13 448 405 954 
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- iat — |15 402 . 
3 .710 63804 720 848/16 — 345 67000 
4‘ ag — 695 877|17 372 
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7 626 887 637 — |20 299 264 
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on e 22 270 231 
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im 551 24 174 
12 528 — 25217 
13 516 513 — 
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17 416 437 ii — — |10 — _ 67870 
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> - 386 ij. = 1 Oh 
21-352 365 15 583 — |14 939 
22 335 }¢ $37 — |1s 965 
23 7 Se « 16 989 
24 8 479 67832 
25 
26 





(general case) a separate P, Q and R branch for each value of K. 





In this general case, for each value of K we have a 
sub-band which is entirely analogous to a com- 
plete three-branch band of a diatomic molecule. 
These sub-bands would be separated from one 
another only if A’ and C’ had very different 
values, since A” is about 1/2 of C” (from infra- 
red measurements). 

Obviously, if C’ and C”, as well as A’ and C” 
are nearly but not exactly equal, case IV goes 
over into case I, with the result that the indi- 
vidual lines are more or less broad, because of the 
near coincidence of lines with the same J but 
different K values. This is equivalent to saying 
that the sub-bands almost coincide. 


DISCUSSION OF THE RESULTS 


We next see how far the observed structure of 
the bands is in agreement with the structure 
expected on the basis of the preceding theory. 
It is at once apparent that the observed structure 
is much more simple than case IV would predict. 
Calculations of the relative intensities of the 
sub-bands show that the intensities of sub-bands 
K=3n(n=0,1,2---) have an intensity only 
slightly higher than the sub-bands for which K 
has the value 3n+1, so that there seems no 
reason why all sub-bands (at least with low 
values of K) do not appear. Since only one sub- 
band is observed for each vibrational transition, 
we must conclude that case IV has practically 
degenerated into case I. This happens, as we 
have seen, when A’ is equal to A”, and C’ toC”’. 

® P. Lueg and K. Hedfeld, Zeits. f. Physik 75, 599 (1932), 
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There is additional evidence that A’=A”. 
The formula for each gP and gR branch, con- 
sidered separately, must depend on J, A’ and A” 
but not on C’, C” and K, and have the following 


form: 
gP= Vom 2B(J)+C(J*), 
qR=vp+2B(J+1)+C(J+1)’, 


where vp=the electronic+ vibrational term 


h 1 1 h 1 1 
: ( + ) and c= enol ). 
8xc\A’ A” 8r2c\A’ A” 


The branch lines have a linear rather than a 
quadratic dependence on J even at high J values, 
which points to the fact that C is very small. 
The average interval between lines should be 
then very nearly equal to 2B. Since A” is well 
known from infrared data, A’ may be determined 
from the above relation, and comes out very close 
to A”. 

If A‘=A” in a diatomic molecule, the head is 
formed at high values of J and the interval be- 
tween head and origin varies considerably with 
the vibrational quantum number.’ This is exactly 
what is found to be the case here. 

There is unfortunately no independent con- 
firmation that C’=C” but this assumption ap- 
pears to be necessary to explain the entire ab- 
sence of other sub-bands. 

All the bands studied come from the normal 
(nonvibrating) state of the molecule and com- 
bination differences can be formed which involve 
constants of this lower state only. These con- 
stants are well known from work in the infrared. 
For example, R(J—1)—P(J+1) =4B"(J+1/2), 
where B’’ =9.94 cm~. From this relation a table 
of differences for successive values of J may be 
constructed, and a search made for successive 
pairs of lines separated by the theoretical differ- 
ences. When this is done it is immediately evident 


2B= 


7™W. Jevons, Report on Band-Spectra of Diatomic Mole- 
cules (University Press, Cambridge, 1932). 
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that many lines are missing. This is not surpris- 
ing when the difficulty of finding weak lines js 
considered. For this reason the numbering of the 
lines had to be made in the following way. 

We may be certain of the origins of five of the 
seven bands. This gives, of course the positions of 
the Qo(=Po) lines. If we assume that A’=A” 
and 2B’=20, for simplicity of calculation, we 
can calculate the position of every line in these 
bands. Then every line found experimentally 
must coincide (within the limit of error) with a 
calculated line. This method is possible only 
when the intervals between lines are much greater 
than the experimental error. 

These calculations will not be reproduced here 
but the general results only given. A total of 109 
lines (exclusive of the Qo lines) were measured. 
From the number of blank spaces, this turns out 
to be a little more than 75 percent of all the lines 
expected. Of the measured lines 60 percent fall 
within 5 cm of a calculated line and 82 percent 
within 8 cm~. There is no trend of differences 
between observed and calculated even at high J 
values. This should certainly be observed if the 
difference between A’ and A” were larger than 
the experimental error. The missing lines appear 
to be scattered entirely at random in the bands 
and thus no theory of their absence can be given. 

From this discussion we may be reasonably 
sure that A’ is not far different from A’, and C’ 
is probably nearly equal to C”’. However A’ 
must be somewhat larger than A” since a head, 
however indefinite, is formed in the R branch 
If C’ were exactly equal to C”, the pyramid would 
be slightly higher in the excited state. A variation 
undetectable in our measurements could, how- 
ever, keep the height the same or make the 
pyramid slightly flatter. Much more accurate 
measurements would be necessary to determine 
unambiguously the dimensions in the excited 
states. 

One of us (A. B. F. D.) expresses his apprecia- 
tion to the Warren Fund of the American 
Academy of Arts and Sciences for financial aid 
during the progress of this work. 
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The Infrared Absorption Spectra of HCl in Benzene 


E. K. PLYLER AND DupDLEY WILLIAMS, Department of Physics, University of North Carolina 
(Received November 30, 1935) 


The infrared absorption of benzene and a solution of 
HCI in benzene has been measured in the region from 2.54 
to 4.2u. The absorption band of benzene at about 3.38, is 
increased in intensity and shifted to longer wave-lengths 
in the solution for the same cell thickness. Other benzene 
bands in the region from 24 to 3u did not show an appreci- 


HE spectra of hydrogen chloride in the 

gaseous state have been studied by a num- 
ber of investigators. Under the resolution ob- 
tained with the grating spectrometer there ap- 
pears a strong band at 3.46y consisting of the 
lines of the positive and negative branches, the 
zero branch being absent. The first harmonic of 
this band has also been observed at 1.76% and 
it is found to have the same general character- 
istics. In solid hydrogen chloride? the funda- 
mental band appears at 3.74u with a series of 
broad lines equally spaced about a zero line. In 
aqueous solution hydrogen chloride is highly 
ionized, whereas in substances like benzene 
ionization does not take place. In this investiga- 
tion the infrared absorption spectra of both types 
of solutions were studied to observe the effects of 
distortions of HCI molecules in the benzene and 
to detect the possible presence of undissociated 
HCI molecules in water. 

A large type Hilger infrared spectrometer with 
a fluorite prism was used as the resolving instru- 
ment. The effective slit width at 3.4u was 0.03,. 
Fluorite windows were used in preparing the 
cells, their thickness being 0.03 mm. The cells 
were made by placing mica washers of this thick- 
ness between the fluorite plates. 

In the case of aqueous solutions of HCI no 
absorption other than that due to water was 
observed in the region from 3 to 3.6u, the con- 
centration being 5 normal. When thicker cells 
were used, the 3u water band was greatly broad- 
ened for the HC! solution, but no separate ab- 


1 W. Burmeister, Vesh. d. D. Phys. Ges. 15, 589 (1913); 
E. v. Bahr, Vesh. d. D. Phys. Ges. 15, 1154 (1913); H. M. 
Randall and E. S. Imes, Phys. Rev. 15, 152 (1920); 
E. S. Imes, Astrophys. J. 50, 25 (1919). 

*P. E. Shearin, Phys. Rev. 48, 299 (1935). 


able change in intensity or position. It was concluded that 
the absorption in the region of 3.44 was due to the HCI in 
solution. When the ratio of the absorption of the solution 
to the absorption of the benzene is found, there appears a 
band with center at about 3.494. The relation of the 
absorption to the dipole moment is discussed. 


sorption band was observed. This change in the 
3u water band is probably due to hydration 
effects. However, there was no absorption band 
in this region which could be attributed to un- 
dissociated HCl molecules, showing that if 
molecules are present in aqueous solution their 
number is insufficient for the cell thickness used 
to produce an observable band. 

In the investigation of the benzene solution it 
was necessary to obtain water free materials. A 
sample of Baker’s chemically pure benzene was 















z 
° 
” 
2g 
= 
77) 
z 
a 
ig 
- 
= HCI in 
= 5 ANHYDROUS 
° BENZENE 
# sor 

30Fr 

10r 

iL i = = i i 1 i 





3.4 38 42 
WAVE - LENGTH Cp) 


Fic. 1. The transmission of anhydrous benzene and HCl in 
benzene, the cell thickness being 0.03 mm. 


215 











216 gE. &. 
BENZENE 

70 
60 

gro 

z 

2 HCI IN BENZENE 

« 

Lan 

SOF 

3 

rm L 
30h 











i 4. i. i 
22 24 26 28 
WAVE -LENGTH () 
Fic. 2. The transmission of HCl in benzene and benzene 
in the region from 2u to 3u, the cell thickness being 
0.6 mm, 


dried over fresh sodium chips for three days be- 
fore the solution was prepared. The hydrogen 
chloride gas was passed through three drying 
towers containing concentrated sulphuric acid 
before it was bubbled into the benzene. The 
benzene solution was saturated at 25°C. 

Since the strong absorption band of HCl in 
the gaseous state occurs at 3.46u, the solutions 
were studied over the range from 2.5u to 4.2y. 
The bands observed in this region are shown in 
Fig. 1. The upper curve represents the transmis- 
sion of anhydrous benzene for a cell thickness of 
0.03 mm. The second curve shows the transmis- 
sion of hydrogen chloride in the benzene solution 
for the same cell thickness. It is seen that the 
widening of the observed absorption band in the 
lower curve takes place entirely on the long wave- 
length side of the band. In order to ascertain if 
the presence of HCI in benzene produced any 
appreciable change in the other parts of the 
benzene spectrum, the absorption was studied in 
the region from 2y to 3u with a cell thickness of 
0.6 mm. The results are shown in Fig. 2. It is ob- 
served that the two transmission curves are very 
similar, indicating that no modification is pro- 
duced by the presence of HCI in the solution. 

The characteristic band in the sample of an- 
hydrous benzene used in this work appears at 
3.384. The position of this broad band has been 
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shown by Barnes’ to vary with cell thickness. He 
observed that the benzene absorption band con- 
sisted of eight maxima lying between 3.23u and 
3.51. 

In Fig. 3 is shown the ratio of the transmission 
of the HCI solution to that of pure benzene for 
various wave-lengths. This curve, which repre- 
sents the absorption of the hydrogen chloride, 
has its maximum at about 3.49u, this being in 
the same general region as the band for gaseous 
HCl. The position of the HCI absorption cannot 
be determined very accurately from this curve on 
account of the fact that benzene varies rapidly 
in its absorption in this region. This change in 
absorption of benzene in this region may produce 
a distortion and a shift in the band when calcu- 
lating the ratio of transmission. The band as 
observed in this work shows a slight shift to the 
longer wave-lengths as compared with the gase- 
ous state. This diminution of the vibrational 
frequency indicates a slight decrease in the bind- 
ing energy between the hydrogen and chlorine 
nuclei. This change may be connected with the 
increase of electric moment of HCI from 1.08 
< 10-"* e.s.u. in the gas to about 1.3 x 10-"* e.s.u. 
in benzene, as reported by Fairbrother.‘ It is in- 
teresting to note that this intense absorption 
obtained in the infrared spectra does not appear 


*R. B. Barnes, Phys. Rev. 36, 296 (1930). 
4 F. Fairbrother, J. Chem. Soc., Part I, 43 (1932). 
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BAND SYSTEM 


in the Raman spectra. Woodward? has reported 
his inability to observe the Raman spectrum of 
HC! dissolved in benzene, pointing out that the 
negative result was probably due to the low 
solubility of HCI in this solvent. In the case of 
certain other organic solvents,* West and Arthur 
have overcome this difficulty by studying the 
Raman spectra at a temperature of — 65°C. They 
were able to observe the Raman spectra of HCl 
in solution at this temperature for several sol- 
vents. Evidently the coefficient of absorption in 
the infrared of HCI in solution is greater than the 
coefficient for HCI in the gaseous state, as can be 
seen from the magnitude of the observed band as 
shown in Fig. 3. An approximate calculation 
reveals that the coefficient of absorption of HCl 
in solution is three to four times greater than in 
the gaseous state. On account of the variation in 


5L. A. Woodward, Physik. Zeits. 32, 777 (1931). 
6 W. West and P. Arthur, J. Chem. Phys. 2, 215 (1934). 
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concentration with temperature and errors in 
cell thickness the coefficient of absorption could 
not be determined very accurately. This indicates 
that the change in the electric moment causing 
absorption in solution is greater than the corre- 
sponding change in the gaseous state. These re- 
sults are in accord with the observation of Fair- 
brother for the increase of the total dipole 
moment of HCI when in benzene solution. This 
increase in change of dipole moment may explain 
the difficulty of observing the Raman spectra at 
room temperature. 

On account of intense bands of benzene occur- 
ring in the regions of 1.764 and 3.46uy it is difficult 
to observe the absorption bands of HC! in a 
benzene solution accurately. However, there are 
other organic solvents which do not have bands 
in these regions that may be used. Then the posi- 
tion of the HCI solution band can be determined 
with greater accuracy than in the present work. 


REVIEW VOLUME 4 


Note on a Band System of Caesium 


P. Kuscu anp F. W. Loomis, University of Illinois 
(Received December 9, 1935) 


A further attempt is made to analyze a band system of Cs: at 6250A. It has been possible to 
arrange a number of bands, degraded to the red, into a square array and to assign quantum 
numbers to the lower state. An assignment reported for this system by Matuyama is shown to 


be incorrect. 


HE molecule Cse gives rise to several band 

systems in absorption. An analysis has 
been given! for several of these systems, but an 
extremely complex system whose maximum in- 
tensity lies at about 6250A has not been satis- 
factorily analyzed. Rompe? has observed this 
system in absorption and Matuyama’ has re- 
ported an analysis on the basis of Rompe’s data. 
While the data of Rompe agree in a rough way 
with those of Loomis and Kusch, many close 
groups of bands are not resolved by him, nor 
does he make any reference to the blue degraded 


‘ Loomis and Kusch, Phys. Rev. 46, 292 (1934). 
* Rompe, Zeits. f. Physik 74, 175 (1932). 
use. Tohoku Imp. Univ. Sci. Rep. 23, 322 


bands which appear in this system and were re- 
ported by Loomis and Kusch. 

It has been possible to find several progressions 
among the bands degraded to the red and to 
assign the quantum numbers of the lower state 
to the bands in these progressions by reference to 
the accurately known ground state of caesium. 
These progressions have been fitted into a square 
array. As was the case with the previously re- 
ported! bands in this system, degraded to the 
blue, it has been impossible to follow the system 
to its origin and therefore to assign quantum 
numbers to the upper state. Matuyama reports 
an origin of the system at 15,801 cm~ but this 
assignment is believed to be incorrect. For, 
certain of Matuyama’s progressions may be iden- 


r. &2UY 


TABLE I. Band heads in the system at 6250A. 

















m cn | ; . es — differences constant within experimental error. 

| obs n —Vealc | obs n —Vealc e 

a ae s The bands, degraded to the red, which have 

2 16054.05 5 6 —0.39/ 5 15890.70 8 12 —0.14 been assigned into a square array are accurately 

2 16029.17 1 4 -—0.16| 7 15878.13 6 11 —0.23 , . 

1 16027.87 4 6 +0.03| 3 15877.33 9 13 +0.29 represented by the formula 

3 16025.26 7 8 -—044!| 9 1586449 7 12 —0.19 

2 1601543 2 § +0.24| 5 15851.60 5 11 —0.30 y=16,175.80+27.34(n’ +4) —0.0733(m’+}4)? 

4 16013.27 5 7 -—0.33|} 7 15850.82 8 13 —0.20 my. 51 (oy? 11)\2 

4 16011.81 8 9 -—0.46| 8 15837.47 9 12 +0.09 —[41.990(0’’ + >) — 0.080051 (0 +3) 

4 16001.78 0 4 -—0.35|10 15824.70 7 13 —0.15 — 0.000164266(0"+3)*], (1) 

2 16000.47 3 6 —0.61| 3 15824.05 10 15 +0.26 

4 1599903 6 8 -—0.35| 4 15811.85 5 12 —O@5 ow :. ' : a ’ differs 

5° 1508830 1 § 4015] 6 1881106 8 14 -030 ” 3S the true quantum number, and differs 

4 15986.82 4 7 -—0.18| 5 15797.92 9 15 +0.01 from the true v’ by some unknown constant. The 

3 15985.01 7 9 —0.18} 4 15785.21 7 14 +0.01 a ae ee > > — a oe 

1 15974.73 2 6 +0.55| 4 1578450 10 16 +0.02 ground case is identical to the one previously 

4 1597256 5 8 rer 9 15771.83 5 13 —0.25 reported.' The observed band heads, their in- 

8 15971.19 8 10 +0.18 8 15 —0.06 —— cine ae on = 

5 1596055 0 5 —040| 6 15758.96 6 14 +0.0g ‘ensities, assignments and the differences be 
3 7 +0.31| 7 15758.54 9 16 —0.04 tween the observed values and those given by 

6 1594698 1 6 —0.16|} 6 15745.44 7 15 —0.29 > oe a. 5 

7 15946.33 4 8 +0.01/ 10 17 +0.08 ©4. (1) are given in Table I. 

4 15945.01 7 10 +0.16] 8 1573242 8 16 —0.15 The bands in this system, degraded to the blue, 

9 15932.13 5 9 -—0.28} 7 1571945 9 17 -—0.03 .. a . ay he renresente : 

4 1593082 8 11 —0.02| 6 1570645 10 18 4.0.04 and previously reported may be represented by 

1 15919.24 3 8 -—0.32| 4 1569340 8 17 —0.06 the formula 

5 15918.18 6 10 —0.35} 6 15680.52 9 18 —0.01 = 29 38(y’ +1) — 7 ’41)\2 

0 15905.99 4 9 +0.19| 5 15667.79 10 19 40.15 ”~ 16,066.03-+2 pe +3) —0.07 O(n +3) 

4 15904.46 7 11 -—0.22| 4 15654.75 11 20 —0.03 —[41.990(v +4) —0.080051(v" +3)? 

4 15892.00 5 10 -—0.07| 4 15641.78 9 19 +0.02 fo ,) g//+1)3 (2 

| 4 1562931 10 20 +0.26 0.000164266(v' +3) }. 2) 


tified with ours, but the quantum numbers of the 
ground state which he assigns to the bands in the 
progressions differ from the true quantum num- 
bers as determined from the known ground state. 
Furthermore, a large number of the bands as- 
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signed by Matuyama cannot be assigned into 
such an array, while keeping the combination 


The formula previously reported' contained an 
error. 

This band system is extraordinarily complex 
and a great many observed features have not 
been explained, even though the two systems, 
degraded to the red and to the blue, do include 
the more obvious features of the system. 
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The magnetic rotation spectrum of Rbz has been observed and a vibrational analysis made. 
The vibrational frequency in the ' ground state is 57.31 cm™ and that in the excited ‘II state, 
48.05 cm~'. The magnetic rotation spectrum of NaRb has also been found and has the vibra- 
tional frequencies 106.6 cm™ and 61.5 cm™ in the 'S ground and excited 'II states, respectively: 
A band system due to RbCs has been identified. 


HE purpose of this investigation is toextend have been made! of the spectra of Lie, Nae and 





our knowledge of the band spectra of the 
alkali metals to Rbe, the only one of the group 
not previously investigated by means of accurate 
methods. Rather complete vibrational analyses 


Ke, and their energies of dissociation determined 
by means of their magnetic rotation spectra. 


1 Loomis and Nusbaum, Phys. Rev. (a) 38, 1447 (1931); 
(b) 39, 89 (1932); (c) 40, 380 (1932). 
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SPECTRA OF 


The absorption spectrum of Cs: has been in- 
vestigated? and an accurate value of its energy 
of dissociation found. Matuyama® has reported 
an analysis of the spectrum of Rbe under low 
dispersion but has not found a value of the 
energy of dissociation. 

In the present investigation the magnetic 
rotation spectrum and the absorption spectrum 
of Rbe in the region 6500-7100A have been 
observed and measured; the two spectra have 
been correlated and a vibrational analysis made, 
in the usual way. The magnetic rotation spec- 
trum extends further than the absorption spec- 
trum, and is more completely developed. This is 
in accord with the usual experience with mag- 
netic rotation spectra, which have, for this 
reason, been found to be ideal for extending a 
band system to the high quantum numbers 
necessary for the determination of an accurate 
energy of dissociation. The data are, however, 
in this case, still not extensive enough to yield a 
reliable value of the energy of dissociation. 

The same apparatus was used for the observa- 
tion of the magnetic rotation spectrum that 
Loomis and Nusbaum!* used for their investiga- 
tions of the other alkali metals. The nickel tube 
which they used was replaced by a glass tube to 
which windows had been fused. The ends of the 
tube were carefully annealed to relieve the 
strains introduced in the windows by the process 
of fusing the windows to the tube. The success 
of the method seems to depend on a careful 
annealing, for the presence of strains in the 
windows makes it impossible to obtain extinction 
when the tube is placed between crossed Nicols 
in the absence of a magnetic field. 

Rubidium, prepared by heating dry rubidium 
chloride with metallic calcium, in vacuum, was 
distilled into the system. Temperatures and 
vapor pressures corresponding to temperatures of 
about 320°C could be applied to the tube for 
long intervals of time without appreciably 
blackening the tube. At higher temperatures the 
inside of the tube quickly blackened and became 
opaque. The impossibility of working with higher 
pressures of rubidium with this technique 


* Loomis and Kusch, Phys. Rev. 46, 292 (1934). 


*Matuyama, Tohoku Imp. Univ. Sci. Rep. 23, 308 
(1934), 
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accounts for the difficulty of extending the band 
system to very high quantum numbers. 

Rubidium has two isotopes, Rb*®* and Rb*’ 
which occur in about the ratio 3 : 1. The three 
molecules Rb*Rb*, Rb®Rb*? and Rb*’Rb*’ 
should therefore have the relative concentrations 
9:6:1. The bands due to Rb*®Rb® and 
Rb*Rb*? may be expected to be of comparable 
intensity. Now in these magnetic rotation spectra 
the most prominent features are the strong lines 
just to the violet of each head, the rest of the 
spectrum being obscured by reabsorption. Hence 
presumably the heads measured belong to the 
molecule Rb*®Rb®* or Rb*Rb*’ according to 
which has the shorter wave-length. In the case 
of rubidium, then, .the observed lines in magnetic 
rotation should correspond to the band heads 
due to the molecule Rb**Rb* at the short wave- 
length side of the origin, and to the band heads 
due to Rb**Rb*? at the long wave-length side 
of the origin. 

Equations to represent the frequencies of the 
assigned bands have been found by a method of 
least squares. The equation 


v= 14,662.6+48.05(v' +4) —0.191(0' +4)? 
—[57.31(v"” +4) —0.105(v"+4)2] (1) 


has been found to represent adequately the fre- 
quencies of the bands attributed to the molecule 


Rb®*Rb*. The equation 


= 14,662.6+47.78(v' +4) —0.188(v’ +4)? 
—[56.98(0" +4) —0.103(v"+4)?] (2) 


then applies to bands attributed to the molecule 
Rb*Rb*’ and has been obtained from the first 
equation by multiplying the constants by suit- 
able powers of p (p=0.9943). The quantum 
numbers, the measured frequencies, estimated 
intensities and the differences between the ob- 
served values and those calculated from Eqs. (1) 
and (2) are given in Table I. Due to a slightly 
misplaced comparison spectrum the values given 
in Table I may be in error by a small constant 
amount, but these errors probably lie well within 
the error of measurement. A great many of the 
observed band heads have also been measured 
in absorption on a 21-ft. grating. The positions 
of these band heads agree essentially with those 
observed in magnetic rotation and are not 
recorded except in the few cases where a band 
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TABLE I. Band heads of rubidium observed in magnetic rotation. 








: 
Yobs | Mac. Yobs 











Mac. #¢ v Yobs | MaG. # v i v 
Rot. Ass. OBSERVED v’ y’’ —"cak Rot. ABs. OBSERVED v’ vo” —"calc | Ror. Ass. OBSERVED t t —"cak 
3 15358.6 8 149456 3 10 —1.2 5 5 14527.1 4 2 +0.5 
3 15341.0 2 18 —2.6 + 4 14936.9 0 6 —1.4 5 3 14515.6 » « —1.6 
3 15322.7 3 19 —5.0 10 5 14925.7 1 7 —0.8 3 14508.6 6 4 +0.8 
3 3 15305.5 4 20 —6.1 7 6 14913.9 2 8 —0.7 10 10 14489.3 3.60 +1.1 
3 4 15288.1 3 18 +1.2 4 14902.5 3 9 0.0| 3 3 14483.0 4 1 +3.4 
1 4 15268.9 4 19 —2.3 10 6 14892.0 0 5 —0.5|; 4 14464.8 
2 15262.6 7 23 +0.2; 10 5 14881.0 1 6 —0.2 5 14450.3 7 #4 —1.9 
a 4 15252.7 5 20 —2.7 3 14868.9 2 7 —-0.7| 2 14441.6 8 § —1.2 
2 15245.6 8 24 0.0 2 14853.5 | 8 10 14432.4 4 0 +0.3 
4 15232.0 4 18 +1.6 10 10 14846.8 0 4 +04; 8 4 14422.7 5 1 —0.9 
2 15238.2 6 21 —1.2 6 5 14835.3 1 5 -0.1 4 3 14414.2 6 2 —0.7 
4 3 15215.9 5 19 +0.9 5 3 14822.3 2 6 —2.0 2 14395.9 8 4 —1.0 
2 15207.0 8 23 —0.7) 4 14815.5 3 8 14386.3 9 5 —1.4 
4 $ 15200.9 6 20 +1.6 3 14808.0 6 4 14376.0 5 0 —0.1 
4 3 15190.6 4 17 +1.4 12 10 14800.9 0 3 +1.1 8 4 14367.1 6 1 —0.8 
4 3 15183.2 7 21 —0.3 4 4 14790.3 1 4 +1.0 5 3 14358.5 7 2 —0.8 
5 2 15176.0 5 18 +1.8; 5§ 14778.8 2 5 +0.3 3 14343.3 9 4 +1.5 
3 15168.3 8 22 +0.9; 4 14766.6 3 6 —1.0) 1 4 14320.6 6 O +02 
7 4 15160.8 6 19 +1.9 8 10 147545 0 2 +1.6 3 3 14312.2 7 1 —0.1 
6 3 15147.7 + 16 +0.1} 2 14746.7 1 3 +4.0 4 3 14303.1 =. —0.9 
6 3 15134.3 5 17 +13; 5 7 14732.7 2 4 +03; 3 2 14294.4 9 3 —1.1 
8 3 15119.9 3 14 0.0 5 14727.5 1 14285.2 19 4 —1.7 
8 3 15107.2 4 15 +1.5 6 4 147206 3 5 —1.2 1 14256.7 8 1 —0.3 
8 3 15093.4 5 16 +2.0 8 10 14706.5 0 1 +0.8 4 14248.5 9 2 —0.4 
9 3 15078.1 3 13 +0.9 | 7 14697.7 1 2 +1.9 5 14239.3 10 3 —1.3 
8 3 15064.2 4 14 +9.8 7 5 14685.3 2 3 —0.5 5 14231.4 11 4 —0.7 
7 3 15050.3 5 15 +0.8 3 14676.0 3 4 +0.3 4 14223.4 12 5 —0.1 
5 15038.7 6 16 +3.4| 2 14665.6 4 5 +0.3 1 142140 13 6 —0.7 
6 3 15032.8 3 12 —1.3 4 4 14660.2 0 0 +2.3 2 14193.5 10 2 —0.5 
7 4 15021.9 4 13 +1.2 9 10 14649.7 1 1 +1.0 3 14184.8 i1 3 —1.0 
7 15016.3 1 9 +0.2 3 14629.2 3 3 —0.1 3 14176.4 12 4 —1.2 
7 15008.4 3 14 +12} 4 5 14618.7 4 4 —0.8 3 14168.4 13. 5 —0.8 
7 3 15002.6 2 10 —0.9| 6 8 14603.4 1 O +2.2 2 14160.0 14 6 —0.6 
8 4 14989.8 3 11 —0.8 5 14593.0 2 1 +0.9 | 1 14153.4 15 7 +1.4 
4 14983.8 0 7 +0.2 7 5 14582.8 3 2 +0.1 1 141156 14 5 +0.5 
6 14977.8 4 12 +0.2 | 1 14551.7 6 5 —2.0) 1 14108.4 15 6 +1.6 
7 14970.3 1 8 —1.2 | 8 10 14545.9 2 0 +41.3 1 14099.8 16 7 +1.5 
3 14964.3 5 13 —0.2;); 3 14538.9 3 1 +3.2 | | 14092.0 17 8 +2.4 
s @ 14958.7 2 5 | 


v=) 
| 

| 
| so 
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°P3/2 state. The ground state dissociates into two 


head was observed in absorption only. Intensities 
normal atoms. If the product of dissociation of 


are given both for magnetic rotation and ab- 
sorption, and where an intensity is given for the upper state is a ?P\/2 atom the energy of 
absorption only, the frequency is that observed dissociation of the lower state becomes 5100 cm™, 
in absorption. and if it is a *P 3/2 atom the energy of dissociation 

A Franck-Condon diagram for this system is becomes 4870 cm. Direct extrapolation to 
given in Fig. 1. Since the system studied lies on dissociation of the lower state yields 7820 cm™ 
the violet side of the resonance doublet of the Rb as its energy of dissociation. Since the upper 
atom, and since the magnetic rotation lines state could be followed nearer to dissociation 
correlate immediately with the absorption heads, extrapolation of the upper state probably yields 
one may conclude, by analogy with Liz and Naz, a better value. The energy of dissociation of 





- 





that the transition is 'IIl'®. 

Extrapolation to dissociation of the upper 
state of the system according to Eq. (1) yields 
17,685 cm~ as the energy of dissociation of the 
upper state. The upper state dissociates into a 
normal 2S atom and one in the ?P state, but it is 
not known whether this atom is in the *P4;2 or 


rubidium has been calculated’ to be 0.47 volt 
(3800 cm") by interpolation between the known 
energies of dissociation of potassium and caesium. 
Experience with the other alkali metals has 
shown that the energy of dissociation obtained 
by direct extrapolation is usually too high, so 
that the value, 5100 cm~, is probably an upper 
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Fic. 1. Franck-Condon diagram for the 'IIl—'Z system of Rb». 


limit to the energy of dissociation of Rbe. Since 
the energies of dissociation of Kz and Cs» are 
known to a high degree of accuracy the value, 
0.47 volt, obtained by interpolation between 
the known energies of dissociation of Kz and Cs 
is still the most reliable value of the energy of 
dissociation of Rb». 

Due to the method of preparing rubidium 
which was used, sodium was present as a rather 
important impurity. Consequently we attribute 
to NaRb the new band system which occurred 
on our plates in the region 5990-6336A. It is 
true that three bands of this system at 16,513 
em~, 16,571 cm~ and 16,684 cm~' were reported 
by Walter and Barratt‘ and attributed to Rb, 
but the vibrational intervals in this system con- 
~ ‘Walter and Barratt, Proc. Roy. Soc. A119, 257 (1928). 


firm our conclusion, for firstly, the ground 
interval, 107 cm~', agrees well enough with one 
which is prominent in the green system which 
Walter and Barratt do attribute to NaRb; 
and secondly the vibrational intervals are ap- 
proximately those which are obtained by interpo- 
lation between the known intervals of Naz and 
Rbe. Moreover since the new band system was 
observed in magnetic rotation it is due to a 
'TIl'> transition. By comparing the vibrational 
frequencies of 107 cm~ and 61 cm™ in the 'Z 
and 'II states of this molecule with the corre- 
sponding vibrational frequencies of the alkali 
metal molecules given in Table II, it is evident 
that the new band system is not due to any of 
the alkali metal molecules composed of similar 
atoms. All of the possible intermetallic com- 
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; T tem) TABLE III. Band heads of NaRb observed in magnetic 
al eo 200 o 400 500 20 rotation. 
° ° ° — ——————————o——————* 
v vobs v vobs 
cs b © ° ° i OBSERVED 1’ —v¥ealc # OBSERVED wv ov ~vcalc 
lem") e a a 3 166783 5 O +0.1/5 16248.1 > = 0.0 
egoo|. 3 166254 4 0 -08/2 16189.7 0 2 +412 
6 “ 3 16572.2 3 0 0.0/4 161435 1 3 -07 
3 165174 2 0 +1.0/4 160831 0 3 —-15 
um © 3 164669 3 1 +04/4 159823 0 4 +4046 
y 5 164106 2 1 —0.1/3 15879.7 0 5 +02 
. 5 16351.7 1 1 1.3/3 157781 0 6 -03 
2 16305.8 2 2 0.0 | 
16800 


Fic. 2. Franck-Condon diagram for the 'IIl<!'Y system 
of NaRb. 


pounds may be excluded which do not have one 
atom of rubidium. Then since the vibrational 
frequencies of a molecule composed of a rubidium 
atom and an atom of some other alkali, may be 
expected to lie between the vibrational fre- 
quencies of Rbz and those of the molecule com- 
posed of two atoms of the other alkali, the new 
band system could be due only to NaRb and 
LiRb. But since Li was not observed .as an 
impurity in any of the observations, the molecule 
is almost certainly NaRb. 
The formula 


9 


v= 16,421.8+61.49(v' +4) —0.945(v’ + 4)? 
—[106.64(0" +4) ~0.455(0"+3)*] (3) 


was found to represent accurately the frequencies 
of the observed band heads. The observed fre- 
quencies, intensities, assigned quantum numbers 
and the differences between the observed fre- 
quencies and those calculated from Eq. (3) are 
given in Table III. A Franck-Condon diagram 
is given in Fig. 2 and is seen to be normal. 

The system between 7230—-7400A which was 
found by Loomis and Kusch? in caesium absorp- 


TABLE II. Vibrational intervals in the 'Z and 'Il states for the 
alkali metals. 








4 


wi 





Li, 351.6 269.7 
Nae 159.2 123.8 
Ke 92.6 75.0 
Rb» 57.3 48.1 
Cs, 42.0 34.2 














rubidium absorption 
found by 


tion also occurs on our 
plates. This system has 
Matuyama’ and attributed to Cse. However, as 
Loomis and Kusch pointed out, the vibrational 
intervals do not correspond to any known 
intervals in Cs». Neither do they correspond to 
any known intervals in Rb». But the observed 
intervals of 49.4 cm™ and 38.3 cm™ do fit nicely 
with intervals for RbCs estimated by interpola- 
tion between known intervals for Rbs and Csy». 
Since Cs is known to be present as an impurity 
in the Rb, and vice versa, it seems safe to attribute 
this system to RbCs. Although the progressions 
in this system are not long enough to fix the 
vibrational assignments uniquely, nevertheless 
the assignments given by Loomis and Kusch are 
very plausible, since, in addition to yielding 
vibrational intervals, which, on the above view 
are entirely reasonable, they also lead to a normal 
Franck-Condon diagram, which Matuyama’s 
assignments do not. The arbitrary numbers 
previously assigned to the bands are undoubtedly 
the true quantum numbers v’ and v’’. 

The band heads in this system are then given 
by the formula 


v= 13,747.21+38.46(0' +3) —49.41(0'"+ 3). (4) 


also been 


The data given by Matuyama are more extensive 
than those given by Loomis and Kusch. 

The author wishes to express his appreciation 
to Professor F. W. Loomis for his interest and 
helpful suggestions. 


> Matuyama, Tohoku, Imp. Univ. Sci. Rep. 23, 322 
(1934). 
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The Disintegration of Nitrogen by Neutrons 


T. W. Bonner* ano W. M. Brusaker, Kellogg Radiation Laboratory, @alifornia Institute of Technology 


(Received November 25, 1935)” 


Nitrogen disintegrations produced by neutrons have 
with nitrogen at reduced pressures. All of the obser 


the following reactions: 


A 
A 







en observed in a cloud chamber filled 
d disintegrations are attributed to one of 


7N'4+ on'—,;B"+.He!, 1N'* + on!—7O* + iH}, 7N'4+ on'—,Het +.Het+ Li’. 


The first of these reactions is found to go for “‘slow neutrons” as well as fast ones. The value of 
the energy of disintegration for this reaction is 2.33+0.26 MEV. For fast neutrons, the results 
are similar to those of Kurie, with nearly a constant amount of energy (3.1 MEV) appearing 


in the disintegration forks. The second and third reactions occur less frequently than the 


first one. 


HE disintegration of nitrogen by fast neu- 
trons was first observed by Feather.' He 
obtained photographs of forked tracks in a 
cloud chamber which was traversed by neutrons, 
and interpreted the two branches of the forks as 
tracks of the two disintegration products. 
Further experiments of this kind have been 
done by Harkins, Gans and Newson,? Meitner 
and Philipp* and by Kurie.* 
Most of the nitrogen disintegrations were 
attributed to the reaction 


>N4+ on'—;B"+.He* (1) 


but some forks were observed where one of the 
disintegration particles seemed to be singly 
charged. The two possible reactions proposed to 
explain this type of fork were: 


7N“+ on'—,C®+,H?, (2) 
>N“4+ on'—,C"+,H" (3) 


In the experiments described below we have 
shown that reaction (1) takes place for slow 
neutrons as well as for fast ones. A preliminary 
report of this work has been given in a Letter 
to the Editor.’ We have also found that reaction 
(3) is responsible for some of the observed dis- 
integration forks. One trident was observed, 
which we have attributed to the reaction 


7N'“+ on!—.He!+.He!+3Li’. (4) 


* National Research Fellow. 
(1933) Feather, Proc. Roy. Soc. A136, 709 (1932); 142, 689 
53). 
* Harkins, Gans and Newson, Phys. Rev. 44, 945 (1933). 
* Meitner and Philipp, Naturwiss. 20, 929 (1932); Zeits. 
. Physik 87, 484 (1933). 
*F.N.D. Kurie, Phys. Rev. 45, 904 (1934) ; 47, 97 (1935). 
* Bonner and Brubaker, Phys. Rey, 48, 469 (1935). 
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EXPERIMENTAL PROCEDURE 


In these experiments the Pyrex glass ring of 
the cloud chamber was replaced by a thin brass 
one, as the boron in the Pyrex absorbs slow 
neutrons strongly. The chamber was filled with 
nitrogen and saturated water vapor at reduced 
pressure. The calculated stopping power of the 
gas in the chamber (0.60) was checked by 
measuring the ranges of polonium alpha-par- 
ticles. By comparing their mean lengths in the 
chamber to their known mean range of 3.805 cm 
at standard conditions, we experimentally de- 
termined the stopping power to be 0.592. Three 
stereoscopic images of the tracks were pho- 
tographed through an f 3.5 lens; illumination was 
provided by a 2000-watt movie flood lamp. 
Examples of nitrogen disintegrations photo- 
graphs are given in Figs. 1 and 2. 


DISINTEGRATIONS BY SLOW NEUTRONS 


From Bethe’s® values of the masses, reaction 
(1) appears to be exothermic by 1.5 MEV, and 
hence should be expected to take place for slow 
neutrons in a manner similar to that of lithium 
and boron. In fact, Chadwick and Goldhaber’ 
reported a small but definite increase in the 
number of counts from a nitrogen filled ioniza- 
tion chamber when the source of neutrons and 
the chamber were surrounded by paraffin. 

If the disintegration is produced by a neutron 
with nearly zero energy, the ;B" and Het 
particles recoil in opposite directions and appear 


© H. Bethe, Phys. Rev. 47, 633 (1935). 
7 Chadwick and Goldhaber, Nature 135, 65 (1935). 


223 





224 T. W. BONNER AND W. M. BRUBAKER 
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D 


Fic. 1. A, B—The short straight tracks on both pictures are those of nitrogen disintegrations by slow neutrons 
C, D—The disintegration fork on C is the ordinary type due to the reaction N'+n!—B"-+ He‘. The fork on D is 


due tothe reaction N4“+n'—~C'+ Hl}, 


as a single straight track in a cloud chamber. 
These tracks should all have the same length 
(neglecting straggling) as the energy contributed 
by the slow neutrons is negligible. 

In our experiments we have used three dif- 
ferent sources of neutrons: Be+H?+paraffin, 
Be+H?, and Li+H?’. In the first experiment a 
paraffin sphere with a radius of 7.5 cm sur- 
rounded the source and additional paraffin was 
placed around the chamber, so that a consider- 


able number of slow neutrons were present. 
We have taken 9000 stereoscopic photographs of 
the chamber and on these have observed many 
straight tracks beginning and ending in the 
chamber, as well as several hundred forked 
nitrogen disintegration tracks. The range dis- 
tribution of the straight tracks for the different 
neutron sources is given in Fig. 3. All tracks 
were measured regardless of their orientation in 
the chamber. Three different types of tracks are 
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DISINTEGRATION OF NITROGEN 





Fic. 2. A } 
One of the particles must have been scattered in order for ¢ to appear so small. B 
of nitrogen into twoalpha-particlesand a Li’ atom. This disintegration was produced by a 12.9 MEV 


neutron. 


represented: those due to recoil protons, recoil 
nitrogen nuclei, and those due to the disintegra- 
tion of nitrogen. The recoil protons result from 
the hydrogen in the water vapor. For all of the 
sources, the recoil protons will be responsible for 
a small number of tracks of all lengths in the in- 
terval measured. These will appear as a general 
background. The nitrogen recoils will be quite 
short, and their maximum range will depend on 
the maximum energy of the neutrons present. 


to 
i) 
w A) 


\ stereoscopic picture of a nitrogen disintegration where the measured @ was 51 


A disintegration 


The neutrons from beryllium have a maximum 
energy of 4.42 MEV,’ and should project nitrogen 
nuclei with a maximum range of 0.28 cm, which 
is in good agreement with the experiments. The 
neutrons from lithium have a maximum energy 
of 13.3 MEV,° and should project nitrogen nuclei 
with a maximum range of 0.51 cm. As only a 
small portion of the neutrons from lithium have 


* Bonner and Brubaker, Phys. Rev. 47, 910 (1935). 
’ Bonner and Brubaker, Phys. Rev. 48, 742 (1935). 
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Fic. 3. Range distribution of the tracks in the nitrogen 
filled cloud chamber. 


energies near the maximum, one can only esti- 
mate the maximum range of the recoil nitrogen 
nuclei ; this estimate from the curve is consistent 
with the above-mentioned 0.51 cm. 

The group of tracks with a mean range of 
1.06+0.02 cm (15°C, 76 cm Hg) is due to the 
combined tracks of the alpha-particles and ;B" 
nuclei produced in the disintegration of nitrogen 
by slow neutrons. The curves show that very few 
tracks of this range were observed except when 
the neutrons were slowed down in the paraffin, 
indicating that very few (if any) ‘slow neutrons”’ 
are emitted from either of the sources.'® 

Table I gives the mean ranges, velocities, and 
energies of the ;B"™ and .He?* particles. We have 
used Bethe’s range-velocity curve for alpha- 
particles." 


TABLE I. 
ENERGY 
NUCLEUS RANGE VELocity X 107° IN MEV 
sBu 0.22 0.329 0.62 
2Het 0.84 0.907 1.71 
1.06 2.33 


The range velocity curve for B" has been com- 
puted from the nitrogen range velocity data” by 
means of Blackett’s relation 


10This agrees with our previous results with indium 
(Phys. Rev. 48, 470 (1935)) which indicated that there 
were no slow neutrons emitted by a Li+H? source. 

11 We wish to thank Professor Bethe for communication 
of the values used in Table I which are based on his new 
range velocity curve for low energy alpha-particles. 

2 P, M.S. Blackett, Proc. Roy. Soc. A135, 132 (1932). 
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R=k(M/Z')f(o). 


The probable experimental error of 0.02 cm 
in the determination of the range would intro- 
duce an error of 0.06 MEV in the energy. The 
probable error which is introduced by the in- 
accuracy of the range velocity curves for »He* 
and ;B" particles of this energy is difficult to 
estimate, but perhaps is about 0.2 MEV 

Using the disintegration energy found in this 
experiment and Bethe’s values for the masses of 
the other atoms appearing in reaction (1) we 
calculate the mass of ;N™ to be 14.0085. 

Nitrogen disintegrations of this type (with 
slow neutrons) are no doubt largely responsible 
for the nitrogen cross section of 11.310-* 
sq. cm for slow neutrons as observed by Dunning, 
Pegram, Fink and Mitchell.“ This is approxi- 
mately 100 times as large as the cross section for 
disintegration by fast neutrons. 


DISINTEGRATION BY FAST NEUTRONS 


Several hundred disintegration forks similar to 
those shown in Fig. 1C were photographed on 
the films exposed in the experiments described 
above. From this number 91 forks could be 
accurately measured both in regard to the 
lengths of the tines and to the angle between 
them. Table II gives the results of such measure- 
ments. The measured lengths of the tracks were 
1.69 times as long as the range because the 
stopping power was 0.592. 

The proper disintegration energy, (Eu.e+£p 
—E,/15), which we shall call E, is the dis- 
integration energy which an observer at rest 
relative to the center of mass would measure. 
Fig. 4 gives the distribution of E. The maximum 
at 3.1 MEV is in agreement with the results of 
Kurie.* We have shown that the Q for slow 
neutrons is 2.06 MEV ;“ hence the most probable 
energy which the neutrons contribute is approxi- 
mately 1 MEV. From the curve it is evident 
that very few neutrons contribute more than 
2 MEV. However, from the calculated energies 
of the neutrons producing the disintegrations, 
we find that most of the disintegrations were 


13 Dunning, Pegram, Fink and Mitchell, Phys. Rev. 48, 
265 (1935). 

44 We obtained this lower value of Q by using the range- 
velocity curve for alpha-particles which was used in com- 
puting alpha-particle energies in Table II. 
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DISINTEGRATION OF NITROGEN 


TABLE II. Data on nitrogen disintegration forks. Rye=range alpha-particles; Rg=range ;B" particles; Eye=energy of 
alpha-particles in MEV;** E,=energy of 5B" particles in MEV; ¢=angle between two tines; EZ, =calculated energy of 
neutron in MEV. 




















: : d Ete EXe+Es ‘ : Ege Ene+Es 
Rue Rp Exe Ep ry En +Ep —E,n/15 RHe Rp Ene Ep ry En +EBp —E£,/15 
NEUTRON Source: Be +:H? NEUTRON Source: Be +:H? 
0.53 0.35 0.86 1.52 173 4.9 2.38 2.08 1.54 0.30 2.75 1.15 143 4.8 3.90 3.60 
0.95 0.36 1.72 1.56 155 4.4 3.28 2.99 2.12 0.25 3.57 0.83 141 5.7 4.40 4.10 
0.83 0.19 1.47 0.50 51* 18.4* 1.97 a 1.08 0.23 1.96 0.68 148 2.4 2.64 2.48 
0.94 0.37 1.70 1.66 155 4.8 3.36 3. 0.84 0.41 1.49 1.98 161 6.2 3.47 3.17 
1.16 0.13 2.11 0.26 146 3.1 2.37 2.07 | 1.46 0.38 2.63 1.77 146 6.1 4.40 4.10 
0.81 0.30 1.43 1.11 153 3.0 2.54 2.34 0.87 0.44 1.55 2.24 163 7.4 3.79 3.49 
1.47 0.20 2.64 0.53 149 2.9 3.17 2.98 1.14 0.46 2.07 2.48 165 6.6 4.55 4.25 
0.82 0.39 1.45 1.82 154 6.4 3.27 2.97 1.96 0.23 3.36 0.71 144 4.7 4.07 3.77 
0.91 0.30 1.64 1.15 148 3.6 2.79 2.55 1.11 0.43 2.02 2.19 158 6.3 4.21 3.91 
1.81 0.18 3.14 0.48 150 3.8 3.62 3.37 . . 
0.93 0.22 1.68 0.65 148 2.1 2.33 2.19 NeEuTRON Source: Li +:H? 
1.08 0.12 1.96 0.22 143 3.3 2.18 1.96 1.42 0.40 2.56 1.87 133 11.3 4.43 3.67 
1.78 0.17 3.10 0.39 141 5.3 3.49 3.19 0.83 0.37 1.47 1.61 162 4.4 3.08 2.79 
1.25 0.23 2.28 0.71 146 2.9 2.99 2.80 0.43 0.38 0.66 1.77 156 8.7 2.43 1.85 
1.48 0.34 2.65 1.44 142 5.9 4.09 3.79 1.13 0.27 2.06 0.95 144 3.6 3.01 2.77 
0.53 0.29 0.86 1.07 133 6.5 1.93 1.63 1.19 0.20 2.16 0.53 128 5.8 2.69 2.31 
1.35 0.13 2.45 0.26 145 4.0 2.71 2.44 1.01 0.34 1.84 1.44 154 3.7 3.28 3.04 
1.23 0.33 2.24 1.36 148 4.1 3.60 3.32 1,19 0.16 2.16 0.36 148 2.7 2.52 2.34 
2.12 0.24 3.61 0.77 142 5.5 4.38 4.08 1.17 0.50 2.13 2.87 136 16.5 5.00 4.10 
2.22 0.23 3.72 0.68 135 7.5 4.40 4.10 | 1.44 0.37 2.60 1.66 157 3.3 4.26 4.04 
0.96 0.13 1.74 0.26 156 1.7 2.00 1.89 1.28 0.49 2.33 2.74 158 8.7 5.07 4.49 
1.11 0.30 2.02 1.15 144 4.2 3.17 2.89 0.98 0.57 1.78 3.45 136 21.8 5.23 4.33 
2.13 0.24 3.61 0.74 146 4.6 4.35 4.05 2.25 0.23 3.75 0.68 137 7.1 4.43 3.95 
0.90 0.33 1.62 1.36 150 4.2 2.98 2.70 1.00 0.43 1.82 2.16 136 12.3 3.98 3.16 
1.07 0.38 1.95 1.77 152 4.4 3.72 3.43 0.63 0.47 1.06 2.61 130 18.5 3.67 2.77 
1.83 0.14 3.19 0.28 175 3.3 3.37 3.15 1.05 0.17 1.91 0.42 156 1.5 2.33 2.23 
0.90 0.18 1.62 0.45 160 0.8 2.07 2.02 1.05 0.29 1.91 1.07 140 ., 48 2.98 2.66 
1.78 0.17 3.10 0.42 140 5.5 3.52 3.22 0.87 0.35 1.55 1.48 157 4.1 3.03 2.75 
1.26 0.21 2.30 0.62 145 3.0 2.92 2.72 0.69 0.40 1.18 1.87 141 10.3 3.05 2.36 
1.54 0.21 2.75 0.62 143 3.9 3.37 3.11 0.€3 0.28 1.06 1.03 147 3.7 2.09 1.85 
2.48 0.14 4.01 0.30 143 7.7 4.31 4.01 1.31 0.56 2.39 3.40 172 9.5 5.79 5.16 
141 0.40 2.55 1.93 150 5.8 4.48 4.18 0.77 0.43 1.35 2.19 168 7.2 3.54 3.06 
0.83 0.46 1.47 2.48 162 8.9 3.95 3.65 1.57 0.24 2.80 0.74 136 5.6 3.54 3.16 
1.64 0.23 2.90 0.68 130 7.1 3.58 3.28 1.41 0.17 2.55 0.39 141 4.1 2.94 2.66 
0.95 0.40 1.72 1.93 149 7.2 3.65 3.35 0.94 0.12 1.70 0.22 137 3.3 1,92 1.72 
1.41 0.29 2.55 1.07 138 5.6 3.62 3.32 2.13 0.34 3.58 1.44 147 4.9 5.02 4.70 
1.12 0.30 2.04 1.11 130 7.6 3.15 2.85 0.83 0.38 1.47 1.77 139 8.9 3.24 2.65 
2.17 0.14 3.66 0.28 146 6.6 3.94 3.64 0.59 0.30 0.98 1.11 142 5.2 2.09 1.75 
0.83 0.43 1.47 2.14 175 6.3 3.61 3.31 1.53 0.64 2.74 4.13 168 13.0 6.87 6.00 
0.68 0.17 1.16 0.39 130 3.3 1.55 1.33 0.94 0.12 1.70 0.22 143 2.7 1,92 1.74 
1.19 0.17 2.16 0.39 173 0.8 2.55 2.50 0.81 0.40 1.43 1.93 175 4.7 3.36 3.05 
1.00 0.36 1.83 1.56 142 6.8 3.38 3.08 0 84 0.40 1.49 1.89 149 7.7 3.38 2.86 
1.64 0.12 2.90 0.24 140 5.9 3.14 2.84 | 0.41 0.42 0.62 2.09 163 il.1 2.71 1.97 
1.58 0.23 2.81 0.71 148 3.2 3.52 3.30 | 0.71 0.59 1.22 3.68 156 19.8 4.90 4.00 
1,17 0.35 2.13 1.48 145 5.6 3.61 3.31 | 2.55 0.17 4.10 0.39 167 4.4 4.49 4.20 
1.60 0.09 2.84 0.13 138 6.5 2.97 2.67 0.46 0.33 0.72 1.36 148 6.4 2.08 1.66 
' 
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* One of the particles must have been scattered. 
* We used the Cavendish Laboratory range-velocity curve for alpha-particles in these computations since Bethe’s new curve which was used in 
the last section was not available. 


caused by neutrons with energies above 2 MEV." will include errors due to three causes: (1) 
Thus we find that some kinetic energy dis-. experimental errors in measurement of track 
appears. This is in agreement with the work of lengths; (2) errors due to straggling in ranges of 
all the others'~* who have studied this reaction. tracks with equal energies and (3) errors in the 
To account for this lost energy, Kurie has pro- range-velocity curves for the particles. All three 
posed the emission of a gamma-ray in aradiative of these possible errors are more important for 
capture of the neutron. Feather has suggested 
that this energy goes into the excitation of the . 
5B" nucleus. a | 

Before trying to decide this question let us 7 
first consider how accurately we can calculate 
both the energy appearing in the forks and the 7 \ 
energy of the neutron producing the disintegra- 
tion. The energies of the ;B" and He‘ particles Vi 

° 42 24 36 = eo 2 

% This is consistent with an average neutron energy of PROPER DISINTEGRATION ENERGY 
3.5 MEV for Be+H? and 3.9 MEV for Li+H?, and a cross 
section for disintegration which varies only slowly with Fic. 4. The distribution of the proper energies (Eye+E yx 
energy. —E,/15) in MEV of the nitrogen disintegration forks. 
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the shorter range ;B" tracks than for the .He* 
tracks, and so when the amount of energy 
appearing in the ;B" particles is relatively large, 
the uncertainty in the proper energy is also 
large. For example, in the case where Ey.=3.72 
MEV, Es =0.68 MEV and £y.+ £3 =4.40 MEV, 
the probable error is only about 0.2 or 0.3 MEV, 
but in the case where Eq. =1.78 MEV, Ep =3.45 
and Exy.+£g=5.23 MEV, the probable error is 
of the order of 1 MEV. 

The calculated energies of the neutrons are 
affected by the errors listed above as well as 
any error in the determination of ¢. Errors in the 
determination of ¢ are due both to experimental 
errors in measuring the angle between the tines 
and to errors introduced when one of the particles 
is scattered through a small angle at a point so 
near the point of disintegration that it is un- 
observable. An example of this sort is given in 
the fork of Fig. 2A, where the angle ¢@ was 
measured to be 51°, but from momentum con- 
siderations it must have been at least 134°. 
The neutron energy calculated from this fork 
was 18.6 MEV, which obviously was in error as 
there were no neutrons present with energies 
over 4.5 MEV. A check on the accuracy with 
which neutron energies can be calculated is 
obtained by knowing that the maximum energy 
of the neutrons from beryllium and lithium are 
4.42 and 13.3 MEV, respectively. Quite a number 
of the calculated energies of the neutrons from 
beryllium were 2 or 3 MEV too high. This shows 
that an accurate determination of the neutron 
energy cannot be made from such disintegration 
data. However, it can be determined well enough 
to show that, with the lithium source, a con- 
siderable fraction of the disintegrations are pro- 
duced by neutrons with energy above 8 MEV. 
In these cases, the disappearance of kinetic 
energy cannot be questioned. 

An explanation of the distribution of E can be 
made by assuming that the excess energy which 
is carried into the reaction by the neutron goes 
into the excitation of the B"™ nucleus. Cockcroft 
and Walton" have shown that B" has excitation 
levels at 2.2, 4.5, and 6.8 MEV. These with 
perhaps a higher level, can account for the ob- 
served distribution if we assume that the boron 
nucleus is excited to the highest state possible, 


#6 J. D. Cockcroft, London Conference Report, 1934. 


BONNER AND 


W. M. BRUBAKER 

consistent with the alpha-particle having enough 
energy (>1.5 MEV) to get through the barrier 
of the N*. This interpretation is clearly very 
artificial and unconvincing. 

An alternative explanation of the distribution 
of proper energies has been proposed by Kurie, 
according to which the neutron is first captured 
by a radiative transition to a virtual level of N™. 
This level corresponds to an energy greater by 
1 MEV than that of N“+ !' and greater by 
about 3.1 MEV than B'"+Het; the N® is not 
stable and disintegrates into B''+He*. That part 
of the breadth of the energy distribution which 
is not conditioned by experimental errors in the 
determination of EZ, may be reasonably attrib- 
uted to the short life of the N". This explanation 
is quite consistent with our results, and has 
only theoretical arguments against it. No reason- 
able model can be constructed for which the 
radiation emitted from the recoiling N" is of 
the right order of magnitude to explain the ob- 
served cross sections for the process; on this 
theoretical basis a cross section of 1/1000 of that 
observed would have to be regarded as exception- 
ally large.” 

The disintegrations by slow neutrons, which 
are presumably not of a radiative character, 
present no difficulties in interpretation, and are 
in fact analogous to those of lithium and boron 
by slow neutrons. The fairly large cross section 
observed for this process may here too be 
ascribed to the proximity of a virtual level of N™. 


OTHER MopbEs OF DISINTEGRATION 


Feather! and Kurie‘ reported another mode of 
disintegration in which one of the disintegration 
particles was singly charged. They suggested 
that these disintegrations took place according 
to reaction (2) or (3). With both the Be+,H? 
and the Li+,H? sources of neutrons we have 
observed forks of this type, an example of which 
is given in Fig. 1D. The data on the forks of 
this type, which were observed when the 


beryllium source was used, are given in Table III. 


1 Dr. Robert Serber has looked into this question theo- 
retically, and we are indebted to him for telling us of his 
results. One essential reason for the small theoretically 
predicted cross sections is the impossibility, with any 
reasonable model, of having approximate resonance with a 
virtual level of N® for the neutron, both before and after 
radiation. 
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TABLE II]. Ranges and energies of particles. 


RANGE 

SINGLY RANGE ENERGIES OF PARTICLES OF THESE RANGES 
CHARGED HEAVY En? Ex) 
PARTICLE PARTICLE Ew Ex: Ecu Ecu +Ecu +Ecu 








NEUTRON Source: Be+:H? 


2.42 cm 0.12 cm 1.04 1.4 0.2 0.3 1.6 1.3 
>1L8 0.32 >0.87 >1.1 1.3 1.4 >2.4 >2.3 
>3.4 0.20 > 1.28 > 1.6 0.6 0.6 >2.2 >1.9 
$2.1 0.22 >0.95 >1.2 0.7 OF >1.9 >17 

1.30 0.11 0.68 0.8 0.2 0.2 1.0 0.9 

2.26 0.09 1.00 1.3 0.2 0.2 1.5 1.2 
>1.4 0.18 >0.73 >0.8 0.5 0.5 >1.3 >1.2 

0.25 0.35 0.26 0.2 1.6 1.7 1.8 2.0 








From these data we cannot determine the 
maximum energy which appears in such forks, 
but we can say that at least 2 MEV appears 
when neutrons with energies up to 4.42 MEV 
are used. According to calculations from Bethe’s 
masses, reaction (2) is endothermic by 4.7 MEV. 
Thus we must turn to reaction (3) to explain 
the singly charged particles. 


7N4+ on'—,C“4+,H" (3) 


The C™ would probably be radioactive, going 
into N“ with the emission of an electron. How- 
ever, such a radioactive C™ has not been ob- 
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served. The upper limit of the mass of C™, 
assuming that 2.0 MEV of energy appears in 
reaction (3) when 4.42 MEV neutrons are used, 
is 2.8 MEV more than that of N". Thus the 
maximum energy of the beta-particles from C™ 
is less than 2.8 MEV. 

We have found one trident which we have 
attributed to the disintegration of nitrogen 
according to reaction (4). Fig. 2B shows a 
stereoscopic pair of photographs of the single 
disintegration of this type which we observed 
when we bombarded nitrogen with the high 
energy neutrons from Li+H*. The energy ap- 
pearing in the three forks is 4.8 MEV and the 
calculated energy of the neutron which produced 
the disintegration is 12.9 MEV. From a calcu- 
lation of the masses involved, this reaction is 
endothermic by 7.0 MEV, which is consistent 
with the observed data. : 

We are indebted to Professor C. C. Lauritsen 
and Professor J. R. Oppenheimer for valuable 
discussions and to the Seeley W. Mudd fund for 
financial support. 
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The Energy Distribution of Neutrons Slowed by Elastic Impacts 


E. U. CONDON AND G. Breit, Palmer Physical Laboratory, Princeton University 


(Received December 9, 1935) 


The problem of finding the distribution in energy of 
particles of mass m, initially. of the same energy, which 
have made m impacts with particles of mass M all initially 
at rest, is solved. It is supposed the impacts are elastic and 
the distribution in angle isotropic in a coordinate system 
in which the center of mass is at rest. If x is the ratio of 


N this note we work out the energy distribu- 

tion of neutrons which, starting with the 
same initial energy, have made m impacts with 
other nuclei all initially at rest. We suppose the 
impacts are elastic and the scattering isotropic 
in a coordinate system in which the center of 
mass is at rest. The result is of some interest in 
connection with current researches on “‘slow”’ 
neutrons. The work grew out of a desire to 
understand a statement due to Fermi! that “It 


1Amaldi, D’Agostino, Fermi, Pontecorvo, Rasetti and 
Segré, Proc. Roy. Soc. A149, 522 (1935). See p. 524. 


the energy after m impacts to the initial energy then the 
chance that x lie in dx at x is (log 1/x)""/(m—1)! for 
m = M. For unequal masses the expression is more compli- 
cated but easy to calculate. The results have some interest 
in connection with the slowing of neutrons by elastic 
impacts with other nuclei, especially with hydrogen nuclei. 


is easily shown that an impact of a neutron 
against a proton reduces, on the average, the 
neutron energy by a factor 1/e.” 

Let the nuclei of the medium be all at rest and 
of mass M while the incident neutron is of mass 
m and energy, Eo. Then by a simple application 
of the conservation laws it is found that the 
neutron energy after an impact is given by 
FE, = Eo(1 — ax) 


where a=4mM/(m+M)? and cos g=1—2x, 


y being the angle of scattering of the neutron in 
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a coordinate system in which the center of mass 
of m and M is at rest. For elastic spheres on 
classical mechanics or for short range forces of 
any type in quantum mechanics the chance that 
¢g lie between ¢ and ¢+dg¢ is proportional to 
d(cos g) so we see that x may take each value 
between 0 and 1 with equal probability. 

After collisions the energy will be EZ, where 


E, = Ey(1 —ax,)(1—axe)---(1—ax,). 


Let x=E,/E». We have to find the probability 
that x lie in dx at x being given that each x; 
has equal chance of having any value between 
zero and unity. Hence x may vary between 
(1—a)" and unity. Write 


(l—ax,)=e"“* and u=>ou; so x=e. 


Each u; varies between 0 and a=log (1—a)~". 
The chance of a given set of u; values is 


dx dxo- + -dx,=e~“du,dus: + «dup. 


The chance that the sum yw, lie in du at wu is 
therefore the factor e~“ multiplied by the chance 
that u lie in du at u assuming each u; to be 
equally likely to have any value between 0 and a. 
The evaluation of this turns out to be a very 
old problem, apparently first considered by 
Laplace,’ and the result is 


1 n 
n(u)du = ———— ju"! — (u—a)"—! 
f 
a"(n—1)! 1 


n ™ 
+( )(w—20) met (— pei ) 
2 n—1 


xX[u—(n—-1 a} 


with the understanding that the term involving 
(x—ka)"~' is to be assigned the value zero for 
x<ka. Thus the distribution function has a 
discontinuous (m—1)st derivative with discon- 
tinuities at the places x=ka with k=0, 1,2 --- n. 
A modern derivation of the result has been 
given by Rietz.* 


Laplace, Théorie Analytique des Probabilités (1820), 
pp. 257-263. 

*Rietz, Proc. Int. Math. Congress, Toronto 2, 795 
(1924). The problem is also discussed by H. P. Lawther, Jr., 
Annals of Math. Statistics 4, 241 (1933) whose Fig. 1 
graphs the distribution function for nm =1, 2, 4, 8, 16, 32. 
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Before we looked it up in the library we had 
worked out a solution of the problem which is 
enough different from those we have seen to 
make it worth communicating briefly. The first 
step is to consider the recursion relation con- 
necting the distribution function for »+1 with 


n 


that for m. Let u=}lu; and v=u+u,,4,. Let 
=1 


I,"(u) be the expression for f,(u) in the range 
(k—1)a<u<ka. The chance of u being in du at 
u and Un4; iN duy4; at Mas; is then J,"(u)dudu,,,. 
Changing the variables to v and u and summing 
over the values of u which lead to a fixed value 
of v we find the recursion relation 


(k-l)a ° 
Lea) = f Iria(u)dut f Ti"(u)du 
(k-—l)a 


v—a 


with the understanding that J)"(u) =0. 

Now let us consider the geometrical situation. 
We have to find the (n—1) dimensional measure 
of the intersection of the hypercube 0=u;Sa 
with the two hyperplanes }-u;=u and ou; 
=u-+du. Plainly it will be proportional to u"~'du 
for 0<u<a by a dimensional argument. As 1 
becomes a little greater than a the hyperplane 
passes m corners and so the expression u”~'du has 
to be corrected by subtracting off n(u—a)"~'du 
to allow for the part of the hyperplane that is 
outside the hypercube. When u becomes just 


n 
greater than 2a the hyperplane passes (*) 


, n aon 
corners so correction by { , }(u—2a)""'du is 


necessary. That the correction has to be added 
this time is readily seen by inspection in the 
three-dimensional case where the figure is easily 
visualized. This argument suggests the general 
form for J,"(u) and it is easily verified that it 
satisfies the recursion relation and is therefore 
correct when properly normalized so that its 
integral from 0 to ma will give the correct 
volume, a”. 

For u>na, I,,"(u) is a sum of terms in u 
(u—a)""!--- which is easily seen to vanish 
identically. In fact this polynomial may be 
written as [1—e~? }"u"-!/(m—1)! which must 
vanish because [1—e~?]" expanded in D con- 
tains only D* and higher powers. Here D is the 
differentiation operator. The ratio of the coeff- 
cients of (w—ka)"~' is fixed because the 


n—l 
’ 
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Fic. 1. Showing the probability as ordinate that a neu- 
tron have a fraction of its initial energy less than or equal 
to the abscissa after a number of impacts with protons that 
is marked on each curve. 


independent ratios must satisfy »+1 conditions. 
It is thus not even necessary to verify the 
recurrence relation for J,"(u). 

The case of most interest experimentally is 
that corresponding to the slowing of neutrons by 
protons. Here m= M and a=1 so a= ~= and the 
complications due to passing corners of the 
hypercube do not arise. All finite values of u 
are now less than a so only /,"(u) plays a role. 
The distribution function for u is therefore 
e~“u"—'du/(n—1)! or for x, 


F,,(x)dx = (log 1/x)""'dx/(n—1)! 


Thus the probability that x have a value between 
0 and ¢ is given by 


é 
P,(§) = f F,,(x)dx. 
0 


The integral involved here is an incomplete 
gamma function for which very complete tables 
exist.‘ We have prepared Fig. 1 in which P,(£) 
is plotted against a logarithmic scale of & for 
several values of m. There is not much point in 
extending the curves below §=10~-* since, for 
ordinary initial energies of neutrons, a value of 
§ much less than this brings the neutron energies 
down to thermal values where the assumption 
that the protons are initially at rest is no longer 
valid. The curves give a good idea of the way 
in which the energy is rapidly reduced by a 
moderate number of collisions. 

Although one can calculate the average of any 
power of x from the distribution function, it is 
simpler to do it from the expression for x as a 





* Karl Pearson, ed. Tables of the Incomplete Gamma 
Function (London, H. M. Stationery Office, 1922). 
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Fic. 2. Same as Fig. 1 except that it is for collisions with 
carbon nuclei instead of with protons. 


product of the m factors (1—ax,). Then we have 


, * fi—(l1—a)**'}" 
#-| f (1—ax,)rds =|) . |. 
0 a(s+1) 


In particular for the case of equal masses, a=1, 
and the ordinary arithmetic mean, s=1 we have 
%=2-". The statement of Fermi referred to 
above is verified for the logarithmic mean, 
exp (log x). 

The function f,(u) can be represented around 
its maximum approximately by e~*(“—"¢/ * as 
may be found by considering x*— # according to 
the above formulas or by empirical fitting of nu- 
merical computations. The latter indicate that 5.7 
gives a somewhat better approximation than 6. 

Calculation with the distribution function for 
unequal masses is quite simple. In Fig. 2 we 
give, as an illustrative example, the integrated 
energy distribution of neutrons which have made 
10 and 20 collisions with carbon nuclei (M=12). 
Comparison of Figs. 1 and 2 affords a striking 
indication of how little the carbon nuclei in 
paraffin contribute to the slowing down of the 
neutrons. 

The above discussion is not intended to give 
the distribution of neutrons slowed down by 
passing through a given thickness of paraffin— 
some of the emerging neutrons obviously perform 
more collisions than others. The distribution 
function considered here is nevertheless useful 
for approximate estimates when most neutrons 
can be considered to have performed the same 
number of collisions. 
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Nonorthogonality and Ferromagnetism 
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The calculations in Heisenberg’s theory of ferromagnet- 
ism have been questioned by Inglis and others on the 
ground that the error resulting from the nonorthogonality 
of the wave functions may possibly increase without limit 
when the number of atoms becomes arbitrarily large. In 
the present paper it is proved that this difficulty does not 
really arise. Semiquantitative formulas are given to correct 
for the error due to nonorthogonality, which is shown to 
be of the order 228? relative to unity, where z is the number 
of neighbors and 6 is the overlap integral (1). A supple- 


INTRODUCTION 


N an interesting paper, Inglis' has stressed the 

fact that all existing secular calculations by 
the Heitler-London method in a system with a 
very large number of electrons are subject to 
question because they assume that the wave 
functions of the different atoms are mutually 
orthogonal. The most notable calculations open 
to this objection are those in Heisenberg’s? 
theory of ferromagnetism. The apparent diffi- 
culty arises regardless of how small is the 
“overlap” or ‘“‘nonorthogonality” integral 


6=fSS ¥i(1)¥,(1)d0. (1) 


Indeed Inglis shows that if the number » of 
atoms is large, the secular determinant is so 
prolific in terms involving 6 as to suggest that 
in the final result the correction factor due to 
nonorthogonality may be of the order 1+6? 
rather than 1+6° as in the familiar two-particle 
case. (Cf. also our discussion following Eq. (8).) 
If this is really the case, the ordinary calculations 
on ferromagnetism, etc., become devoid of all 
meaning, since regardless of how small is 6, the 
error can be made arbitrarily large by taking 
sufficiently great. In a crystal, the number 1 is, 
of course, to be considered as enormous. 

To be sure, the success of Heisenberg’s theory 
of ferromagnetism is indirect evidence that the 
nonorthogonality catastrophe is only an apparent 
one. However, until a mathematical proof has 
been given, the theory can scarcely be said to be 


1D. R. Inglis, Phys. Rev. 46, 135 (1934). 
2 W. Heisenberg, Zeits. f. Physik 49, 619 (1928). 


mentary note is included on a new method of approxi- 
mating the partition function in Heisenberg’s theory. 
This approximation should be somewhat better than the 
assumption of a Gaussian distribution, but agrees even 
worse with experiment, provided one assumes orthogo- 
nality. Actually, the influence of nonorthogonality is 
sufficiently large to render uncertain any attempt to 
deduce exactly the critical conditions (minimum number 
of neighbors, etc.) necessary for ferromagnetism. 


on a rigorous basis. In the present paper we aim 
to supply such a proof. 

Let us consider a system of m electrons in n 
distinct states. Each state is to relate to a 
different atom. We shall assume that we are 
solving only the problem of exchange degeneracy. 
Every atom must then be in an s state, as 
otherwise there would be the spatial degeneracy 
to contend with. The secular equation associated 
with the exchange degeneracy is 


det. K=0, 
where K is the matrix 
K=> (V—WA®)P™, (2) 


Here W is the desired energy constant, P‘” is a 
Dirac permutation matrix,* and V“?, A“ are 
the integrals 


VO=S-- f(POv*)Hvdy,: - -dv,, (3) 
AW= SL - f(POW*)Wdo;- + -do,. (4) 


When preceding a wave function ¥*, the symbol 
P, of course, denotes a permutation operator 
rather than a permutation matrix, while H is 
the Hamiltonian operator. 

The sum in (2) is over the various possible 
permutations, and embraces only a finite number 
of terms since m is finite. Written out more 
explicitly, with the constituent cycles of each 
P“® specified, (2) is 


* For fuller explanation of the meaning of permutation 
matrices in connection with exchange degeneracy see P. A. 
M. Dirac, Proc. Roy. Soc. A123, 714 (1929), or Chap. XI of 
his Principles of Quantum Mechanics, first edition. 
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K= —W14+ Dail Vin —AG)W Pui 
+ ECV iin —Agiin WIP ijn) 
+ DLViincan —AgaanW Papuan 
+ YC Veiien —AGinn WIP jer 
+permutations of fifth and higher order. (5) 


A notation such as (ij)(klm) means i—j—i; 

k-l—m-k, and any given permutation is to be 

counted only once in the summation. Thus 

Pix; is not to be regarded as distinct from 
(Uk) (17) 

P.;) «uy, etc. The wave functions V to be used in 
(17) (AO) 

(3)-(4) are those of the simple product form 


W=y1(1)2(2)-+-vn(n), (6) 


since before solution of the secular problem, the 
proper linear combinations are not known. 
Because of the invariance of H under permu- 
tations, it matters not which particular WV is 
selected for use in (6). Thus the first two factors 
in (6) could equally well be y,(2)¥2(1). In 
writing (5) we have assumed that the additive 
constant in H (and hence in W) is determined in 
such a way that (3) vanishes for the identical 
permutation, and the coefficient of 1 in (5) is 
consequently simply—W. This point is discussed 
further in the appendix, and must be kept in 
mind when the exchange integrals (3) are com- 
puted, as they are not invariant of the origin 
for H. 

To show that it is not obvious that the 
nonorthogonality or A terms are of subordinate 
importance, let us try the approximation some- 
times made in the literature of stopping with 
second order permutations. Then (2) becomes 


K=>) il Vien —AgGinp W JP — W1. (7) 


To see how things go, we note that one root of 
(7) is easy to find, viz., that corresponding to the 
antisymmetric orbital solution. Here the P,;;) all 
have the characteristic value —1. Let us assume 
that the nonvanishing A;;;), V¢;;) all have common 
values 6°, J, respectively, and are }nz in number, 
as is the case when each atom is coupled to z 
neighbors. Then according to (7) 


W= —3nzJ/(1—4nz6?). (8) 


The nonorthogonality thus spoils the answer 
completely for large values of m. In fact, (8) 
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would give such an absurdity as W==< if 
6*=2/nz. We shall see later that in reality it is a 
poor approximation to stop with only second 
order permutations and that for large m, rather 
than use (8), it is a better approximation to 
disregard nonorthogonality completely, and em- 
ploy simply 


K=Dp€VioPun]-W1 (9) 


as the secular matrix. 


FACTORIZATION OF THE EXCHANGE INTEGRALS 
AND SECULAR DETERMINANT 


Let us consider an exchange integral such as 
Vinay Which is associated with a permutation 
which can be factored into simple permutations, 
with each cycle involving different atoms. In 
the computation of this integral, the most 
important portions of the Hamiltonian function 
are the part which involves the coordinates of 
the electrons on atoms i and j, and the part 
which involves the coordinates of the electrons 
on atoms & and /. This point is analyzed more 
fully in the appendix—it seems best not to 
consider it more fully at present since it is 
perhaps fairly obvious that if, for example, we 
are permuting electrons i and j, it is their part 
of the Hamiltonian function which gives the 
most important contribution to the exchange 
integral. This assumption provides a simplifi- 
cation needed for our later work by permitting 
factorization of the higher order exchange inte- 
grals into exchange integrals of lower order and 
nonorthogonality integrals. Thus it enables us to 


write ; 
Vinay = Vendant+ Vandda, (10) 


similarly 
J (ij) (KD (mn) = J p)AcG@nAcmna + V end Acmn) 


| 


+ VimmAcinAcen, r (I 1) 





V (i (imay = J (i) ACimny + Vimm Aci, etc. J 


The trick is now to note that if we accept (10), 
(11), then under certain conditions (2) or (5) is 
nearly the same, for large n, as 


[1+ DS AcpPuntDAuiwP aie 
+ DA unPunant: ++ LX ViaPua 
+> ViiinP ain +X VieiienPeijen —W1), (12) 
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for when we multiply out (12) and employ (10)- 
(11), we obtain all the terms of (5) and also 
certain extra terms. No permutations of the type 
(ij)(kl), etc., occur in the second factor in (12). 
Use of (12) will be warranted when these extra 
terms are not too abundant. They occur when- 
ever one or more subscripts is repeated in both 
factors in the multiplication. For instance, 


AcinP ua V ieyP ny = Ati V (ieyP (ajay (13) 


is a term included in (12) but not in (5). 

If (12) were rigorously true, the secular 
problem would be greatly simplified. Since the 
determinant of the product of two matrices is 
equal to the product of the determinants, we 
could equate the determinant of the second 
factor in (12) to zero. It is particularly to be 
noted that the only P’s above the second order in 
this factor are those involving permutations 
which cannot be factored into two or more 
distinct permutation cycles. Such permutations 
we shall call unfactorable permutations. On the 
other hand, all types of permutations are included 
in the first factor of (12). If it is allowable to 
discard the V’s associated with nonfactorable 
permutations, i.e., to set 


V aje = Vejen = °° > =0, (14) 


then the second factor of (12) reduces to (9), 
and the energy levels W are the same as those 
associated with a secular problem (9) obtained 
by completely disregarding nonorthogonality and 
higher order permutations. The whole demon- 
stration thus resolves itself into the consideration 
of the amount of error involved in the substi- 
tution of (12) for (5) and in the assumption (14). 
We shall discuss the latter first. It is to be 
emphasized that while (14) discards nonfac- 
torable permutations, it is very necessary in 
the proof to retain the V’s associated with 
factorable permutations, for it is their very 
existence which makes possible the cancelation 
of most of the effect of nonorthogonality. 
Although they do not appear in the second sum 
in (12), their effect has not been discarded, as 
they are yielded on multiplication of the two 
factors. 

Examination of the approximation (14). We 
may distinguish between two cases: first, one in 
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which every atom is coupled to every other 
atom, i.e., where V,;;) is of about the same 
order of magnitude regardless of the location of 
t and j. The second case is that where the 
important coupling is between adjacent atoms, 
so that we may take V;;;)=0 unless z and j are 
neighbors. Fortunately it is the second case 
which occurs in nature. Of course it is an 
idealization to consider only the coupling be- 
tween adjacent atoms, but the coupling energy 
decreases exponentially with the interatomic 
distance, and this fact makes the convergence 
essentially that characteristic of the second case. 

In the first case, the number of nonfactorable 
permutations is comparable with the factorable 
and (14) surely is an unwarranted approxi- 
mation. For instance, the number of permu- 
tations of type V (ix. is of the same order m‘ as 
the number of the form Vi; qn. The terms of 
type V¢ij.n are presumably individually of the 
order of magnitude V(;)4q. like Voijun, and 
the abundance of the former is so great that 
regardless of the smallness of A,:, they over- 
shadow the retained ordinary second order 
terms V;;, which are individually much larger, 
but only 3m? in number. 

In the second case, the nonfactorable permu- 
tations are relatively much less abundant than 
in the first. For instance, in linear chains, or in 
simple, face-centered, or body-centered cubic 
arrangements, there are no terms of the third 
order type V(;;«), as long as we are including only 
coupling between adjacent atoms. The number 
of terms of form V,;;,. is at most‘ of the order 
nz’, where z is the number of neighbors possessed 
by a given atom. On the other hand, the number 
of factorable fourth order permutations, such as 
Vinay, Which we retain, is of the order nz’. The 
number of ordinary second order permutations 
involved in (5) is mz. If Vi: «un is comparable in 
magnitude with V;;; Aun, the ratio of the effect 
of the neglected unfactorable permutations to 


‘Unfortunately it does not appear feasible to derive 
exact expressions for the number of discarded nonfactorable 
permutations appropriate to the various types of spatial 
arrangements. In the case of the linear chain, such permu- 
tations are entirely absent if we consider only coupling 
between adjacent atoms. For simple cubic gratings (z=6) 
there are 6n relevant permutations of type (ijkl). So we have 
probably overestimated the number of nonfactorable 
permutations of mth degree connecting adjacent atoms 
when we state that their abundance is of the order nz™”. 
However, it seems better to be on the safe side of things. 
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that of the ordinary V;;; terms is presumably of 
the order 2A,;. Similarly, the number of un- 
factorable terms of the mth order is roughly 
nz”, or less,* and the size of each term 
Vupdua®*. The abundance increases with 
order only as a power of z rather than m, and is 
offset by the occurrence of a higher power of 
Au. Hence we conclude that (14) is warranted 
if Ac;;) is small compared to 1/z. 


CALCULATION OF ERROR INVOLVED IN SUBSTI- 
TUTION OF (12) FoR (5) 


When the multiplication in (12) is performed, 
it is easily seen that the superfluous terms as 
compared with (5) are relatively small in number 
provided each atom is coupled only to its 
neighbors. For instance, there are about nz? 
terms of the unwanted structure (13), whereas 
there are nearly n*z*/4 terms of the desired type 
Aci VanPunPan, individually of the same order 
of magnitude V;;A;,. Hence it is reasonable that 
(12) is an allowable approximation. This argu- 
ment cannot, however, be regarded as entirely 
rigorous, for the extra terms are still exceedingly 
numerous in absolute numbers. A similar diffi- 
culty did not arise in connection with (14), as 
(14) is an approximation internal to the second 
factor of (12), whereas the first factor is much 
more prolific in higher order terms. Hence it will 
be safer and more illuminating if we examine 
quantitatively the difference between (5) and 
(12) by studying some specific examples in 
which the energy can be calculated explicitly 
both with and without the approximation (12). 

In all these examples we shall make the 
approximation (14), and also assume that the 
4’s associated with nonfactorable higher order 
permutations vanish, so that 


Acijey = Acijen = +++ =O. (15) 


Eq. (15) is the natural analog of (14), and re- 
quires that A“ vanish for any permutation 
which is not factorable into simple permutations, 
since the A“ for any differently factorable 
permutation have at least one factor of the form 
(15). For instance, Acij)jim)=AcijpAaim. The 
suppositions (14) and (15) can be regarded as 
slight specializations of the model as far as the 
comparison of (5) and (12) is concerned. Because 
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only a small percentage of the V's and A’s are of 
the form (14) or (15), the error involved in 
passing from (5) to (12) would very probably 
not be materially altered if different values than 
(14) and (15) were assumed. We shall suppose 
that all nonvanishing A;;;) and V;;;) have common 
values 6* and J, respectively. The notation 6* 
rather than 6 is used because a simple permu- 
tation introduces in (4) two one-electron inte- 
grals of the type (1), namely, one for each of the 
two electrons involved. The first and second 
expressions of (11) now become 3/Jé6* and 0, 
while more generally 


AMs§*, YOO weJ§5**-?, (16) 
if P‘ is expressible as the product of & simple 
permutations of neighboring atoms, with different 
atoms involved in each constituent cycle. Other- 
wise A“ and V“ vanish. 

Even with these approximations, the secular 
problem is one of hopeless complexity except as 
we confine our attention to particularly simple 
special solutions. One root of this character is 
that which corresponds to the antisymmetric 
orbital state or in other words to completely 
parallel alignment of spin. Here P“ has the 
characteristic value +1 or —1 according as it is 
an even or odd permutation. The solution of (5) 
thus becomes 

JT — vi +2026? —3v354+--- ] 
W= . (17) 
[1—v16?+9254— +--+] 


where », is the total number of permutations of 
degree 2k which can be factored into k simple 
permutations. The gist of Inglis’s observations, 
when applied to (17), is that since 1, v2, «++ are 
very large in a crystal, the numerator and 
denominator of (17) certainly have vastly 
different values than when one stops with the 
first term, so that it is not obvious that (17) has 
nearly the same value as when 6=0. The essence 
of our answering argument is that the terms 
involving 6 enter nearly as a common factor in 
numerator and denominator, so that the value of 
the ratio (17) is, for small 6 but arbitrarily large 
mn, not much different from that obtained by 
stopping with the first term. To see that this is 
really so, we now consider the following particu- 
lar models. 
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(a) Open linear chain. Here each atom has two neighbors, except the end atoms, which have only 


one. The », may be evaluated rigorously, and one finds that the denominator D of (17) is 
D=1—(n—1)6?+3(n—2)(n—3)64+---+(—1)*(n—k)---(m—2R4+1)5%*/R!I+---. 
This finite series is exactly summable’® to 


D=6" sinh (n+1)¢/sinh y, where cosh g=1/26. (18) 


It is not necessary, for any of the models, to give the explicit form, or mode of evaluating the 
numerator of (17), as, except for a factor J, the numerator is obtained merely by differentiating the 


denominator with respect to 6? and so presents no trouble. 
With the aid of (18), one finds that (17) becomes 


cosh* gf(m+1) cosh (m+1)¢ cosh¢g _— sinh ¢ 
W=-2)— Reeoeenes wien } (19) 
sinh ¢ sinh (n+1)¢ sinh ¢ cosh¢ 
For very large m, (19) is the same as 
W= —}385-*J[n(1 —46?)-!—n+ (1 —46?)-! — (1 —46?)-" ] = —(n—1)J(14+362+---). (20) 


The dots denote omitted terms of order 5‘ and higher. Now —(n—1)J is the value of the energy 
which would be obtained with neglect of nonorthogonality, i.e., by using (12) rather than (5). So the 
correction factor due to nonorthogonality is (1+36?+---), and not (1+é6?+---) as (8) would 
suggest. 

(b) Closed linear chain. A closely allied problem is that of the cyclic linear chain, which differs 
from (a) only in that there are no ‘“‘end”’ atoms. The series in the denominator of (17) is now® 


1—né?+43n(n—3)54+---+(—1)*n(n—k—-1)---(n—2k4+1)5*%*/k!+---=265" cosh ng, 


with ¢ as in (18). So (17) becomes 


W = —2J(cosh* ¢/sinh ¢)(n tanh ng—n tanh ¢), 
which for large n is 
W= —}n(J/6*)((1 —46?)-!—1 ] = —n V[14+-36?+--- ]. (21) 
The correction 1+36?+ --- due to including nonorthogonality is exactly the same as in (20), as one 


would expect since it is physically reasonably apparent that it doesn’t make much difference whether 
or not the chain is closed. The fact that (21) involves m and (20) contains m—1 is, of course, because 
the closed chain involves one more pair than the open chain. 

(c) Three-dimensional configurations. We now tackle the more physical problem where the atom has 
z rather than 2 neighbors. (The simple cubic arrangement has z= 6, the body-centered z= 8, the face- 
centered z= 12, etc.) We shall neglect the ‘‘edge effect’’ resulting from the fact that atoms at the edge 
of the crystal do not have their full quota of neighbors. The similarity of results with (a) and (b) 
makes it clear that this is legitimate. Even then, it would be exceedingly difficult to obtain, to say 
nothing of sum the series with, an accurate general formula for »,. However, a sufficiently good 
approximation is 


View = vel 3uz—2ko+k+O(Rk/n) |/(kR+1), (22) 


where the symbol O(k/n) means that terms involving 1/n and other negative powers of m have been 
neglected. 


®> For evaluation of the series see Bromwich, Theory of Infinite Series, pp. 177-178. I am indebted to my father, 
E. B. Van Vleck, for calling my attention to this reference. 
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NONORTHOGONALITY AND FERROMAGNETISM 


When we use (22) the denominator of (17) becomes 


D=1—4nz6?+ ---+(—1)**"(4nz)(3nz—2s+1)---(3ns—2ko+hk)6***?/(R+1) !+---. 


This series can be evaluated by the binomial theorem, and is merely [1—(2z—1)é? }"*/4=-®), 
Consequently (17) is 


W = —}nzJ/[1—(22—1)8? ]= —}nsJ[1+(22—1)6?+---]. (23) 


The correction factor for nonorthogonality again involves only z and not n. As a valuable check on 
the validity of the approximation (22), we note that if z=2, as in the linear chain, this correction 
becomes 1+36?+---, agreeing perfectly with (20) or (21). (The agreement, however, is lost when 
higher powers of 6? are considered, as further development of (21) gives 1+36?+ 1064+ 358°+ - -- 
whereas (23) yields 1+36*+ 964+ 276°+ - - -.) 

Note particularly that the overlap integral 6 enters in (20), (21), (23) only in its square. Conse- 
quently if, say, we have 6=0.1, a value not unreasonable in ferromagnetic applications, and if z= 
or 12, as in the ferromagnetic lattices, the correction due to nonorthogonality modifies the answer by 
about 20 or 30 percent—not enough to spoil the general trend of the Heisenberg theory. 

It is, of course, reasonable that the correction due to nonorthogonality should depend on the 
number z of neighbors rather than on the total number of atoms in the system. If this correction 
kept increasing with m, it would presumably be impossible to isolate different portions of the universe 
for even approximate calculations, and it is doubtful whether quantum mechanics would make sense, 
for regardless of how large were the mean distances between atoms, the results would keep depending 
on the total size assumed for the complete system. 

(d) Coupling between every atom. It is interesting to note that (17) can be evaluated accurately here. 
The denominator becomes 


D=1—3n(n—1)6?+---+(—3)*n(m—1)---(n—2k4+1)6**/k!4+--- =(2/5)"77,(1/8), 


where //,(x) is the Hermitian polynomial of order n. We are interested in the asymptotic value of 
H,(x) for large values of n, and for our purposes it is sufficient to use the approximation //,,(x) ~const. 
Xe!" cos [(2m)!x], which may be obtained by applying the W.K.B. expansion to the differential 
equation satisfied by H,(x). (If m is odd rather than even, the sine rather than cosine should be used.) 
The expression (17) becomes 

W=}J[ —né-*— 6-*+-(2n)'d-3 tan (2n/6*)!). (24) 


The result (24) fluctuates with m in an irregular, meaningless way, and is entirely different from the 
formula obtained with disregard of nonorthogonality. It is thus fortunate that the present case (d) 
is one of only academic interest. It is barely possible that the trouble would be less acute if calcu- 
lations could be made without the assumptions (14)—(15), which are entirely unwarranted in (d). 





Discussion of other roots than the antisymmetric 
orbital solution. We have apparently made our 
discussion very special in the preceding examples 
by considering only the solution which is 
antisymmetric in the orbits. Conceivably one 
might object that the nonorthogonality cor- 
rections are more important under other condi- 
tions. To show that this is not so, we may note 
that there is one more case which can be solved 
accurately,—viz., that obtained by using the 
symmetrical orbital solution. The latter obvi- 


ously doesn't satisfy the exclusion principle, as 
we can’t have all spins mutually antiparallel at 
once. However, for formal purposes we can 
imagine ourselves in a universe not governed by 
Pauli. One obtains the symmetrical from the 
antisymmetrical problem merely by changing all 
minus to plus signs in (17). Hence all our 
previous calculations (a), (b), (c), (d) can be 
immediately adapted if we simply change the 
sign of 6? and J. Thus the correction factor due to 
nonorthogonality becomes [1—(2z—1)é?+---] 
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rather than [1+(2z—1)6?+---] in case (c) and 
remains of the same order of magnitude as 
before. These two cases, the symmetric and 
antisymmetric, must be regarded as the two 
extremes, and in actual problems the error is 
probably intermediate between that in two 
limits. If one makes an approximation sim- 
ilar to that of Inglis' of keeping only diagonal 
elements in an m, system of representation (i.e., 
one in which each spin is space quantized 
separately), the calculations can be carried 
through as before even in the general case. The 
only difference is that J and 6? are replaced by 
qgJ and qgé*, where g is the fraction of spins that 
are parallel. The correction factor due to 
nonorthogonality is then [1+9(2z—1)6?+--- ]. 
Because of the rejection of nondiagonal elements, 
this formula is surely quantitatively wrong, but 
seems to indicate that the error due to neglecting 
nonorthogonality is intermediate between that 
in the two extremes. 

A somewhat better approximation for the 
general intermediate case is probably the follow- 
ing: assume that in each term (—1)*é** in the 
denominator and (—1)*6?*~* in the numerator of 
(17) the factor (—1)* is to be replaced by the 
mean value of the product Pye) ¢34)...¢...2%) of R 
simple permutations (each involving different 
cycles) for a state of spin S’ for the entire system 
(crystal). This mean can be calculated with the 
vector model, or otherwise.’ One finds that to 
a degree of approximation comparable with (22) 


a 
X[1+2n-'k(2+3x-!'+x-*) ], (25) 
where x= —}—2(S’/n)?, (26) 
and that 
W =4nsxJ[1+f(x) +--+] (27) 
with f(x) = (32+sx7!+-2x) 6. (28) 


ry 


(27) agrees with (23) for the saturation case 
S’=}n, x=—1. The factor $nzxJ is the usual 
Heisenberg mean value of W for a state of spin 
S’. The remaining factor 1+f/(x) cannot be 
regarded as an accurate correction for non- 





®* We omit details of the calculation; for the general 
method see pp. 340-342 of my Electric and Magnetic 
Susceptibilities. 
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orthogonality, as it is certainly not correct to 
replace the mean value of the quotient of two 
matrices by the quotient of the means. Never- 
theless, a crude quantitative correction for 
nonorthogonality is probably obtained if we first 
calculate an eigenvalue W by means of (9), 
and then take the corrected eigenvalue to be 
WL1+f(y)] with y= W/3nzJ. It is better to use 
y rather than x as an argument in (28), since then 
even in the correction term we allow somewhat 
for the spread in energy for states of a given S’, in 
virtue of which W is not the same for a particular 
solution as the mean value. 

The error involved in the approximation (14), 
which is probably less important than the 
correction studied in the preceding paragraph, 
can be avoided by using the secular determinant 
obtained from the second factor of (12) rather 
than from (9). After a solution of the resulting 
secular problem is obtained, the factor 1+/f(y) 
explained in the preceding paragraph is to be 
inserted. Even then, there still remains the error 
due to the approximation (10-11), which seems 
to be comparable with the other errors (cf. end of 
appendix) and for which a quantitative allowance 
does not appear easy. 


NOTE ON THE EVALUATION OF THE PARTITION 
FUNCTION IN HEISENBERG'S THEORY 
OF FERROMAGNETISM 


In closing, we shall digress briefly to consider 
another question than nonorthogonality. In 
Heisenberg’s theory of ferromagnetism, one is 
interested in the mean value of 


exp (—H/kT)=>0>,(—H/kT)*/qg! = (29) 


for a state of given spin. Here the Hamiltonian 
matrix is H= >> V;;;)P,;;) with the usual neglect of 
nonorthogonality and higher order permutations. 
Even with these approximations, the mean value 
of (29) cannot be accurately evaluated. It occurs 
to us that instead of using Heisenberg’s assumed 
Gaussian distribution,? a possible procedure is 
the following: To a degree of approximation 
comparable with (25) and (27), it can be shown 
by a calculation with the vector model,*® that 


He! =xJ[3ns+qa+O(k/n)1H*, (30) 


where a=1+3x~-'+2x-* and where x is defined 
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NONORTHOGONALITY 


as in (26). Here O(¢/n) denotes error of order g/n. 
If we employ (30), the statistical average of the 
series (29) may be evaluated by use of the 
binomial theorem. Then the mean value of (29) 
becomes 


exp (—H/kT) =exp (—nQ) (31) 
with 


Q=}2a~"' log [1+xak'T-'J(1+f(x)) ]. (32) 


In writing Eq. (32), we have multiplied J by 
the factor 1+/f(x) defined by (28) in order to 
make a crude correction for nonorthogonality 
suggested by the preceding section. If we set 
f(x) =0 and if we expand the logarithm in (32) in 
a Taylor’s series in J, carrying the development 
only to terms of order J, the result is the same 
as when the spread of energies for a given S” is 
disregarded. If terms in J® are included, the 
formula is the same as with Heisenberg’s Gaussian 
distribution. The magnetization associated with 
(31)-(32) is determined from the transcendental 
equation’ 


2S’/n=tanh [(8H/kT)+(2S’/n)dQ/dx]. (33) 


Here H is the applied magnetic field, 8 is the 
Bohr magneton he/4rmc, and 2.S’/n is the ratio 
of the intensity of magnetization 2.S’8 at temper- 
ature J to the true saturation magnetization 8 
achigved at 7=0. In order for the body to 
exhibit ferromagnetic behavior, i.e., in order for 
(33) to have a real positive root in 2.S’/n of the 


AND FERROMAGNETISM 239 


proper character when //=0, it is necessary for 
dQ/dx to be greater than or equal to unity at 
x=—}4. If we set f(x) =0, i.e., assume orthogo- 
nality, this condition can be fulfilled only if 
x=16, an impossibly high value. This result is 
disappointing, for (33) is presumably a better 
approximation than Heisenberg’s formula based 
on a Gaussian distribution, where ferromagnetism 
is obtained if z=8. Hence we give (33) only 
passing mention. However, it must be noted that 
the critical condition for ferromagnetism is 
sensitive to the correction for nonorthogonality. 
If, for instance, in (28) we take (s—}3)é?=0.13 
then (33) admits a ferromagnetic solution when 
z=10. Ferromagnetism is allowed for z=8 if 
6°=0.028: such values for 6? seem rather high, 
but are by no means inconceivable. Thus the 
nonorthogonality effects may explain away the 
dilemma, but this statement is not beyond 
question, since use of the factor 1+/(x) is only a 
partial and inaccurate correction for non- 
orthogonality, as explained in the preceding 
section. More likely, the trouble arises from 
omission of the O(k/n) term in (30).8 The one 
safe conclusion is that the influence of non- 
orthogonality may be sufficiently important to 
render very uncertain any attempt to deduce the 
minimum number of neighbors necessary for 
ferromagnetism. 

The writer has had valuable discussions with 
Professor J. C. Slater, Professor D. R. Inglis and 
Dr. H. M. James, for which he thanks them. 


ApPENDIX—Basis FoR Eg. (10) 


The Hamiltonian function for the m electrons is 
H =2,[(—h?/8x2m)V 2 —e2/ri'J4+2>,0k, j)+C, (34) 
with Q(k, 7) = —e?/ryi —e?/rj*§+2/rj,+0?/Ri*, 


where r;/ is the distance from nucleus j to the electron 
which is located on atom &. It is convenient for our pur- 
poses to label electrons according to the atom with which 
they are identified in the wave function W defined in (6) 
to which H is applied. The terms e*/r;, are the inter- 
electronic and e?/R?* the internuclear energy. We have 
assumed that each atom-ion exerts a Coulomb field on the 
electrons, but it is no trouble to adapt the analysis to the 
case of variable screening, wherein the atomic fields are 
central but not Coulomb. The bracketed part of (34) can 
immediately be dropped from consideration, as each factor 
in (6) may be assumed to satisfy an equation of the form 

™The passage from (31) to (33) is made most easily by 


expanding Q in a Taylor's series about the most probable 
spin, as explained on p. 330 of the preceding reference. 


[(—h?/8x2m)V 2 — e*/ri'— Wi i(i) =0, 


and so the effect of this part of (34) can immediately be 
included by shifting the origin for the energy by an 
additive constant = W;. This shift, and others, is covered by 
the additive constant C in (34) which, as stated after Eq. 
(6), is to be so chosen that the integral (3) vanishes for the 
identical permutation. 

As a typical case, consider the value of (3) for a permu- 
tation (12)(34). The portion of (34) which we desire to 
retain to get (10) is Q(1, 2)+0Q(3, 4). We must show that 
the effect of the remaining Q terms in (34) is unimportant. 
First we note that there are an enormous number of the 
form Q(k, 7), with k, 7>4. The contribution of the totality 
of these terms to the integral (3) is Aqs2)Acs4) U where 





5 The omission of the term O(k/m) in (22) is not nearly as 
serious as the corresponding omission in (30), for the series 
(17) encountered in connection with nonorthogonality are 
concerned with a correction effect, and need not be summed 
as accurately as (29). 
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Usa fies S2u>j>dk, j)dvs---don. 


Now the coefficient of Pi:2):3s) in the part of (5) involving 
W explicitly is — Aci2)(s4) W = — Ag2)Aca) W. Hence it is seen 
that the effect of the terms involving k>j>4 is covered 
simply by displacing the origin for the energy by an 
amount U, which we suppose done. (This is not quite the 
appropriate displacement in the case of higher order 
permutations. For instance, in the case of Py2)(s4)(58) one 
would wish to absorb in W only terms corresponding to 
k>j>6. However, except in the case of very high order 
permutations, which presumably involve very small inte- 
grals, the difference in the number of terms included is 
negligible in view of the fact that the number » of atoms 
is very large.) 
It remains to be shown that we can disregard terms of 
the form 
Q(3,1), Q(4,1), Q(3,2), Q¢4, 2) (35) 


and Q(k,j) where k>4 and j=1,2,3,4. (36) 


The effect of each individual term in (35)—(36) is small in 
absolute magnitude compared with that of Q(1, 2) or 
Q(3, 4) and presumably differs from the latter by a factor 
of the order 6, where 64 is defined as in (1). Namely, 
there are 2n—2 factors in (3) which are large in regions 
where Q(1, 2) or Q(3,4) are important, as is seen on 
substituting (6) in (3), specialized to P‘? =P) (94) and 
remembering that each y function is large only on its own 
atom. On the other hand, only 2n—3 factors are large 
where any term of the type (35) or (36) is important. 
Furthermore, if each atom is coupled only to its neighbors, 
the terms (36) outnumber Q(1, 2), Q(3, 4) only by a factor 
2s. Since z, the number of neighbors, is small (6, 8, or 12) 
for the various cubic arrangements, we are justified in 
neglecting (35), (36) if the interatomic distance is suff- 


FEBRUARY 1, 1936 


PHYSICAL 


ciently large. From the foregoing considerations, it appears 
that the error due to dropping (35), (36) is of the order 224 
when measured relative to unity. This fact at first sight 
seems rather disquieting, as the errors we estimate due to 
other causes are of the order 226° (cf., for instance, Eq. 
(23)). Thus the approximations (10) and (11) would 
appear to be the dominant causes of error, causing trouble 
unless 6 is exceedingly small. However, one fortunately 
finds that the order of the error can be reduced to 2:8? if 
in computing Va»), for instance, one takes the potential in 
the integrand of (3) to be not just K(1, 2)+2Wi+U+C 
but rather that 

K(i, 2)+2Wi+U+CS--+ Sf Zj>2LK(j, 1)+K(j, 2)]| ¥s!? 

»++!y/2dvg-+-dvn (37 


inclusive of contributions from the time exposure charge 
clouds of atoms other than 1,2. Such a procedure is 
essentially the analog of the Hartree method. The error 
resulting from (35) may seem to still be of order 6 rather 
than 6°, but as a matter of fact is entirely inconsequential. 
Indeed for a given permutation Py) ;34) it is altogether 
unlikely that 3 or 4 be a neighbor of 1 or 2, meaning that 
for the great bulk of permutations, terms of the type (35) 
are nonexistent so long as we are including only coupling 
between adjacent atoms. Hence the approximations (10) 
and (11) are no more serious sources of error than the other 
causes which we consider. 

It is interesting to note that no internuclear terms e?/ Ri 
contribute to (37), as their contribution to C just cancels 
their effect on the rest of (37). As already stated, the terms 
=W; and U also have offsets in C. 

In the preceding discussion we have considered for 
simplicity only fourth order permutations, but the argu- 
ment can be readily extended to others. 
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Quantum-Mechanical Description 


HENRY MARGENAU, Sloane Physics Laboratory, Yale University 
(Received November 16, 1935) 


It is pointed out that, by removing a single quantum-mechanical postulate commonly 
accepted, several conceptual difficulties in quantum-mechanical description can be eliminated. 


HE discussion of a recent paper by Einstein, 

Podolski and Rosen! has brought to light 
an interesting divergence of opinions as to the 
meaning of reality. To the writer it seems, 
however, that the issue at stake is more than 
the merit or demerit of any particular conception 
of reality, and that a very fundamental point of 


1 A. Einstein, B. Podolski and N. Rosen, Phys. Rev. 47, 
777 (1935); E. C. Kemble, ibid. 47, 973 (1935); A. E. Ruark, 
ibid. 48, 466 (1935); N. Bohr, ibid. 48, 696 (1935). 


quantum-mechanical axiology is involved. For 
Einstein, Podolski and Rosen have shown, by 
employing correctly the usual quantum-mechan- 
ical postulates, that the state of system 1 which, 
by hypothesis, is isolated from system 2, depends 
on the type of measurement performed on 
system 2. This, if true, is a most awkward 
physical situation, aside from any monstrous 
philosophical consequences it may have. We wish 
to show that, by the removal of a single postulate 
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commonly accepted, the real difficulty inherent in 
Einstein-Podolski-Rosen’s conclusion disappears. 


I 


The postulate in question may be stated thus: 
When a (single) measurement is performed on a 
physical system, then immediately after the 
measurement the state of the system is known 
with certainty. That is to say, if the system is in 
a state yg, the measurement causes the state 
function to transform itself into y, where y is an 
eigenstate of the operator belonging to the 
measured observable. This assumption is clearly 
involved in Einstein-Podolski-Rosen’s interpre- 
tation of the reduction of the wave packet, and 
if it be denied that in general a measurement 
produces an eigenstate, their conclusion fails, 
and the dilemma disappears. We shall refer to 
the postulate stated as I. 


II 


I is not self-evident, or a priori true; for there 
is hardly more justification for supposing that a 
single measurement determines completely the 
state after, than to suppose it to determine 
completely the state of the system before the act 
of measurement. The latter is usually denied; 
to admit the former supposition would introduce 
a peculiar asymmetry into quantum-mechanical 
description. But there are stronger arguments 
against I. 


Ill 


I is definitely contradictory to another, in a 
sense more fundamental, postulate, namely, 
Schrédinger’s time equation. To supply sub- 
stance to the following simple but formal con- 
sideration the reader may think of the system as 
a particle, of the state function as referring to a 
pure case with respect to momentum, of M as a 
coordinate measurement, and therefore of y as 
a 6-function. I implies that, under the act of 
measurement, g—y, or, to use customary ter- 
minology, the operator M converts ¢ into y: 


Me=y. (1) 
We observe that in this equation M cannot be a 
unique operator such as those with which we 
are usually confronted in quantum mechanics, 
for ¥y depends on the outcome of the measurement 


and cannot be predicted. That is to say, while we 

know, in our example, that y is a 6-function, 

the position of its peak is not fixed. ((1) must not, 
of course, be mistaken for a characteristic value 
equation, since ¢ and y are different functions.) 

Eq. (1) is sometimes interpreted differently. 

It is supposed that M selects from an assemblage 

of systems ¢ a subassemblage y and that it is 

not the total assemblage ¢ which undergoes a 

transition from a state to an eigenstate.? But 

there are difficulties in such a view, of which two 
may be mentioned. 

1. If ¢ truly referred to a larger collection of physical 
systems than y, the two functions would have to be 
constructed in different configuration spaces, and it 
would be difficult to give mathematical meaning to (1). 

2. It is certainly proper to make a measurement on a single 
system. In that case the above interpretation, which 
involves the splitting off of a subassemblage, is meaning- 
less. 


We thus conclude that (1), if it has significance 
at all, must imply the actual transformation of 
the state function for a single system into an 
eigenstate function for the same single system. 

On the other hand, the measurement is cer- 
tainly a physical operation, describable in the 
ordinary manner as an interaction between 
physical systems. As such it is subject to 
Schrédinger’s equation. Let the “‘natural’’ Hamil- 
tonian operator (the form valid for the isolated 
state) of the system be Ho, and let Hy represent 
the interaction with the measuring device (in our 
example coupling terms with the radiation field), 
so that 7=H,+Hy. Then 


He=(h/271)(d¢/dt). 


Thus if the measurement takes place in a small 
time At, we have 


Ag=(22t/h) AtH ¢. 


Since g+Ag is the function into which the 
measurement converts ¢g, and which we pre- 
viously called y, 


This relation is of the same form as (1). But it 
is unlikely to make y a 6-function even if the 
correct form for Hy were known. The essential 
contradiction, however, arises from the fact that, 


2 1 owe this criticism to Professor Kemble. 
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if we put 
M=1+(2ri/h) At(Ho+Hw) 


as (1) and (2) require, the left-hand side of this 
equation is an unspecifiable, nonunique operator, 
while the right-hand side is unique.* Since both 
(1) and (2) cannot be correct, we abandon (1). 


IV 


The assumption that a single measurement 
fixes a state either before or after the act of 
observation is untenable from the point of view 
of any reasonable probability theory. In contrast 
with classical physics, quantum mechanics de- 
fines its states in terms of functions, gy, in con- 
figuration space in such a way that there exists 
a unique (aside from arbitrary phases) corre- 
spondence between ¢ and a sequence of prob- 
ability amplitudes. More explicitly, ¢, and hence 
the state which it represents, is entirely equiva- 
lent to a set of numbers a; associated with a 
certain operator P. This set of numbers defines 
a probability distribution, whose elements are the 
observations on the observable p which belongs 
to the operator P, and whose properties are the 
eigenvalues of P.* A quantum-mechanical state 
is therefore synonymous with a probability dis- 
tribution. A probability distribution, however, 
cannot be fixed by a single measurement, but 
requires a very large number of observations. 
This is true even if the distribution consists of 
zero’s for all properties but one (quantum- 
mechanical pure state). 

To illustrate: In throwing a die, the probability 
distribution, which corresponds to the quantum- 
mechanical definition of a state, is the sequence 
of numbers {---3 for the properties 1---6. 
If a 5 turns up in a single throw, I can conclude 
that I have thrown a 5, but cannot infer the 
distribution }---%. This can be done only after a 
large number of throws. The fact that I have 
thrown a 5, while it has converted my ignorance 
with respect to the outcome of this observation 


* This situation is of course known. Von Neumann (Math- 
ematische Grundlagen der Quantenmechantk) deals with it 
most thoroughly. We find it difficult to accept his ex- 
planation for reasons to be discussed elsewhere. 

*For terminology see V. Mises, Wahrscheinlichkeits- 
rechnung (Deuticke, 1931). The present point of view has 
previously been presented in greater detail by the author in 
The Monist 42, 161 (1932). 
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into certain knowledge, has not changed the dis- 
tribution of the probability aggregate from which 
I started. 


V 


The removal of I renders impossible the 
preparation of states in the manner advocated by 
Dirac and others. But the situation is perhaps not 
altogether unsatisfactory on that account, for it 
is still possible to ascertain a state, namely, by 
performing a large number of observations. 
A typical instance of a complete set of simul- 
taneous observations which does define a state 
is the photographic record of a spectrum, where 
the intensities represent at once the distribution 
of probabilities among the various pure energy 
states. It must be admitted, however, that the 
relation between abstract states and experience 
is a source of difficulties which weaken the pro- 
posal here made and which must be subjected to 
closer scrutiny. 


VI 


While the previous arguments seem to show 
that I is in several respects undesirable and con- 
flicting with other axioms, one may also convince 
himself that it is unnecessary. For no significant 
quantum-mechanical calculation requires I. Space 
does not permit an elaboration of this assertion 
here.® ; 


VII 


Returning to the argument of Ejinstein- 
Podolski-Rosen, the solution of their problem in 
the light of the present considerations is clear. 
If a great number of measurements of A had been 
performed on the state of system (1) after the 
interaction, all the coefficients y, in Eq. (7) of 
Einstein-Podolski-Rosen’s paper would have 
been determined. Conversely, if many measure- 
ments had been made on B, all the ¢, in their Eq. 
(8) would be known. A difficulty can arise only if 
the expansion of the right of (7), with the experi- 
mentally determined y,’s, yielded a different 
function that does the expansion in (8) with 
experimentally determined g,’s. There is no 
indication that this will be the case. 


5 The author hopes to return to the entire problem ina 
paper to be published in J. Phil. Sci. 
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On the Significance of Slopes and Other Parameters Estimated by Least Squares 


W. Epwarps DeminG, Bureau of Chemistry and Soils, U. S. Department of Agriculture, Washington, D. C. 


(Received December 2, 1935) 


On the assumption that the errors of observation in both 
x and y coordinates are normally distributed at each point, 
and that the errors ¢q, «, €- in the estimated parameters a, b, 
cof the empirical formula F(x, y; a, b, c) =O are also normal, 
it is found that the ratio e,w,'/¢ is distributed on a Student 
(Pearson type VII) curve. W, is the weight of a, and ¢? is 
the weighted sum of the squares of the residuals. The 
general case is found to be the same problem as the simple 
one discussed by Deming and Birge and others. When a is 
known, the test of significance is direct with the normal 
probability integral. When o@ is not known, the classical 
procedure is to use an estimate derived from the fit of the 


1. INTRODUCTION 


HIEFLY through the work of J. Neyman 
and Egon Pearson! the distinctions between 
various statistical testing criteria are becoming 
resolved. One important point in the application 
of such tests concerns the significance of the 
parameters that have been estimated by ‘‘fitting”’ 
a formula to a set of observations. For instance, 
suppose that a set of observed points in the xy 
plane gives a good fit to the linear relation y=a 
+bx, a being the intercept, and b the slope, as 
estimated from the data. At the same time, let us 
suppose, there may be good reasons for thinking 
that the slope should be something else, say B. 
We wonder if the discrepancy between B and } 
is significant. As an aid in forming a judgment 
we may apply statistical tests, some of which will 
now be discussed. 


2. THE -DISTRIBUTION OF THE ERROR IN THE 
ESTIMATED PARAMETERS 


The observations, if devoid of error, would 
satisfy some relation such as 


F(x, y; a, B, y) =0, (1) 


a, 8, y being the true and unknown values of the 
parameters.” Tentatively at least, the form of the 
function F is supposed known, and it is called the 


1J. Neyman and Egon S. Pearson, Biometrika 20a, 263- 
294 (1928); Proc. Camb. Phil. Soc. 29, 492-510 (1933); 
Phil. Trans. Roy. Soc. A231, 289-337 (1933). 

? The problem will be restricted here to two dimensions, 
the xy plane, and to three parameters. The extension to 
more general conditions is immediate. 


curve, and, as Student and Fisher and others have insisted 
in other work, it is better to use the Student integral (the z 
test) rather than the normal integral (the u test). However, 
as is now pointed out, the classical procedure involves the 
assumption that the fit of the curve is an average fit, on 
which basis the numerical difference between the two 
integrals is not apt to be serious. It is important to remem- 
ber nevertheless that the test thus made is the z test and 
makes use of the diverging s contours on the e, s plane. Some 
distinctions between the u and z tests are drawn and 
discussed. 


“empirical” formula. Whether in the end this 
particular form F is to be rejected or retained 
may depend to some extent on how well it fits 
the data at hand. Quantitative tests for goodness 
of fit are provided by the chi-test and the analysis 
of variance.* Unfortunately, however, in many 
problems the circumstances will not permit either 
of these tests to be made. 

Now of course the observations are not devoid 
of error, so the values of the parameters as 
estimated from the data will differ more or less 
from the true values a, 8, y. If the estimates of a, 
8, y made from the data by least squares are 
denoted by a, b, c, then the errors €,, &, € in 
these estimates will be 


fa, & €=a—a, b—B, c—y. (2) 


It will be sufficient for the purpose to make the 
usual assumption that the errors €., @, € are 
distributed normally; then‘ 


Ddeq = [wa'/ao(2x)! Je~“a*s*/2#*de, (3) 


3A history of the chi-test — its application to curve 
fitting are given in two papers by the writer; J. Am. Stat. 
Assoc. 29, 372-382 (1934); Phil. ‘Mag. 19, 389-402 (1935). 
The analysis of variance will be discussed in forthcoming 
papers. R. A. Fisher first recognized in the x*-distribution 
an instrument for testing the form of an empirical formula, 
Proc. Roy. Stat. Soc. 85, 597-612 (1922); therein he also 
applied the analysis of variance to the same end. The term 
“analysis of variance”’ is due to Fisher. 

* Probability curves will be written in differential form, 
and D (connoting density) will be used indiscriminately for 
the ordinates of all curves. The quantity under the differ- 
ential sign specifies the abscissa of the probability curve 
that D is the ordinate of, and the whole expression (i.e., 
either side of the equation) is the probability in a differential 
range. 
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will be the chance that the error in a lies in the 
interval ¢,+3de,. Similar equations will hold for 
«, and ¢«.. The symbols wa, w#, w- will be used for 
the weights of the estimates a, }, c; the symbol 
o for the r.m.s. error (called by some writers the 
quadratic mean error) of observations of unit 
weight. 

The weights wa, w», w. of the parameters are 
determined in the process of solving the normal 
equations. It is of course necessary to have the 
weights of the observations at the outset; other- 
wise the least-squares procedure cannot be 
commenced.* What to use for ¢, the r.m.s. error 
of an observation of unit weight, is always a 
problem. Its value can, of course, be estimated 
from the fit of the curve to the data under con- 
sideration, that is, by external consistency, as 


Birge’ would say; however, as will be seen later, : 


such an estimate should be used with discretion. 
Other estimates of ¢ can sometimes be made, for 
example, from the fit of the same formula to the 
data of other but supposedly similar experiments; 
moreover, it can be estimated also from the con- 
sistency of observations of known weight made 
on single coordinates, that is, by internal con- 
sistency.” By an “‘a priori known value of o’’ we 
shall mean a statistically reliable estimate of o 
made from sources other than the fit of the curve 
to the data under consideration. Since weights are 
purely relative, o is not required for the least- 
squares solution; it is, however, needed for 
certain tests to be described later. 


3. THE u TEST FOR THE SIGNIFICANCE OF AN 
ESTIMATED PARAMETER 


If o is known a priori, we derive at once from 
Eq. (3) the probability with which an estimated 
parameter will be in error by an amount greater 
than some proposed amount, uv. Thus, suppose 
that there is good reason for thinking that the 
parameter a in Eq. (1) has the value A, but that 
upon fitting the formula to a set of data, we find 
the value a for our estimate of a. Is the discrep- 
ancy between A and a significant? 

5 See, for example, O. M. Leland, Practical Least Squares 
(McGraw-Hill, 1921), Article 146; A. De Forest Palmer, 
Theory ef Measurements (McGraw-Hill, 1912), Chapt. 9. 

*W. Edwards Deming, Phil. Mag. 11, 146-158 (1931); 
17, 804-829 (1934). 

7 Raymond T. Birge, Phys. Rev. 40 207-227; 228-261 
(1932). 
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By proposing the value A for a we in effect 
propose 
u=a—A (4) 


for the error in a. The chance P,, that an error in 
a will exceed the magnitude |x| is that portion of 
the area under the normal curve of Eq. (3) lying 
external to the abscissas +u. P,, is easily evalu- 
ated from tables by noting that it is also the 
area under the normal curve of unit standard 
deviation lying external to the abscissas +uw,!/c. 
If |u| is high, P,, is low, and vice versa. 

Considerations of the numerical value of P, 
will be called the ‘‘u test.”’ If P, turns out to be 
very low, e.g., 0.01, then in the long run only 
once in a hundred such sets of observations 
would the estimate a of a be in error by as much 
as or more than |u| =|a—A| by the operations of 
chance alone. In such an event we must admit 
that 


(i) this is an unusual set of data (one set in a hundred), 
and the error ¢«, in a just happens to be unusually large; 
or that 

(#%) the true value of the parameter a@ does not differ from 
a by as much as |a—A |—in other words, A as a proposed 
value for a must be regarded with suspicion; or that 

(tit) the assumed value of ¢ is too small. 


A high value of P,, simply means that there is 
no indication that there is anything wrong with 
the proposed value A for a; it does not mean that 
A is necessarily the correct value a. The other 
parameters 5 and ¢ are tested in an exactly 
similar manner. 


4. THE CuHI-TEST FOR GOODNESS OF FIT 


For the development and interpretation of the 
z test we need first to know how the sum of the 
squares of the residuals is distributed. In dealing 
with this problem we shall be led into the chi-test. 
By an x-residual V, is meant the x-component of 
the line segment joining an observed point with 
the calculated position of this point. A similar 
definition describes the y-residual V,. The situa- 
tion is depicted by figures in two papers by the 
writer. The method of least squares yields not 
only the estimates a, b, c of the parameters a, 8, 7, 
but also the residuals V, and V, at every point;* 
in fact by a well-known systematic scheme of 
solving the normal equations,* the quantity 


8 See for example O. M. Leland, reference 1, Article 147. 
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¢* = r(w.V,"+w,V,’) (5) 


can be computed directly without the inter- 
mediate calculation of the individual residuals. 
The summation for ¢? is to be taken over all 
points. But if at any point the y coordinate (e.g.) 
is alone subject to error, the term w, V, contrib- 
utes nothing to ¢*; and an analogous statement 
holds if the x coordinate is alone subject to error. 
If now x? be defined by the equation 


2= 7/07 =(1/07)=(w.V,*+w,V,’), (6) 


and if k denotes the number of observed points 
diminished by the number of adjustable param- 
eters (k is called the ‘‘number of degrees of 
freedom,”’ following R. A. Fisher), it is a fact 
that on the basis of the normal law of error, and 
on the assumption that the method of least 
squares is employed (the method of least squares 
is essentially a procedure for minimizing ¢? or x’, 
or of maximizing P,, soon to be introduced), the 
accidental fluctuations in x? follow the distribu- 
tion law* 


Ddx?=[1/1(3k)24*](x?)!4-Me“dy?,— (7) 


The mean of this curve falls at k; hence though 
x? may by the operations of chance take any 
value from 0 to ~, its mean value in a large 
number of experiments is k. Moreover, x? is 
confined statistically closer and closer to the 
neighborhood of & as & increases (i.e., as the 
number of observed points increases). 

The chance P, that a certain value of x? will 
be exceeded is that portion of the area under the 
chi-square curve of Eq. (7) lying beyond that 
abscissa x’. Given x’ and k, P, is easily evaluated 
from tables. If x? is high, P, is low, and vice versa. 
In case P, is unusually low, then it may be worth 
while to consider the factors that might con- 
tribute singly or severally to produce a high 
value of x. Thus if P, turned out in any problem 
to be 0.01 we should have to admit that 


(4) this is an unusual set of data (one in a hundred), and 
x’ just happened to be large; or that 

(tt) the assumed form F of the empirical curve is not the 
correct relation between x and y; or that 

(tt) constant or systematic errors are present in some of 
the points; or that 

(tv) the assumed value of o is too small. 


The chi-test is thus a statistical tool for testing 


goodness of fit and for detecting constant and 
systematic errors. 

A high value of P, (e.g., 0.2 or higher) simply 
means that there is no indication of trouble. It 
is not a definite indication that the form F of the 
empirical formula is the correct one, for just as 
by chance P, can be small when the formula does 
fit, so also by chance P, can be large when the 
formula does not fit. 


5. THe zs TEST FOR THE SIGNIFICANCE OF AN 
ESTIMATED PARAMETER 
It will be helpful now to introduce the symbols 
S, and 2,, defined by the equalities 


9 9 9 9 


= 0° / Wa = 0°X*/ Wa 


=(1/w.)=(w.V2+w,V,*), (8) 


re 


i 
Sa = €a/Sa=€aWa'/o. (9) 


The subscript a on e, w, s, and z is of course to be 
replaced by 6 when the parameter D is under con- 
sideration; similarly for c. However, to facilitate 
writing the equations that follow I shall fre- 
quently attach to these quantities no subscript at 
all, trusting it to be understood. 

From Eggs. (3) and (7) it is found upon multi- 
plication and replacing x? by ws?/o* that 


1 
993 


Daeds =| — ¢ we? “de 
) 


o(27r 


wik Ss k—1 
x| = (-) e7we? “as| (10) 
T'(3k) 28 De to 


must be the simultaneous distribution of the 
independent statistics « and s. Now from Eqs. 
(8) and (9), deds can be replaced by sdsdz. Upon 
making this replacement and allowing s (or ¢) 
to be unspecified by the device of integrating 
over s from 0 to ©, we find that 


Ddz=[1/B(4k, 3) }(1+327)'@+Pds = (11) 





must be the distribution of z. 

This is a ‘“Student”’ or Pearson type VII curve. 
It is symmetrical about the origin z=0. It be- 
comes normal and increasingly concentrated 
about the origin as k increases. More important, 
o has disappeared; the distribution of z is inde- 
pendent of o. 
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The chance P, that a certain value of |z| (i.e., 
of |€.\wa'!/¢) will be exceeded is that portion of 
the area under the Student curve of Eq. (11) lying 
external to the abscissas +s. Like P, and P,, P. 
is also easily read from tables. If |z| is high, P, is 
low, and vice versa. In the situation that P, is 
unusually low (e.g., 0.01), it must be that one or 
more of the following factors are contributory: 
(t) this is an unusual set of data; the error e, in a@ just 
happens to be unusually large, or (and this is important) 
@ happens to be unusually small—that is to say, the 
formula happens to give an unusually good fit; or that 
(it) the true value of the parameter a does not differ from 
a by as much as |a—A|—in other words, A as a proposed 
value of a must be regarded with suspicion. 


Thus the z distribution has led to another 
criterion for testing the significance of an esti- 
mated parameter. The z test, unlike the u and x 
tests, does not depend upon ¢. But another con- 
sideration turns up to take its place: among the 
possible explanations of a low value of P, there is 
the chance of an unusually good fit (small ¢), 
which did not need to be considered in examining 
small values of P,. 

The z test is called by Fisher the ¢ test; his 

= sk. 


6. DISCUSSION 


Apparently, from what has preceded, both the 
u and ¢ tests are criteria for the significance of 
an estimated parameter. Are the two tests equiv- 
alent? If P, is low, must P, also be low? The 
answer is no. They can be compared and studied 
by a simple graphical representation. 

The values of ¢, and s, that result from fitting 
a single set of data can be represented by a point 
in a two-dimensional diagram, ¢«, for one co- 
ordinate and s,{[=¢/w,' by Eq. (8)] for the 
other. If N sets of data are taken there will be N 
points on the diagram, and if N is very large the 
relative density of these points (number per unit 
area divided by JN) will be just the D of Eq. (10). 
An experimentally obtained plot of 2043 points 
is shown in a recent paper by Treloar and 
Wilder.® 

Let us pass for a moment to the simplest 
problem in curve fitting, namely, a set of n 
observations presumably all made on the same 


® Alan E. Treloar and Marian A. Wilder, Ann. Math. 
Stat. 5 324-341 (1934). 
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unknown magnitude yu. The » observations x, 
X2, **+, Xn, all of unit weight, are to be fitted with 
the simplest of all curves, x =a. It is not difficult 
to prove that the value of a which minimizes the 
sum of the squares of the residuals, =(x;—a)? or 
¢*, is a=Z, the arithmetic mean. Hence x= is 
the result of fitting the curve x=a to the n 
observations by least squares. The weight w, of a, 
i.e., of Z, turns out to be m, the number of obser- 
vations. In this simple problem, s, (or simply s) as 
defined in Eq. (8) is merely the r.m.s. residual, 
commonly called the standard deviation, of the 
n observations; moreover x? is then merely 
ns*/o?. In the notation of Deming and Birge,'® u 
denoted the true but unknown error —y in the 
estimated parameter a (i.e., in Z); hence their 
Eq. (11) for the simultaneous distribution of u 
and s must be just the density of points in the e, 
s plane for the simple case. We should therefore 
expect to find Eq. (10) in ¢« and s to be identical 
in form with Deming and Birge’s'® Eq. (11) in u 
and s, also of course with Neyman and Pearson's 
Eq. (v) in m and s in their 1928 paper,' and such 
is in fact the case. 

This means that a comparison of the u and z 
tests for the significance of parameters is no 
different in the general case from what it is in the 
simple case, which fortunately has been clearly 
presented by Neyman and Pearson.' The subject 
has also been discussed by others, among whom 
Treloar and Wilder® and Deming and Birge"’ 
have already been mentioned. It will therefore be 
sufficient here to give a brief summary of the 
theory. 

”W. Edwards Deming and Raymond T. Birge, Rev. 
Mod. Phys. 6, 119-161 (1934); Phys Rev. 46, 1027 (1934). 


The u, s surface and the tests resulting therefrom are dis- 
cussed on pages 130-139 of the former. 
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The least-squares solution of a single problem 
in curve fitting provides us with ¢ and w and 
hence with the s coordinate of a point in the e, s 
plane. The corresponding ¢ coordinate is of 
course not known; if it were there would have 
been no problem in the first place. We plot there- 
fore the “test point,’’ whose coordinates are 
¢/w' on thes axis and |u| = |a—A| on the « axis. 

Through the test point and its image there can 
be drawn the three straight contours 


AA and BB: |u| = const. 
EE: s (or ¢ or x?)=const. (12) 
CO and OD: 's| = |u| /s =const. 


If the test point is moved along the u contour 
BB, u and P, remain constant, but P, and P, 
vary from 0 to 1. If the test point is moved along 
the s contour OD, z and P, remain constant while 
P,, and P, vary from 0 to 1. If it is moved along 
the s contour EE, s (and hence also ¢ and x?) 
and P, remain constant while P, and P, both 
vary between 0 and 1. P, and P, are both unity 
when the test point is on the axis Os, but they 
will differ seriously at moderate values of |u 
except when the test point is near the particular 
s contour ¢?=ko*, which represents an average 
fit." Hence P, and P, will agree substantially 
only when the curve gives an average fit. 

P,, P,, and P, are, respectively, the fractional 
number of points that in the long run fall external 
to AA and BB, above EE, and external to CO 
and OD. To the question, ‘‘What is the probabil- 
ity that the true error |e,| is greater than |u| 
=|a—A|?” the answer is obviously P,, since it 
names the fraction of points that in the long run 
will have error greater than |ja— A}. 

A difficulty arises, however, from the fact that 
P, requires a knowledge of ¢, which is often 
lacking or unreliable. In such a situation it is 
necessary to select some contour other than 
|u| =const., and for this there are many possi- 
bilities, for besides the three straight contours 
shown in the figure a number of curved contours 

4 Approximately, but not exactly, half the points in the 
¢, s plane will in the long run lie above the s contour 
¢*=ko*. The mean value of ¢? is ko? because the mean 
value of x? is k; but the median value of ¢? is not exactly 
ko*, since the x? distribution, Eq. (7), is not symmetrical. 


However, the asymmetry disappears and the mean and 
median approach each other as & increases. 


could be drawn. But Neyman and Pearson,' by 
reducing the problem to one in the calculus of 
variations, have shown that of all these possi- 
bilities, the z contours and the sz test, in the 
absence of o, give the greatest consideration to 
alternative values that might be proposed for a. 

The classical procedure of examining estimated 
parameters—to describe it in the present ter- 
minology—is to replace o? by ¢?/k in the normal 
distribution of ¢« given in Eq. (3), which pro- 
cedure is equivalent to assuming that the fit of 
the curve is average. The quantity we®/o* in the 
normal curve Eq. (3) is thereby replaced by 
wek/¢ or 2*k (= Fisher’s #?), which from Eq. (11) 
is distributed, not normally, but in a Student 
curve; hence the classical procedure of using the 
normal probability integral should, in strictness, 
be replaced by the use of Student’s integral. It is 
to be remarked, however, that as far as numerical 
values are concerned, the distinction is not of 
great moment, since, as was remarked earlier, P,, 
and P, are not greatly different numerically along 
the s contour ¢?= ko’, which is precisely the rela- 
tion that is forced by the classical procedure. 

The point to be emphasized is not so much 
that Student’s integral should replace the normal 
integral in the classical procedure, but that the 
results therein obtained are values of P., not Py. 
It is not a question of accuracy in the integral, as 
seems to be the main contention in a recent paper 
by T. E. Sterne,” who worked out the distribu- 
tion of z (Sterne’s #/(m—m) is equivalent to the 
present 2”) for a certain type of empirical formula. 
The integral that he calls P (page 569 of reference 
12) is the present P,, but he has discussed it as 
though he had obtained a refined value of P,. 

In closing I should like to say that although 
Eq. (11), being the distribution of the generalized 
z defined by Eq. (9), probably includes all previ- 
ously related results as special cases, I do not 
regard it as anything fundamentally new over 
Fisher’s paper* of 1922. The main purpose at 
present is interpretation. 

It is a pleasure to express my indebtedness to 
Professor Raymond T. Birge of the University of 
California for his stimulating interest and for 
keen criticisms of an earlier draft of this paper. 


2 T. E. Sterne, Proc. Nat. Acad Sci. 20, 565-571 (1934). 
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The Tension Coefficients of Resistance of the Hexagonal Crystals Zinc and Cadmium 


MILDRED ALLEN, Physics Department, Mount Holyoke College 
(Received November 27, 1935) 


The tension coefficients of resistance of the hexagonal 
crystals, zinc and cadmium, have been measured. These 
coefficients have been found to be independent of the 
secondary orientation. This is in agreement with a theory 
set forth by P. W. Bridgman, as corrected by J. W. 
Cookson; furthermore the experimental points lie within 
experimental error of the curves given by the theory 
representing them as functions of the primary orientation. 


N two previous papers':? a study has been 

made of the tension coefficients of electrical 
resistance of the trigonal crystals, bismuth and 
antimony, as functions of the primary and 
secondary orientations of the crystals with re- 
spect to the cylindrical axis of the casting. The 
findings of these experimental papers have been 
in agreement with a theory of the change of 
resistance produced by stress developed by P. W. 
Bridgman® on the basis of crystal symmetry. 
The experiments presented here were under- 
taken in order to extend this work to crystals of 
a different type of symmetry and so to verify 
the theory further ; hence the hexagonal crystals, 
zinc and cadmium, have been studied. In the 
meantime, a correction of Bridgman’s theory 
has been made by John W. Cookson‘ and it must 
be examined to see how far it invalidates or 
changes the interpretation of the already pub- 
lished results for bismuth and antimony. Inas- 
much as hexagonal crystals may be considered 
as a degenerate case of trigonal symmetry, the 
modification of the theory involved will first be 
shown as it applies to the trigonal case. The 
elastic constants enter into the generalized equa- 
tions of Hooke’s law in the same way in both the 
trigonal and hexagonal cases, but in the hex- 
agonal case one of the six elastic constants of 
the trigonal case is zero. That means that to 
determine the change of resistance with stress 
fewer constants will be necessary for hexagonal 
than for trigonal crystals and that to pass from 
the trigonal to the hexagonal case it is only 
necessary to equate to zero one or more of the 
coefficients involved. 
"1 Mildred Allen, Phys. Rev. 42, 848 (1932). 

* Mildred Allen, Phys. Rev. 43, 569 (1933). 


+P. W. Bridgman, Phys. Rev. 42, 858 (1932). 
* John W. Cookson, Phys. Rev. 47, 194 (1935). 


Cookson’s correction gives the same formal dependence of 
the coefficients on the orientation as Bridgman’s original 
theory, but demands additional constants fully to define 
the effect of deforming forces on the electrical resistance. 
Further experimental work will be necessary completely to 
determine this entire set of constants, and theory is 
sketched to show that torsion experiments would furnish 


the necessary data. 


CORRECTION OF THE THEORY FOR 
TRIGONAL CRYSTALS 


In his theory P. W. Bridgman assumed that 
the coefficients p,, relating the change in re- 
sistance to the stresses should behave exactly 
like the elastic constants relating the strains to 
the stresses, that is, that they should be repre- 
sentable by a symmetrical matrix in which 
Prs= Per. In January, 1935 Cookson published a 
note in the Physical Review showing on the basis 
of the known symmetries of the tensors, tension 
and electrical resistance, that Bridgman’s theory 
was incorrect in that the matrix of the piezo- 
resistive coefficients was not necessarily sym- 
metrical and that as a result more constants are 
necessary to define the change of resistance under 
tension than Bridgman indicated. In the case of 
trigonal crystals, two more constants are neces- 
sary making a total of eight instead of six. 

In comparing the Bridgman and Cookson 
coefficients it is well to note that they are defined 
differently by the two men. The coefficients pj), 
Pi2, pis and p33 are the same in the two cases; 
the Bridgman coefficients py and p44 are twice 
as great as those defined by Cookson. In accord 
with the earlier results reported, those used by 
Bridgman will continue to be used in this work. 

According to the original Bridgman theory the 
tension coefficient 8 depends on the primary 
orientation @ and the secondary orientation ¢ in 
the case of a trigonal crystal in the following 
way, the trigonal symmetry showing up in the 
term in ¢: 


B= (1/pe) {p11 sin* 6+ ps3 cost 0 
+ (213+ p44) sin® 6 cos? @ 


—2p.4 sin’ 8 cos 0cos 3g}. (1) 


248 





am Ah 


19 


im 


e of 
inal 
fine 
nce, 
y to 
r is 
lish 


at 
re- 
tly 





TENSION COEFFICIENTS OF RESISTANCE 249 


The corrected theory gives the modified relation 


p= (1 pe) Pil sin‘ 0+ p33 cos* 6 
+(pis+psi+p44) sin? 6 cos* 6 


—(pist+pai) sin’ @cos @cos3¢}. (2) 


It is to be noted that 6 is the same formal 
function of the orientation angles as before, and 
that therefore the conclusions drawn regarding 
the fact that the tension coefficients show the 
characteristic symmetries of the crystal still 
hold. In my second paper,’ however, an attempt 
was made to evaluate the piezo-resistive coeff- 
cients numerically and only part of the numerical 
results can be accepted under this new interpreta- 
tion. Eqs. (2) on page 571 of that paper are 
empirical data drawn directly from the experi- 
mental points and the functional relation (1) of 
this paper; the functional relation (2) of the 
present paper can be equally well satisfied if 
there be substituted in these equations for 
(2pist+ pas) the new relation (pi3+p3:+ 44) and 
for pis the expression (pi4+4:)/2. Since these 
data for bismuth do not appear explicitly in 
either of the earlier papers, they are given here: 
Uncorrected 
pu =(—7.2+0.16) xX 10-° 
pss = (—6.2+0.41) x 10~° 


(pis t+ pat pss) = (—7.140.99) x 10-9 
(pist+pn)/2 =( 13.540.27) X10~° 


Corrected for strain (3) 
pu = (—7.740.15) X10 9 
pss = (— 6.8+0.40) x 10 9 


(pis t+ psi + pas) =(-—9.3 +0.95) x 10-° 
(pis t+pu)/2=( 13.240.26) x 10°. 


It is obvious that as before the experiments with 
linear tension can be expected to furnish only 
four relations between the coefficients. 

In the previous work two additional relations 
between the constants were obtained from the 
values of the pressure coefficient of resistance for 
the two extreme orientations of the crystals, i.e., 
with the corrected notation, 


1 AR 
(- —) = —(p11+pi2t+pis) ‘Pe, (4) 
Pb RJ 690 


1 AR 
(- ~) = —(2p31+ p33) / pe. (5) 
bP Reno 


Even with these, there is now an insufficient 
number of relations to determine all the coeffi- 
cients involved. Thus the values of the individual 
p’s given in the paper on antimony are invali- 
dated (i.e., Eqs. (6) and (27)). 


FURTHER THEORY 


This insufficiency in the number of relations 
then requires the finding of additional data to 
give further relations between the constants. 
The most obvious procedure is to twist the 
crystal about its longitudinal axis. The resulting 
change in resistance must show the characteristic 
symmetry, if Neumann’s general law is to hold. 
This has not been done experimentally, but 
perhaps it is worth while to sketch the underlying 
mathematics and so to indicate the experimental 
problem that must be solved before these piezo- 
resistive coefficients can be completely known. 

A moment of force N’ is to be applied to twist 
the crystal about its axis of figure, i.e., about the 
Z’ axis. Then N’ will determine the values of 
the shearing forces according to the formulae 
given by Voigt,® 


X'  =2N'y'/Qa’, (6) 
y’ =— 2N’x’ Qa?, (7) 


where Q is the area of the circular cross section 
of the cylindrical casting and a the radius of 
this cross section, and x’ and y’ are the distances 
along the coordinate axes from the shearing 
force to the axis of the cylinder. It is to be noted 
that these shearing forces increase proportion- 
ately to the distance from the axis of the cylinder, 
so that the change in resistance will not be 
constant throughout the crystal, but those fila- 
ments parallel to the axis and’ near to the axis 
will change resistance relatively little compared 
with those far from the axis. The total change in 
resistance will then be the resultant of the 
changes of resistance for these individual fila- 
ments of small cross section. If f:(6, ¢) represents 
the shearing coefficient of resistance (1/X’,-) 
X (AR/R),; for one of the individual filaments of 
length / and of cross section dx’dy’ twisted by a 
force X’, and f2(@, g) the shearing coefficient for 
the same filament acted upon by the force Y’,,, 


5 W. Voigt, Letrbuch der Kristallphysik (1928 edition), 
p. 636. 
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the total relative change in resistance for that 
filament will be 


(AR/R);=fi(0, o)X' 2 +f2(0, oe) Y's. (8) 


The problem then is to add up all the changes in 
resistance for these various filaments which are, 
of course, in parallel with each other. In the un- 
strained state the usual law for the addition of 
resistances in parallel gives for the value of R, 
the resistance of the entire crystal, the relation 


1 1 
o_o — (9) 
R f R; 
and in the strained state there is the similar 
relation 
1 1 
———=}) ——_ (10) 
R+AR *f R;+AR, 


The latter may be expanded by the binomial 
theorem, since AR/R is a small quantity both 
for the individual filaments and for the entire 
crystal, to give the relation 

AR; ‘) 

R, 


1 AR 1 
R R f R; 


which becomes, in the light of (9) and of the 
fundamental definition of specific resistance as 
applied to the filament of length / and of cross 
section dx'dy’, 

AR fR 

= ffiX' et fe” 

R_ / pol 
The substitution in this of the values of the 
shearing forces given by Eqs. (6) and (7) gives 
for the total relative change in resistance the 
relation 


(11) 


)dx'dy’. (12) 


i 


AR 2N’ ih—o%8 
—=f1(4, o -f y'dy'dx'’ 
R Pa? 0° 


0 


2N’ a(a2 —y 2)4 
— f2(0, ¢) —f J x'dx'dy' (13) 


and, when the indicated integration has been 
carried out, the torsion coefficient of resistance 


becomes 


1AR 2 2 
=—/f,(0, ¢) —-—f2(0, ¢), (14) 
3 3r 
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where F=N'/aQ (15) 
is the force per unit area applied tangentially to 
the periphery of the crystal to produce the 
twisting moment N’. 

The problem then resolves itself into finding 
the functions f:(@, ¢) and f2(@, ¢) giving the 
relative changes in resistance of the individual 
rectangular filaments when subjected respec- 
tively to the shearing forces X’, and Y’,,. 
The primed axes are the axes of figure and the 
unprimed axes the crystallographic axes, with Z 
the trigonal axis and X the digonal axis, as 
before. The known laws of the transformation of 
tensors give the components of the shearing 


force X’,, in the unprimed set of axes: 

X,=2X’', sin @sin ¢ cos ¢, 

Y,= —2X’. sin @ sin ¢ cos ¢, 

Z2,;=0, 

: ‘ (16) 

Y.= —X’, cos @sin ¢, 

Z,=X'» cos 8 COS ¢, 

X,=X’', sin 6(cos? g—sin? ¢). 
These values must be substituted in the six 
fundamental equations for the change in re- 
sistance due to the application of stresses to a 


trigonal crystal 

bri =X out Yyor2t+Zepist VipisatZ0+X,0, 
bree=X 2012+ Yyout+Z-p13— Y.pist+Z0O+X,0, 
X p31 + Yyp31 +Z-p23+ YO+Z0+4X,0, 


6733= 


P41 Pal P44 
6fo3 = X,— — Y,—4+Z204+ Y:—+2Z04+X,0, (17 
? ? ? 


6r31= D 4 0+ Y,0+Z.: 0+ Y 04+Z. + +X yp 


6ri2= X,0+ Y,0+Z.0+ YV.0+Z.put+Xy( P11 — Piz). 


The total change in resistance is then to be found 
from the relation 


AR = dr; sin? 6 sin? g+ ree sin? @ cos? ¢ 
+ dr33 cos? 6+ 26re3 sin 8 cos 6 cos ¢ 
+26r3; sin @ cos @ sin ¢ 


+2érj2 sin? @sin gcos g, (18) 
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which gives for the final shearing stress coefficient 
the value 


1 (—) 
(8, ¢) =——{ — 
fi : Ye R Sf 


= | pis +2p4:} sin? 6 cos 6 sin 3¢. (19) 
pe 
This is seen to have the necessary trigonal 
symmetry in ¢. If the same procedure is followed 
through for an applied shearing force Y’,-, the 
shearing stress coefficient comes out to be 


1 AR 1 
fo(, ¢) -(- ) =—[(2p11— 2pis— pas) 
R Sf Pe 


Xsin*® 6 cos 0+(2p3:+ p44) sin 8cos*? @ (20) 
+ | pia(sin® 6—cos* 6) — 2p4; cos? 6} sin® 6 cos 3¢ }. 


Thus the final shearing coefficient of resistance 
due to the application of a twisting moment N’ 
becomes (according to Eq. (10)) 
1 AR 2 

— — =——[ (pu+2p4;) sin? 6 cos @ sin 3¢ 
N'/aQ R= 3xpe 

+ (2pi13+p4s— 211) sin*® 6 cos 0 

“= (2p31+ pas) sin 6 cos*® 6+ }2pa1 cos? @ 

— pi4(sin? 6—cos? @)} sin? @cos3¢]. (21) 


This gives five additional relations between the 
piezo-resistive constants which are more than 
are necessary. Thus for a complete solution of 
the problem of the relation of the change of 
resistance to stress, experiments involving the 
twisting of the crystals may be performed. 

The change of resistance due to the change of 
dimensions will be more complicated to compute 
than in the case of simple longitudinal tension, 
chiefly because the coefficients involved trans- 
form in a much more complicated way. 


EXPERIMENTAL PROCEDURE 


The study of the hexagonal crystals, zinc and 
cadmium, was undertaken for two reasons. The 
first reason, as has been indicated, is that hex- 
agonal crystals may be considered as a de- 
generate case of trigonal crystals in which the 
theory already developed may be applied if one 
puts pi4=p4,;=0. On examination of the above 
theory, it will be seen from Eq. (2) that the 
tension coefficient for hexagonal crystals must be 


bh 
wn 
— 


independent of the secondary orientation ¢. 
The second reason is that all real metals crystal- 
lize either in the hexagonal or the cubic system, 
and the results for these two metals will con- 
sequently be representative of those for a con- 
siderable number of other metals. Bismuth and 
antimony, on the other hand, are only quasi- 
metallic, having a higher resistance than real 
metals and with arsenic the only other element 
having the same type of symmetry. 

The experimental procedure with zinc and 
cadmium was the same as that used with 
antimony, the change in resistance being meas- 
ured by the deflection method. These hexagonal 
crystals required greater care, partly because the 
effect to be measured is so small and partly 
because the crystals bend so very easily. The 
order of magnitude of the tension coefficient 
proved to be the same as for antimony, but since 
the specific resistance of zinc and cadmium is 
only about one-sixth that of antimony, the 
change in the galvanometer reading produced by 
the application of the tension was correspond- 
ingly smaller, and in fact varied between the 
limits 0.2 and 1.2 cm. The uncertainty in the 
measurement of the cross section of the crystal 
cylinder surely introduced as great an error as 
arose from the smallness of these galvanometer 
readings, since the diameters were 1/16 inch and 
the cross section was never accurately circular. 

The zinc used was spectroscopically pure and 
was furnished by the New Jersey Zinc Company. 
The cadmium crystals were made from the pure 
Kahlbaum metal. Both were cast by the writer 
in cylinders 1/16’ in diameter and 1}” or 2” 
long. The melting points of both these metals 
were well below that of Pyrex glass, so that they 
could be crystallized in glass molds in an atmos- 
phere of CO, to prevent oxidation as they were 
lowered through the furnace. The metals would 
have stuck to the molds if these had not been 
washed out previously with a solution of lacquer 
in chloroform. The elastic limit of both is low: 
500 grams tension, approximately 250 g/mm?, 
for many orientations could be safely applied to 
the zinc crystals and 100 grams to the cadmium. 
The allowable force was determined by plotting 
the change in resistance against the force applied 
and finding how far the linear relation between 
the two held. It would be of interest to study 
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what happens beyond the ela:tic limit, but this 
has not been done as yet. 

Zinc and cadmium are both normal as regards 
the sign of the pressure coefficient, it being 
negative in both cases. In this they differed from 
bismuth and antimony for which the pressure 
coefficient of resistance is positive. These ab- 
normal elements showed abnormality also in the 
sign of the tension coefficient of resistance which 
was negative for some orientations and positive 
for others. Cadmium and zinc, however, are 
entirely normal in that the tension coefficient is 
found to be positive for all orientations. 


Zinc RESULTS 

Fifteen zinc crystals were measured; two 
others were measured but discarded, the values 
of the coefficients found being nearly twice that 
to be expected in comparison with the other 
values: such high values result easily from 
slightly cracked crystals. 

According to the theory sketched above the 
tension coefficient 6 of electrical resistance will 
be related to the primary orientation 6 (@ being 
the angle between the hexagonal axis of the 
crystal and the cylindrical axis of the casting) 
by the equation 


B=(1/pe)} p11 sin* 6+ 933 cos* 6 
+(pis+psit+pss) sin? 8 cos? @}. (22) 


(The specific resistance for the orientation @ is 
indicated by the symbol pe, and pi, p32, p13, pst 
and p4, are five of the six constants relating the 
change in resistance to the applied stress.) 
Thus the data from the tension experiments 
furnish three relations between the six constants 
involved. Applying the method of least squares, 
the values of these three constants, using the 
observed experimental data without correcting 
for the change in resistance resulting from the 
change in shape of the crystal, are the following : 


p11 =(0.86+0.26) X10", 
ps3= (14.09+0.45) X10-", (23) 
(p13-+ 31+ pas) = (15.27+1.2) K10-", 


where the tension 7 is expressed in units of 
kg/cm*?. The values of these same constants 
corrected for the change in shape of the crystal are 
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Pi= (0.00+0.29) x 10 a 
p33 = (11.634+0.50) K10-", (24) 
(pis t+psit+pas) = (15.19+1.35) K10-". 


The resulting empirical curves for 8 are shown 
in Fig. 1 together with the individual observed 
points. The difference between the points and 
the curve is in some instances as great as 
2x<10-*, but in view of the small deflections 
which were measured this is well within experi- 
mental error, so that the agreement seems satis- 
factory. The average difference between the ob- 
served points and the curve is about 1.1 X 10~°. 
It is of interest to compare these results with 
those found by other observers for polycrystalline 
zinc. Tomlinson® in 1883 found the value for the 
tension coefficient of zinc to be 2.75 X10~* and 
Rolnick’ quite recently, using a.c. technique, 
found the value 3.26X10~-*. Both these values 
lie within the range of those observed with single 
crystals of differing orientations. If the formula® 
for the tension coefficient of a haphazard arrange- 
ment of crystals of all possible orientations be 
applied, i.e., 
= (25) 


B 1 [== =] 


p 3 Pe=90 Pe=0 


the average tension coefficient comes out to be 
7.8X10-*. p is here the specific resistance of the 
haphazard arrangement computed according to 


the formula 
: tf 32 1 
--| +—| (25a) 
p 3 Pe=—90 Pe=0 


This is far greater than the observed value and 
leads to the conclusion that the orientation of 
the crystals is not absolutely at random, but 
that there is a tendency for the zinc crystals to 
form with their principal axes of symmetry 
perpendicular to the axis of the cylindrical 
casting. This is in agreement with the known 
fact that in making single zinc crystals they tend 
to crystallize with the principal axis of sym- 
metry perpendicular to the length of the crystal. 

Since zinc crystals are easily cleavable, the 
determination of their orientations could be 
readily made. The tension coefficients were found 
~ € Tomlinson, Phil. Trans. 174, 1 (1883). 


7H. Rolnick, Phys. Rev. 36, 506 (1930). 
’ P. W. Bridgman, Proc. Am. Acad. 60, 305_(1925). 
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Fic. 1. Tension coefficient of zinc. 


to be independent of the secondary orientations, 
as demanded by the theory. For instance, three 
crystals of approximately 85° primary orienta- 
tion and of secondary orientations 3°, 10° and 
14° had respectively the tension coefficients 2.71, 
1.58 and 2.91X10-* which agree within the 
experimental error to be expected as indicated by 
the agreement of the experimental points with 
the empirical curve. 


CADMIUM RESULTS 


Sixteen cadmium crystals of various orienta- 
tions were measured. The resulting data were 
again compared with the theory, as given in 
Eq. (2). Fig. 2 then shows the agreement of the 
observed experimental points with the curve 
determined according to this theory. In only one 
case is the discrepancy between an experimental 
point and the curve greater than 2X10~-*. This 
agreement is satisfactory in view of the fact 
that the deflections were again at least as small 
as in the case of zinc. The constants in Eq. (2) 
had the following values: 

(a) when no correction is applied for the change in shape 
of the crystal, 
pu =( 6.08+0.33) x10, 
pss = (16.3 +0.53) X10, (26) 
(pist+psit+ pas) = (11.7 41.7 )X10™, 


(b) when the correction is applied for the change in shape,’ 


pu =( 4.81+0.34) x10, 
pss = (11.6 +0.55) x10, (27) 

(pis + psit pas) =( 6.2241.8 )K10™, 
The tension coefficient for polycrystalline cad- 
mium has not been determined and so the results 


_* These corrections are computed from the strain coeffi- 
cients for cadmium given by P. W. Bridgman, Proc. Am. 
Acad. 60, 305 (1925). 
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Fic. 2. Tension coefficient of cadmium. 


of this experiment cannot be compared with it. 

Inasmuch as cadmium crystals do not cleave, 
the determination of the crystal orientation was 
a matter of some difficulty. Various methods were 
tried. One was to observe under a microscope 
the angles made by the glide planes when the 
crystal was stretched.” Another was to project 
the reflection pattern on a screen, supporting the 
crystal in a spherical flask of benzene which 
acted as a lens." However, neither of these 
methods seemed to work easily and in the end 
there was adopted P. W. Bridgman’s method” 
of observing the reflection pattern with the eye 
and recording it by marking on a wooden sphere 
(through whose center the crystal was rigidly 
attached) the points at which a mirror held 
parallel to the reflection planes of the crystal 
would touch the sphere. Cadmium crystals then 
showed hexagonal symmetry, their reflection 
pattern consisting of twelve dots in each hemi- 
sphere, six near the pole and six near the equator. 
Regular reflections from the glass interfered with 
the crystal reflections if the crystals were ob- 
served before the removal from their glass 
jacket. The crystal reflections could be notice- 
ably increased by etching them in dilute HCl, 
where one part of water was added to two parts 
of acid. The resulting orientations are probably 
good to about 3°. 

I am much indebted to Harvard University for 
permission to carry out these experiments in the 
Research Laboratory of Physics and to Pro- 
fessor P. W. Bridgman for advice in connection 
both with the experiments themselves and their 
interpretation. 


1 Bruce Chalmers, Phil. Mag. 14, 612 (1932). 


1 E. Jakowlewa, Physik. Zeits. Sowjetunion 3,429 (1933). 
2% P, W. Bridgman, Proc. Am. Acad. 60, 305 (1925). 
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Paper I. Fundamental Principles and Transformations Between Accelerated Systems 
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A new approach to the relativity theory, suggested by 
the theory of E. A. Milne, is developed. This approach, 
like Milne’s, dispenses with the concepts of measuring 
rods of undefinable rigidity and clocks of undefinable 
periodicity. A new category of equivalent relatively 
accelerated reference systems with Euclidean geometry 
and constant light-velocity is described, and the space- 


time transformation for such systems is developed. It is 
shown that in an effectively empty world Einstein's 
assumption of an invarient physical interval and an 
absolute four-dimensional space-time is in contradiction 
with the underlying principle of the relativity of motion, 
and therefore either the one or the other must be aban- 
doned. 





HE fundamental assumption underlying 

Einstein’s theory of relativity is that the 
physical interval between two nearby events 
(the square of the element of measured distance 
minus the square of the product of the velocity 
of light by the element of measured time) is an 
invariant having the same value for all reference 
systems. This assumption is a natural inference 
derived from the Minkowskian complex four- 
dimensional space-time representation of the 
Lorentz transformation, and has led to cosmo- 
logical predictions which have been verified by 
observation. Nevertheless, the author of the 
present paper believes that Einstein’s postulate is 
too restricted to include all possible motions of 
material particles. In this paper he will present 
an alternative theory, and will give reasons for 
believing that it, rather than Einstein’s theory, 
represents the proper formulation of relativity 
in an effectively empty world. 

The present investigation was prompted by 
the perusal of a recent book by E. A. Milne,' to 
whom the writer wishes to make due acknowledg- 
ment. In this important work Milne offers an 
approach to the relativity theory which avoids 
the undefinable concepts of rigid measuring rods 
and periodic clocks. In spite of their great 
advantages, the writer believes that Milne’s 
methods are faulty in certain respects, particu- 
larly his definition of equivalence, in that it 
implicitly involves synchronism as well, and 
his apparent belief that physical geometry is 
conventional. The fundamental principles pro- 


1E. A. Milne, Relativity, Gravitation and World-Structure 
(Clarendon Press, Oxford, 1935). 


pounded here, while suggested by Milne’s treat- 
ment, differ from his in many essential particu- 
lars, and the space-time transformations for 
relatively accelerated reference systems are be- 
lieved to be altogether new. 

This contribution is divided into four parts. 
In Part 1 fundamental methods will be outlined 
and the principle of relativity will be stated in 
its general form for the effectively empty world 
in which we are interested; in Part 2 applications 
will be made to one-dimensional systems; in Part 
3 the special theory for a three-dimensional space 
will be shown to follow immediately from the 
fundamental principles; and in Part 4 the new 
transformations for relatively accelerated three- 
dimensional reference systems will be developed 
and contrasted with Einstein’s theory. It is the 
author’s intention to follow this paper shortly 
by another in which the transformation of the 
electromagnetic field between accelerated sys- 
tems will be developed and the necessary revision 
of the fifth or force equation of electromagnetic 
theory, which is demanded by the new relativity 
principle, will be obtained. In the conclusion to 
this paper some qualitative comments on the 
motion of an electron, which are expected to be 
developed quantitatively in the succeeding 
paper, will be presented. 


ParRT I. FUNDAMENTAL PRINCIPLES 


To emphasize the fact that a single observer's 
measurements are confined to the single point 
occupied by himself, we shall designate such an 
observer as a particle-observer. Each observer is 
supposed to possess a temporal intuition, that is 
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to say, if two events E,; and E, occur at himself, 
he can judge without ambiguity whether £, 
takes place before £;, simultaneously with £;, or 
after E;. We shall provide each particle-observer 
with a device for assigning numbers 7), tz --- to 
events occurring at himself in such a way that, 
if event EZ, occurs simultaneously with £;, the 
numbers tz and 7 assigned to the respective 
events are the same, whereas, if EZ, occurs after 
E;, then te>7, and vice versa. This device, 
which may be quite arbitrary in all other 
characteristics than the one specified, we shall 
call a clock, and we shall name r+ the local time 
of the particle-observer under consideration. 

Next we shall adopt certain conventions which 
will enable a particle-observer P to employ 
light-signals, timed by his clock, in such a way as 
to describe quantitatively the motion of any 
moving-element M. Let P dispatch a light-signal 
to M at time 7. On arrival at M the signal is 
immediately sent back to P, whom it reaches at 
time 73. Choosing an arbitrary constant ¢ (a 
constant whose numerical value, once chosen, 
remains the same for subsequent repetitions of 
the experiment) we define the distance r2 of M 
from P when the signal reaches M by 


Y2=}c(t3— 171), (1) 


and we define P’s value of the time t at which 
the signal reaches M by 


le=3(t3+71). (2) 


Since re/(t2— 71) =re/(t2—te) =c, the constant c 
represents the velocity of the light-signal in 
terms of the conventions adopted for measuring 
distance and time at a remote point. 

It should be noted that both rz and & are 
computed by P from coincidences occurring at 
himself. The first represents P’s estimation of 
the distance of M and the second P's estimation 
of the time at M at which the signal reaches M. 
We shall call & the extended time of P at M. Ifa 
second particle-observer is attached to M, the 
local time of the second observer when the signal 
reaches him may be quite different from P's 
extended time f:, and his estimation of the 
distance of P may not agree at all with ro. The 
notation used designates local time measured at 
a particle-observer by the Greek letter 7, the 


computed time at a distant point being denoted 
by the Italic letter ¢. The only measurements 
made are carried out at the particle-observer 
under consideration, and no yardstick of unde- 
fined rigidity nor clock of undefined periodicity 
is assumed. 

Evidently each one of two particle-observers 
P and P’ constitutes a moving-element in the 
experience of the other. Thus P, acting as 
observer, may describe the motion of P’ relative 
to himself, or P’, as observer, may describe the 
motion of P. We shall designate by letters 
without primes local times measured by P or 
quantities computed therefrom, and by corre- 
sponding letters with primes local times meas- 
ured by P’ or quantities computed from these 
times. We attribute to light-signals dispatched 
from one particle-observer to another the fol- 
lowing property: Jf two light-signals are sent from 
one particle-observer to another, the light-signal 
which is dispatched later from the one will be 
received later by the other. This fundamental 
principle underlies all the theory to be developed. 
In effect it is equivalent to limiting our consider- 
ation to particle-observers with relative velocities 
less than the velocity of light. 

Now suppose that a light-signal is dispatched 
from P toward P’ at time 7 and is received by 
the latter at time 72’. Let a second light-signal 
be dispatched from P’ toward P at a time 7)’ 
earlier than 72’ and be received by P at a time 
re later than 7, the time 7,’ being so chosen that 
Tt?’ — 7), =T2—7). Then we say that 7, and 7)’ are 
corresponding times. Evidently this condition can 
always be fulfilled, for if t2’—7;'>re—7, the 
light-signal from P’ can be replaced by one sent 
a little later, which will increase both 7,’ and re 
by virtue of the principle stated in the last 
paragraph, making 7t2’—7;’ smaller and tz2—1 
larger. The pair of light-signals under discussion 
are illustrated by the lower solid lines in Fig. 1, 
the time being plotted vertically and the sepa- 
ration of P and P’ horizontally. 

The statement that 7; and 7,’ are correspond- 
ing times does not necessarily imply that rz and 
t;’ are corresponding times also, for, if the signals 
received by P’ and P at the times 72’ and rz are 
immediately returned and reach P and P’ at 
the times 7; and 73’ respectively, the fact that 


lad ’ 
Te —~ Ti 


T2— 7; does not of necessity lead to the 
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equality of 73’—72’ and r3—72. In conformity 
with our present notation we shall always 
designate corresponding times by identical sub- 
scripts. 

Next consider a pair of corresponding times 
titAr, and 7’+Ar;’ defined by the pair of 
light-signals represented by broken lines on the 
figure, which are dispatched from P and P’ at 
times respectively Ar; and Ar,’ later than the 
signals sent at the corresponding times 7, and 
7;,, and received at times Arz and A7,’ later. As 
the times of dispatch correspond, Are’ —Ar;’ 
=Arz—Ar. Now, if Ar:’=Ar, and hence Ar’ 
=Are, whatever Ar, may be, we say that the clocks 
of P’ and P are equivalent, or that the two particle- 
observers are equivalent. If, in addition, the clocks 
of the two particle-observers are set so that 
7;,=7:, and therefore all corresponding times 
are identical, the clocks of the two observers are 
said to be synchronized. In future we shall deal 
only with particle-observers who are equivalent, 
and, when we are concerned with two particle- 
observers alone, we shall generally suppose their 
equivalent clocks to be synchronized. 

It follows from the definition of equivalence 
that all pairs of corresponding times at two 
equivalent particle-observers differ by the same 
amount. Therefore, if 7; and 7; are corresponding 
times, re’ and rz are also. Conversely, if all 
pairs of corresponding times differ by the same 
amount, the particle-observers are equivalent. 
If the clocks of two particle-observers are 
synchronous, corresponding times are identical. 
Hence synchronism implies equivalence, although 
equivalence may exist without synchronism. 
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Let P and P’ (Fig. 1) be equivalent but not 
necessarily synchronous. Let the signals dis- 
patched at the corresponding times 7; and 7,’ 
and received at the corresponding times r2 and 
te’ be immediately returned toward the particle- 
observers from whom they originated, reaching 
the latter at the corresponding times 7; and 7;’. 
In this case we may say that the signals interlock, 
the signal dispatched from P at the time 7, 
being received by P’ at the time 72’ and immedi- 
ately returned to P whom it reaches at the time 
t3. Evidently rz is some function of 7, which 
could be obtained empirically by observing the 
values of r2 corresponding to different values of 
71. Now, if 7; becomes t2, 72 becomes 73. So 7; 
must be the same function of 72 as 72 is of 7. 

If we are given the law of motion of P’ 
relative to P, that is, if we know 72 as a function 
of P’s extended time ¢ at P’, we can express 1; 
as a function of 7, by (1) and (2). Let this 
functional relation be 73;= F(7:). Since we must 
have 


t3=f(T2), t2=f(11), (3) 


it follows that our problem is to find the function 
f such that 
f(f(m1)) = F(1). (4) 

Not only are relations (3) necessary for equiva- 
lence; they are also sufficient. For all we need do 
is to assign the values 7: +k, re+k, 3+ --> to 
the times 7’, 72’, 73’ at which the various 
signals in Fig. 1 are dispatched from P’, where k 
is a constant. Then the clocks of the two particle- 
observers are equivalent. If k=0 they are 
synchronous as well. 

An alternative condition for equivalence arises 
from the fact that Eqs. (3) imply that dr3/dr. 
is the same function of 73 and 72 as dre/d7; is of 
vz and 7;. Let this function be denoted by g. Then 


dr3/dt2=g(Ts, 72), dte/dr,=g(t2, 71). (5) 


But, from the equation of motion of P’ relative 
to P, we get dr3/dr,;=G(rs, 71), where G is a 
known function. Consequently 


g( 7s, T2)g(72, T1) =G( 13, 71). (6) 


If we can find a function g satisfying (6), we can 
then integrate (5) and determine the common 
constant of integration so as to satisfy (4), thus 
obtaining the relations (3). Consequently this 
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condition, like the preceding one, is both neces- 
sary and sufficient for equivalence. 

Consider two equivalent particle-observers P 
and P’. From (1) the distance of P’ from P at 
time 72’ at P’ is 

— - 

re=9C(T3— 71), (/) 

whereas the distance of P from P’ at the corre- 
sponding time 7 is 


re’ = }4c(73' — 71’). (8) 


But 73’ — Te’ =73— Te as the observers are equiva- 
lent. Hence re’=re, that is, the two distances are 
the same at corresponding times. As regards the 
velocity v2 of P’ relative to P, taken as positive 
if P’ is receding from and negative if P’ is 
approaching P, we find from (1) and (2), 


dra = dr3/dr,;—1 


v2 =— =c— . (9) 
dts dr3 dr,+1 
whereas the velocity v2’ of P relative to P’ is 
dr,’ dr;' dr,’ —1 
v= (10) 





= =C ° 
dts’ dr;' dr;'+1 


As the two particle-observers are equivalent, 
dr3'=drz if dr;’=d7. Therefore ve’ =ve, that is, 
the two velocities are the same at corresponding 
times. Evidently the conclusions reached here 
hold also for accelerations or for higher deriva- 
tives with respect to the time. 

As dr3/dr7, is necessarily positive, we see from 
(9) that v2 can never have an absolute magnitude 
greater than c. For, as drs/dr; increases from 0 
to ©, v increases monotonically from —c to c. 
If we solve Eq. (9) for drs/dr7; we get 


dr3/dr,;= (1+ 2) /(1—Be), B2=2/c, (11) 


a relation we shall find useful later. 

If P and P’ are synchronous as well as equiva- 
lent, corresponding times are identical, and 7’, 
v’, etc., are the same functions of the extended 
time ¢’ of P’ at P as r, v, etc., are of the extended 
time ¢ of P at P’. 

To pass from a pair of equivalent particle- 
observers to a group of such observers, it is 
necessary first to consider three. So let us 
introduce in addition to the two equivalent 
particle-observers P and P’ a third particle- 


~ 


observer P’’. The readings of the clocks of P 
and P’ are fixed to within an arbitrary additive 
constant by the condition that they be equiva- 
lent. If P’’ has a specified motion in terms of P’s 
already assigned time scale such that (4) or (6) 
is satisfied, P’’ can be furnished a clock equiva- 
lent to that of P. Indeed, if P’’s motion relative 
to P is of the same type as P’’s, then the fact 
that P’ is equivalent to P is sufficient to insure 
that P” is also equivalent to P. The equivalence 
of both P’ and P” with P, however, does not 
necessarily imply equivalence with each other, 
which calls for separate investigation in each 
individual case. 

In a space of more than one dimension, motion 
cannot be completely defined by reference to a 
single observer, for, in addition to motion along 
the line of sight, an angular motion about the 
observer may exist. We have recourse, then, to a 
reference system, which is defined as a dense 
assemblage of particle-observers filling all space, 
such that each particle-observer is synchronous 
with and at rest relative to every other particle- 
observer. Let P and P’ be two equivalent 
particle-observers not relatively at rest with each 
of whom a reference system may be associated. 
If these two reference systems have the same 
geometry with respect to P and P’, respectively, 
they are said to be equivalent. If, in addition, we 
may take as P and P’ any pair of particle- 
observers in the two reference systems, the 
reference systems are homogeneous. In this case 
each particle-observer in the one reference system 
is equivalent to every particle-observer in the 
other. The Euclidean inertial systems of the 
special relativity theory are equivalent and 
homogeneous, but the Euclidean reference sys- 
tems with constant relative accelerations, the 
discovery of which is reported in Part 4 of this 
paper, are equivalent but not homogeneous. 

Insofar as electromagnetic theory is con- 
cerned we are interested in an effectively empty 
world. The philosophy underlying the relativity 
principle for such a world is that no preferred 
reference system exists in nature. Hence it is 
impossible to avoid the conclusion that the laws 
of physics must be identical relative to all equivalent 
reference systems, that is, reference systems with the 
same geometry and the same constant light velocity 
associated with equivalent particle-observers. This 
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is our statement? of the principle of relativity for 
an empty world. 

In particular all equivalent reference systems 
with Euclidean geometry and the same constant 
light velocity must be physically indistinguish- 
able. The inertial systems of the special relativity 
constitute such a category, which possesses the 
incidental property that the physical interval 
dx*?+dy?+dz?—c*df# is an invariant. If this were 
the only group of equivalent reference systems 
with Euclidean geometries and constant light 
velocities, there would be no need for a restate- 
ment of the principle of relativity. The signifi- 
cance of the present contribution lies in the 
discovery of a new category of reference systems 
with Euclidean geometries and constant light 
velocities which have constant relative accelera- 
tions (in the relativity sense) and for which the 
physical interval is mot an invariant. In all 
probability there are many other such categories 
as yet unsuspected. 


PART 2. ONE-DIMENSIONAL REFERENCE 
SYSTEMS 


As a space of one dimension has no geometry, 
it is much simpler to treat than a space of three 
dimensions. Consequently we shall confine our 
attention in this Part to equivalent particle- 
observers and equivalent reference systems in 
relative motion in a space of one dimension. 
First we shall consider a one-dimensional refer- 
ence system. 


Reference system 


Let P and P’ be two synchronous particle- 
observers such that P’ is at rest relative to P. 
Then the distance 72 of P’ from P at the time f. 
is not a function of t. Hence re=r (a constant). 
Putting }c(73— 171) for r2, rs —71=2r/c, and Eqs. 
(3) are 


T3=T2+17/C, t= +9r/c. (12) 
Since the distance re’ of P from P’ as computed 
by P’ must be the same function of tf’ as re is of 


? The principle of relativity was stated in almost identical 
terms in the author’s Introduction to Electrodynamics, 
published in 1922, but the full significance of the statement 
was not realized at that time. 
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lo, ro’ =re=r and P is at rest relative to P’. 
Moreover, as the clocks of the two particle- 
observers are synchronous, r.’=7,, the identical 
subscripts indicating corresponding times. As 


3(t3+ 11) = t2= Te’ (13) 


h= 


from (12), the extended time of the one particle- 
observer coincides with the local time of the other. 

If we introduce a third particle-observer P” 
at rest relative to P it is easily proved that P” 
is at rest relative to P’ and may be synchronized 
simultaneously with P and P’. Furthermore the 
addition law of distances is readily obtained. 
This law states that the distance of P” from P 
as calculated by P is equal to the distance of P’ 
from P as calculated by P plus the distance of P” 
from P’ as calculated by P’. 

As the local time of an event at P’ or P”’ is 
identical with P’s extended time of the event, 
there is no need of distinguishing between the 
local time of an event at one of the particle- 
observers and the extended time of the occur- 
rence of that event in the experience of one of 
the other particle-observers. We may time dis- 
tant events at P’ or P”’ by means of the extended 
time ¢ of P, secure in the knowledge that the 
local time of the event is the same as P’s extended 
time of the event. Furthermore, as the distance 
between P’ and P” as computed by either of 
them is the same as the excess of the distance of 
P” from P over that of P’ from P as calculated 
by P we may introduce a coordinate system 
with P as origin and employ only distances as 
computed by P. The aggregate of these distances 
we shall call the extended space of P. 

It follows from the above that we can adjoin 
to any particle-observer P a dense linear assem- 
blage of particle-observers P’, P’’, P’”’ --- at 
rest relative to P and synchronous with him. 
Each one of these particle-observers is at rest 
relative to every other, and synchronous with 
every other. The aggregate of particle-observers, 
therefore, forms a reference system. As all time 
and space measurements made in the extended 
time and space of P are identical with those 
made in the local time of the particle-observer 
concerned, we may refer without ambiguity to 
the extended time and space of P as the time 
and space of the reference system. 








ted 








NEW RELATIVITY 259 


Constant relative velocity 

Next consider two equivalent particle-ob- 
servers P and P’ the second of which has a 
constant velocity relative to the first. Then 
rg=0(tz—to) where fo is a constant. In accord 
with the conventions adopted in Part 1, fr is 
essentially positive, and the velocity v is equal 
to a positive constant for 4>¢» (particles sepa- 
rating) and to the same constant with the 
opposite sign for t<fo (particles approaching). 
1(r3+71) for fe, 


Putting $c(r3—71) for re and 
this equation becomes 


(r3— to) (71 —to) = (14+8) (1-8), B=v/c, (14) 


from which it follows that Eqs. (3) take the form 


T3—lbo 1+, , To—lo 1+ } 
-( ). -( - y (15) 
T2—lo 1-8 Ti: —lbo 1-8, 


which may be written as 


T2—lo re lo 
ee = T3—lbo— = 11—lbo+—. (16) 
(1—?)! c c 
Hence P’s extended time of the event 7’ is 
given by 
T2—lto To’ —ty’ 
to—to=3( 73+ 11)-!2=——__ 3 —— , az) 
(1—6*)! (1—,*)! 
where fo’ =/,) if P and P’ are synchronous. 

In addition to the two equivalent particle- 
observers P and P’ moving with constant relative 
velocity we shall now introduce a third particle- 
observer P’’ moving with constant velocity vp: 
relative to P and therefore equivalent to P. 
We shall show that the velocity vp.’ of P”’ 
relative to P’ is constant and therefore that P’”’ 
is equivalent to P’ as well as to P. Also we 
shall obtain the addition law of velocity. In 
order to make our notation consistent through- 
out, we shall designate here the velocity of P’ 
relative to P by vp, and that of P relative to 
P’ by vp’. As shown earlier, vp’ =vp». 

Consider the interlocking signals r,;—>7;’’ and 
t;'—»7; of Fig. 2. By (11) we have 


dr;/d7,;=(1+6p-)/(1—Bp-), (18) 
and drq /dre’ =(1+B8p’)/(1—Bp’), (19) 


where 8p,,-=vp./c is constant by hypothesis. 





P Pp’ Pp” 
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Now as P and P’ are equivalent, 


dro dre 1+Bp\!' dr; dr; 1+8p-\? 
dry = ie Oe 
from (15). But 
dra dra dt, dt, 1+8p1—Bp: 
joe des dr def 1—Bp 1 +p 
Comparing with (19) we have 
1—Bp,» 1—Bp  1—Bp’ 


= _ (20) 
1+B8p» 14+8p-14+6p’ 


which shows that 8p.’ or vp’ is constant. Hence 
P” is equivalent to P’ as well as to P. It is 
easily shown that the same clock which makes 
him equivalent to P makes him equivalent to P’. 

Eq. (20) is the addition law of velocity obtained 
by Einstein in 1905. It may be put in the more 
usual form 


vp = (vp: +p’) (1+vp-vp-’ c*). (21) 


Let us adjoin reference systems S and 5S’, 
respectively, to two synchronous particle-ob- 
servers P and P’ moving with constant relative 
velocity v. It can be shown very simply that each 
particle-observer in S is equivalent to every 
particle-observer in S’, and hence that the two 
reference systems are homogeneous as well as 
equivalent. Taking P and P’ as origins of axes 
fixed in their respective reference systems and 
making t=r’=0 when P’ passes P, we can 
obtain the relations between P’s and P’’s 
specifications of the position and time of the 
event Q (Fig. 3) by means of the light-signals 
indicated in the figure. From (15) 
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(1—8)'r5;=(14+8)'7r4’, (1—B8)'re’=(14+8)'71, 
as ¢t)=0 and 7,’=7r,, and 


U =3(ry/ +12’) 


1 t—(B/c)x 
= ———§ {}(r5+ 171) —38(75—71)} =— 


(1—s*)! (1B)! 
(22) 
x’ =hce(r4' — 72’) 
1 x—uvt 
=———_| }e(15— 11) —}0(175-+171)} =. 
(1—p*)! (1—,*)! 


These constitute the Lorentz transformation of 
the special relativity theory for one dimension. 
Constant relative acceleration 


We shall now investigate the properties of the 
linear reference systems adjoined to two syn- 
chronous particle-observers P and P’ which 
have a constant relative acceleration @ (in the 
relativity sense). The differential equation of 
motion of P’ relative to P is 


dro, dt? = (1 —v_? c?) ig, 


the integral of which is 


1+ 7ro/c?=(1+¢%to?/c?)!, (23) 


if the particle-observers meet at rest at time 

zero. Expressing 72 and f2 in terms of 7 and 1; 
as usual we find 

1 ™m1—1 t3=0/C, (24) 

and Eqs. (3) become 

1 To—1 T3=0 2c, 1 ™m1—1 T2= 2c. (25) 


First we shall obtain the addition law of 
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acceleration by introducing a third particle- 
observer P’’ who has a constant acceleration 
op relative to P, and who meets P without 
passing at the same time that P’ does. Then 
P” as well as P’ is equivalent to P. We shall 
make P” as well as P’ synchronous with P, with 
t=r’'=7’=0 at the instant of meeting. For 
uniformity of notation we shall here denote the 
constant acceleration of P’ relative to P by ¢p.. 


Considering the interlocking signals 7;—>73"’ and 
t3'—7; of Fig. 4 we have from (24) 
1/71—1/7;=¢p-'/C, 
and from (25) 
1/71—1/72’ =gp-/2c, 1/74 —1/75=@p-/2c, 


as t.=7, since P’ and P are synchronous. 
By combining, 
1 'T9! —1 4 = (op —gp')/C, 

which shows that P’”’ has the constant accelera- 
tion dp.’ =dp— dp: relative to P’ and therefore 
is equivalent to P’. It can be shown very simply 
that the same clock which makes P”’ synchronous 
with P makes him synchronous with P’. The 
addition law of acceleration, is, then, 


op =or tor. (26) 


We can adjoin to each of the synchronous 
particle-observers P and P’ moving with con- 
stant relative acceleration @ a dense linear 
assemblage of synchronous particle-observers 
relatively at rest. The two reference systems S 
and S’ so formed are equivalent but not homo- 
geneous. To find the space-time transformations 
we shall take X and X’ axes in the direction of 
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the acceleration of P’ relative to P with P and 
P’ as origins of S and 5S’, respectively, and 
consider the signals (Fig. 5) necessary to specify 
the event Q. From (25) 


T5 Ti 
; , 


=T aes 


“14(6/2c) 75 1—(¢/2c)r; 





as tr. =7,-. So P’’s extended time of the event Q is 


=3(r4' +72’) 
3(75+171) 


> f{Ts—T1 2 € ftstn\? 
mae) 

2c 2 4c? 2 

ox 2 ¢t’ 
Ee 

2c? 4c? 


and P”’s estimation of the distance x’ of Q is 











; 


x’ =4e(74' — 72’) 


O/ts—T1\”) Of ts+71\? 
bel r 42 ( ) -“( odes ) 
2 2 2 2 


@fts—T1\)? # ftst71\? 
CECE) 
2c 2 4c? 2 


— (28) 





Note that S’ extends only from —2c?/¢+(cf)! 
to ~, 


If we put 
§=1+¢x/2c*, 


t’=1-—¢x'/2c?, T’=¢t'/2c, 


T=9t/2c, 


which amounts to taking a new origin in S at 
— 2c?/¢, and a new origin in S’ at 2c?/¢ combined 
with a change in sense of the axis, then the 
space-time transformation assumes the simpler 
form 


T’=T/(#-T?), T=T'/(¢"-T”), 


P " (29) 
=t/(P—-T*), E=t'/(¢"—T"), 

This transformation gives 
(2—T7°)(¢?—T”) =1, (30) 

and yields the invariants 
T’/?’ =T/é, (31) 

dx” —cdt” dx*—cdt 

and —_— = (32) 


A wt ass a ~ 
f-T* #7 


It is seen from (32) even in this one-dimen- 
sional case that the physical interval dx” —cdt’” 
in S’ is mot in general equal to the physical 
interval dx*—c*dé? in S. Hence the rather firm 
foundations on which the present theory rests 
are quite incompatible with the fundamental 
postulate on which Einstein’s theory is based. 

The relation between the velocity V’ =dx’/di’ 
of a moving point relative to S’ and its velocity 
V =dx/dt relative to S is given by 


T’\ s1-V'/e\} T\ s1-V/ce\} 
(2) (229 (2) (529) 
¢’ 1+-V'/c é 1+V/c 
In particular, the velocity v relative to S of a 

point fixed in S’ is given by 


B=v/c=2&T/(2+T?). (34) 


When ‘=0, then, all points in S’ are simultane- 
ously at rest in S. The acceleration relative to S 
of a particle fixed in S’ is 


f=dv/dt =(1—*)'t’¢. (35) 
A particle-observer in 5S’, therefore, has the 
constant relativity acceleration £’¢ relative to S. 
In terms of the coordinate measure of S, this is 
an acceleration 
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o:= §o/(t—T*?) (36) 


for a particle-observer in S’ at & at time ¢. At 
the instant ¢=0 when all particle observers in S’ 
are at rest relative to S, this reduces to ¢;=¢/E. 

It is of interest to note that even when the 
two reference systems S and S’ are relatively at 
rest time and space measurements do not agree 
except at the common point occupied by P and 
P’. For 

dt’ /dt=dx' /dx =1/(1+x/2c?) =1— 9x’ /2c?. 


Hence it is obvious that S and S’, while equiva- 
lent, are not homogeneous, and that the physical 
interval cannot be invariant. 


PART 3. EQUIVALENT ‘THREE-DIMENSIONAL 
REFERENCE SYSTEMS WITH CONSTANT 
RELATIVE VELOCITIES 


Although we have seen that a linear assem- 
blage of synchronous particle-observers relatively 
at rest can be associated with any particle- 
observer in a one-dimensional space, it does not 
follow that a dense three-dimensional assemblage 
of particle-observers can be associated with an 
arbitrary particle-observer P in such a way that 
each particle-observer in the group is at rest 
relative to and synchronous with every other. 
An analysis of the problem shows that at most 
we can associate three particle-observers satis- 
fying these conditions with an arbitrary observer 
P, such that no three of the four observers lie 
on the same light ray. We cannot, then, associate 
a three-dimensional reference system with an 
arbitrary particle-observer. The existence of a 
three-dimensional reference system in nature is 
a matter which must be investigated empirically. 
Furthermore, if such a reference system is found, 
its space geometry is also a matter for experi- 
mental investigation. For we have adopted a 
definite convention for the measurement of 
distance, and we can use it, once we have found 
a reference system, to determine experimentally 
the ratio of the circumference to the diameter of 
a circle, etc. The geometry of a reference system 
is, therefore, not conventional, but a matter to 
be investigated empirically by means of the 
procedure adopted for measuring distance. 

We take it therefore, as a matter of experi- 
mental knowledge, that there exists, at least in 
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the limited region occupied by the solar system, 
one reference system with effectively Euclidean 
geometry and constant light velocity c. We need 
not be disturbed by the fact that the deflection 
of a ray of light passing near the limb of the 
sun, or the solar red shift, may indicate a slight 
departure from the ideal conditions assumed, for 
these effects are so small that the departures 
represented by them are negligible insofar as 
the electromagnetic theory in which we are 
interested is concerned. 

Once the existence of a single reference system 
S with Euclidean geometry and constant light 
velocity is established, it becomes possible to 
prove or disprove the existence of other equiva- 
lent reference systems with Euclidean geometry 
and the same constant light velocity. In this 
part we shall outline the methods necessary to 
show that any dense three-dimensional assem- 
blage of particle-observers all of whom are 
moving with the same constant velocity v rela- 
tive to S constitute a reference system S’ 
equivalent to S, and therefore having the same 
Euclidean geometry and the same constant light 
velocity as S. As the transformation between 5S’ 
and S is merely the Lorentz transformation of 
the special relativity theory, the discussion, 
which is presented solely for the purpose of 
developing the methods to be employed in 
Part 4, will be compressed as much as possible. 

Let P be particle-observer in S chosen as 
origin of a set of axes X, Y, Z, and let P’ be an 
equivalent particle-observer moving with con- 
stant velocity v relative to S along a line parallel 
to the X axis and distant h from it. Then the 
equation of motion of P’ relative to P is 


ra? = h? +0"(te — bo)*, 
where f) is a constant, and, if we put 
a* = h?(1 — 6?) /v*, 
(a?+(r3—to)*)!+-72—to 148 








@+(i—W)'tn— 1-8 
Hence (3) becomes 
(a?+(13—to)®)'+r3—to | 1+8\' 
rn crema Prt il 


(a?+(r2—to)?)'+ 72—bo (=) 
(a?+(71—¢o)?)'+71—bo 


-te)?))- 1-8 
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which can be written in the more convenient form 


9 


T2—ly ro re 
—- = Tg —lyp —— = 71 —L + —. (39) 
(1—?)} Cc c 
Eq. (39) is identical with (16) for the one-di- 
mensional case. As P may be any particle-observer 
in S, it follows that each particle-observer in S’ 
is equivalent to every particle-observer in S. 

Now suppose that P’ is moving along the X 
axis of S. As proved in the discussion of the 
one-dimensional case we can associate with P’ a 
dense linear assemblage of particle-observers 
distributed along this axis all of which are at rest 
relative to P’ and synchronous with him. All of 
these particle-observers have the same constant 
velocity v relative to S, and, in view of what we 
have just proved, are equivalent to every 
particle-observer in S. We will take their locus 
for the X’ axis of S’. The Lorentz transformation 
(22) applies, then, to the linear assemblages of 
particle-observers constituting the X and X’ 
axes of S and SS’, respectively. 

Next consider an event P,’’ occurring at the 
point x, y, 0 at the time ¢. If the coordinates of 
P’ relative to S at time 0 were x, 0, 0, they are 
x:+vt, 0, 0 at time ¢. To determine the distance 
of the event P,”’ and the time of its occurrence 
in the experience of P’, we must send a light- 
signal from P’ so as to reach P,’’ at time ¢, and 
then send it back to P’. Let the light signal 
leave P’ at time 4; when P’ is at a distance 
from P,’’, and arrive back at P’ at time ¢; when 
P’ is at a distance 7; from P,’’. Then 


re=|x—x,—v(t—ri/c)}*+y", 
rs = |x—x,—v(t+r3/c)}?+y". 


We solve these equations for r; and rs and then 
determine #, and 4;. By means of the Lorentz 
transformations for the linear reference systems 
constituting the X and X’ axes of S and SS’, 
respectively, we find the local times 7,’ and 7;’ 
at P’ of the departure and return of the signal. 
Then from (1) and (2) we have 


t—(B/c)x 
3 (73'+7,’) =——— 7 (40) 
(1—p?)! 


t’ 





(x—x,—vt)? 
1—,? 
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Now suppose that P” represents a particle- 
observer moving relative to S with the same 
velocity v as P’. Then x=x.+v1, y= yo, where x2 
and y are constants. In this case (41) becomes 


; (xe—x,)? . 
r' =| ——+y"] , (42) 


showing that r’ does not change with the time. 
Hence all particle-observers moving relative to S 
with the same constant velocity v as P’ are at 
rest relative to P’. To be synchronous with P’ 
their local times must be equal to P’’s extended 
time ¢’. We observe from (40) that they are then 
synchronous each with each. Finally (42) repre- 
sents the Pythagorean theorem for S’. From it we 
see that we can construct a Euclidean mesh in 
S’, and therefore that the geometry of S’ is 
Euclidean. Consequently S and S’ are equivalent 
homogeneous reference systems with Euclidean 
geometries. The Lorentz space-time transforma- 
tion for the three-dimensional case under discus- 
sion is obtained immediately from (40) and (41). 


PART 4. EQUIVALENT THREE-DIMENSIONAL 
REFERENCE SYSTEMS WITH CONSTANT 
RELATIVE ACCELERATIONS 


We have shown that a dense assemblage of 
particle-observers all moving with the same con- 
stant velocity v relative to a given Euclidean 
reference system S may be synchronized each 
with each so as to constitute an equivalent ref- 
erence system. This category of reference sys- 
tems, however, does not comprise all reference 
systems equivalent to S which have Euclidean 
geometry and equal constant light velocity. We 
shall now show that a three-dimensional reference 
system 5S’, equivalent to S, may be adjoined to a 
particle-observer P’ moving with constant rela- 
tivity acceleration ¢ relative to S, and that the 
geometry of S’ is Euclidean. 

Take as origin of S the synchronous particle- 
observer P with whom P’ coincides when momen- 
tarily at rest in S and orient the X axis in the 
direction of the acceleration of P’ relative to S. 
As proved in Part 2, we can adjoin to P’ a linear 
reference system extending along the X axis of 
S from —2c?/¢+(c#)! to «, the particle- 
observers constituting this reference system hav- 
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ing the constant accelerations relative to S 
specified by (36). The transformation (29) applies 
to the linear assemblages of particle-observers 
constituting the X and X’ axes of S and S’, 
respectively. 

Let Q’ be a particle-observer in the linear 
reference system adjoined to P’ at a distance x’ 
from P’, which serves as the origin of the X’ 
axis in S’. We shall denote by x» the coordinate 
of Q’ in S when t=?/=0. Then the constant ac- 
celeration ¢, of Q’ relative to S is given by 


71) oi = 1 + oxo ‘2c? (43) 


in accord with (36). 

Next consider an event Q,”’ occurring at x, y, 0 
at the time ¢. We dispatch a light signal from Q’ 
at a time ¢, so chosen that the signal will reach 
Q,’ at time ¢, and then send the signal back to Q’ 
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whom it reaches at time /;. Let 7; be the distance 
of Q,’’ from the position x, of Q’ at time ¢,, and r, 
the distance of Q,’’ from the position x; of Q’ 
at time /;. Then 


1 o: n\-7? 
1+ (21a) =| 14 (:- ) |. 
c c c 
$1 or rs\°7 
14 (xx) =| 1+ (+4 )]. 
me a ¢ 


where 7;?=(x—x)?+y", 13?=(x—4x3)*+y". (45) 


From these equations we must find 7 and 73, 
and then 4; =t—r,/c and tj=t+73/c as functions 
of x, y, xo, t. If we put pcos @=1+¢x/2c, 
psin 06=y/2c, &b=1+4x0/2c, T=¢t/2c, after 
expressing ¢; in terms of ¢ by (43), we find as 
the result of a laborious algebraic calculation 


ot; {Eo P+ 2(p°—T *— Sop | cos s 8) —(2e cos 6—£)[ ( (o? — 1%)? — 2(o*— T”) )Eop cos 6+"? |! 


2¢ (2p cos 6— 
ots Thk +2 (o*— T*— Sop cos 8)} + (2p cos 6— 


2c (2p cos @— 


£o) 2472 
$0 ~ ad )?— 2(p?— T*)Eop cos 6+ £97»? |} 


To find the local times 7,’ and 7;’ at Q’ corresponding to ¢, and ¢; we have from (29) 


or’ /2c= —c/bt+(c2/ pt? +1/k,?)!. 


Thus we find 


ori _&t— ~C( p° ?— T*)? — 2(6? “a )Eap COs 9+" 2p? |! 


2c 3 


2)? 2(p?—T? )Sop COs o+8 op}? 


oT3' foT+[(e?—-T 


2c g —T?) 
Finally, we get for Q’’s time #’ and distance r’ of the event Q,”, 
ot! /2c=(o/4c)(13' + 71') =T/(e?—T?), (46) 
or [ (p?— T?)?—2(p?-—T? )Eop cos 6+ £9?p? }! , 
—=—(13'— 1;') =—____ ~ (47) 
2c? 4c Eo(p?: —T?) 


Now suppose that Q” is a particle-observer 
moving along a radial line in the X Y plane of S 
drawn from the point O at —2c*/¢, 0, 0 at an 
angle @ with the direction of P’’s or Q’’s motion, 
as shown in Fig. 6. Let Q” be at rest in S at the 
same instant that P’ and Q’ are, and let Q”’ have 
a constant acceleration ¢2 which is the same func- 





tion of his distance from O as that of P’ or Q’ is. 
Then, if we denote by 7 the distance of Q” from 
an origin in Son the line OQ” at the same distance 
from O as P, 

$/¢2=1+ 70/2’, (48) 


where fo is the value of r at ‘=0, and 
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Fic. 6. 


1+ (¢2/c?)(r—ro) = (1+-¢2F/c*?)!. (49) 


Eliminating ¢: from (49) by (48) and putting po 
for the value of p at t=0, we find 


p?>—T*=ppo. (50) 


Combining (50) with (47) we find that we can 
eliminate both p and T, getting, 
or’ [ po? — 2£Eopo cos O+E,? }}. - 
a (51) 
2c Eopo 


Therefore the distance r' of Q” from Q’ as measured 
by Q’ does not change as the motion progresses. 
The aggregate of uniformly accelerated linear 
reference systems radiating from the point O forms 
a three-dimensional reference system S' each 
particle-observer of which is permanently ai rest 
relative to every other. Moreover, comparison of 
(46) with (29) shows that Q’’s extended time at 
Q” is identical with Q’’’s local time. Hence each 
particle-observer in S' is synchronous with every 
other. 

It remains to show that the geometry of S’ is 
Euclidean. Since the particle-observers in S’ 
form a rigid aggregate in their own distance 
measure, it is sufficient to investigate their 
geometry at the instant ‘=//=0 when S’ is at 
rest relative to S. As the projection of S’ on S is 
symmetric about the X axis of S, we can simplify 
the analysis by first investigating the geometry 
of a section of S’ lying in the X Y plane of S. 

First we shall find the projections on S of all 
straight lines in S’. Put Ro=2c?/¢+1r0, Xo 
=2c?/¢+x 9. Then Ro and Xo are the distances 
of Q” and Q’ from the singular point O at t=0, 
and the Pythagorean theorem (51) becomes 

4c* [Ro?—2RoXo cos 0+ X 7 }! 

‘= —————_————. (52) 
¢ RoXo 





r 


If Q’ and Q” are neighboring points, this gives for 


the element of distance 


4ct (dRP+Rede |} 
dr’ =— — ' 


? R,? 


(53) 


A straight line is defined by 6 fdr’ =0 between 
fixed limits. Minimizing the integral by the usual 
methods we get 


dR, /d?+Ro=0, (54) 
the complete integral of which is 


(Ro cos 0—A)*?+(Ro sin 6—B)?=A?+B?, (55) 


where A and B are constants of integration. 
Hence all circles through O in S are straight lines wn 
S’ at t'=0. In particular all straight lines in S 
radiating from the point O are straight lines in S’. 
The point O in S is the point at infinity of S’, 
and therefore there are an infinite number of 
straight lines joining O to any point in S’. 

Next we shall find the projections on S of all 
spheres in S’. The equation in S of a sphere 
about Q’ in S’ with radius r’ is the relation be- 
tween Ro and @ given by (52) when Xo and r’ are 
held constant. Rearranging this equation, 


Xo 2 


Ry cos @———— 
1-— (¢?r’ 4c*)?X °? 


(g?r’ /4c*)2X of 
+ R,? sin? = , (56) 
11 —(@?r’ /4c*)2X 67}? 


which is a sphere about a point on the OX axis 
distant 

(g?r’ /4c*#)?X 3 _ 

— (57) 

1 —(¢?r' /4c*)?X¢? 

from Q’. The projections on S of the radii of this 

sphere in S’ are, of course, circles through O 

and Q’. 

Now we shall show that angles are preserved in 
passing from S’ to S. First we express the Eq. 
(56) of a sphere in S’ about Q’ as center in terms 
of polar coordinates p, x with Q’ as origin and OX 
as polar axis, getting, 


£? = (¢"r' /Act)2X 2(p2-+2X op cos x+X02). (58) 


To obtain an element of the circumference in the 
XY plane of S we must differentiate (58) 
holding r’ and Xo constant, and substitute in 
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(53), which becomes 


4ct  (dp*+ p*dx’*)! 
d? =— - saarpeeoninnstitanetarie (59) 


$2 pe+2X op cos xX? 


in our new coordinates. Dividing dr’ by the radius 
r’ of the sphere we find the angle dx’ subtended 
in S’. Thus: 
dx 
dy’ =—_—____—_—_—__. (60) 
1+(p X 0) COS x 


As p approaches zero, whatever Xo and x may be, 
dx’ approaches dx. Hence angles are invariant for 
the transformation from S’ to S. 

It is easily shown that the circles through Q’ 
and O in S, which are the projections on S of 
straight lines radiating from Q’ in 5S’, intersect 
orthogonally the projections (56) on S of all 
spheres in S’ with Q’ as center. 

We are ready now to construct a Euclidean 


mesh in SS’. First we show that all circles in S 


tangent to OX at O are the projections of 
straight lines in S’ parallel to the X’ axis. The 
equations of such circles are 

R°? cos” 6+ (Ro sin 6—B)?=B?, (61) 
and the perpendiculars to the X’ axis in S’ are 
the circles 


(Ro cos 0—A)?+ R,? sin? 6=A?. (62) 


These two families of circles intersect orthogon- 
ally, showing that in S’ (61) is perpendicular to a 
straight line at right angles to the X’ axis. 

All that remains is to calculate the lengths of 
the sides of a rectangle in S’ bounded by the 
straight lines (61) and (62). If Q;’ is a point in S’ 
on the perpendicular to the X’ axis through Q’, 
then the coordinates Ro, 6 of Q;’ must satisfy the 
relation Ro=Xpo cos 6, where Xo, O are the co- 
ordinates of Q’. Hence the distance of Q,’ from 
Q’, obtained from (52), is 

4c‘sin@ 4c‘tan 6 
r= — —_ = —_ —_. (63) 
@ Ro @g# Xo 


This is identical with (61), showing that all 
straight lines in S’ represented by (61) are equi- 
distant from the X’ axis. Incidentally, for a 
given Xo, r’= 2% when 0=7/2. 
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Next consider two perpendiculars in S’ to the 
X’ axis passing through the points whose co- 
ordinates are X;, O and X2, O, respectively. The 
equations of their projections in S, obtained from 
(62), are Ro=X; cos 6 and Ro=X-_ cos @. Let Q,' 
and Q,’ be the intersections of these curves with 
(61). Then the coordinates R), 6, of Q;’ and Rs, 6, 
of Q,’ satisfy the relations 

R,=X, cos 6;=2B sin 6,, 
R.= X_ cos 62=2B sin Oe. 


The distance of Q,’ from Q,’ in S’ is, in accord 
with the Pythagorean theorem (52), 


P 4c* [R2?—2R2R, cos (02.—6,;)+R:? }! 


r 
¢? R2R, 
(64) 
4c* Xo—X;, 
go XX. 


which 1s independent of B. 
We have now constructed a Euclidean mesh in 
any plane of S passing through the X’ axis. Lastly 
we have to consider measurements involving a 
change in azimuth y about the X axis. Let Q,’ 
and Q,’ be two points with coordinates Ro, 0, ¥ 
and Ro, 6, ¥+dy, respectively. The angle sub- 
tended at O by radii vectors to Q,’ and Q,’ is 
sin @dy. Hence (53) gives for the distance be- 
tween them in S’, 
4c‘ sin 6dy 


dr’ =- (65) 
? Ro 


Their common perpendicular distance in S’ from 
the X’ axis is specified by (63). Dividing (65) 
by (63) we obtain for the difference in azimuth 
in S 

dy’ =dr'/r' =dy. (66) 


So, if we measure y’ and y from coincident planes 
through the X axis, we have y’=y. 

We have completed the proof that 5S’ is a 
Euclidean reference system with constant light 
velocity c, constituted of synchronous particle- 
observers relatively at rest. We shall now as- 
semble the space-time transformation between S 
and S’. First we note from (63) that the distance 
in S’, measured along the perpendicular to the X’ 
axis, of points on the radial line OQ” of Fig. 6 
becomes less and less as Ry increases. From (29) 
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we can write (63) in the form 
r’ =X’ tan 0 (67) 


at ‘=/'=0, if X’ is the distance of Q’ (Fig. 6) 
from an origin O’ at a distance 2c?/¢ to the right 
of the former origin P’. All radial lines diverging 
from O in S converge at O’ in S’. Just as O is the 
point at infinity of S’, so O’ is the point at in- 
finity of S. 

Now, taking O’ as origin of a set of spherical 
coordinates R’, 6’, y’ in S’, where the polar angle 
@’ is measured from the negative direction of the 
X axis, we have at once from (67) that @’=8. 
We have already shown that y= y, although, if 
we wish to make both sets of coordinates right- 
handed we must measure y’ in the opposite sense 
to y and write ¥y’=—y. The relations between 
R’, R, t' and ¢ are given by (29). So, if we put 


p=oR/2¢, p'=oR'/22, T=¢t/2c, T’=ot'/2c, 


we have for the complete space-time trans- 
formation 


T’'=T/(e-T?), T=T'/(p"-T"), 
p' =p/(p?—T*), 
6’ =0, 0=0', 

y=—-y¥, y=-y. 


The differential invariant of this transforma- 
tion is 


[1/(p?— T*) ]\dp?+(p?— T*)(d@°+sin? 6dy*) —dT?} 
a ct (p?— T”) ]\idp” 
+(p"—T”)(d6"+sin? 6’dy*)—dT"}. (69) 


The physical interval dR?+ R?(d@*+sin* 6dy*) 
—cd between two nearby events as measured 
in S$ is mot equal to the physical interval dR” 
+R"(de”+sin? 6’dy”) —cdt” between the same 
two events as measured in S’. 

As the physical interval is not an invariant 
the present theory is incompatible with Ein- 
stein’s. But the present theory is based solely on 
the assumptions that in an effectively empty 
region (1) there exists at least one Euclidean 
reference system with constant light velocity, 
and (2) all equivalent Euclidean reference sys- 
tems with constant light velocity are physically 


indistinguishable as regards the formulation of 
the laws of nature. The first assumption is gen- 
erally admitted to represent the result of meas- 
urement; it is difficult to see how the second can 
be denied without denying the philosophy under- 
lying the whole idea of the relativity of motion. 
For equivalent particle-observers and equivalent 
reference systems have been defined in such a 
manner that two such particle-observers or two 
such reference-systems stand in precisely the 
same relation to the underlying constant light 
velocity. In fact the space-time of a reference 
system has been constructed, not from yard- 
sticks of undefinable rigidity and clocks of un- 
definable periodicity, but from the concept of a 
universal constant light velocity. Hence the con- 
clusion seems inevitable that the fundamental 
assumption of an invariable physical interval, 
which underlies Einstein's relativity, is unten- 
able. Either the postulate of an absolute four- 
dimensional space-time, or the postulate of the 
relativity of motion in an effectively empty 
world, must be abandoned. 

Let us suppose that no external electromag- 
netic field is present in S at the time t=0 at which 
the relatively accelerated reference system S’ is at 
rest relative to S. Then, clearly, no external elec- 
tromagnetic field is present in S’. Consider an 
electron at rest in both S’ and S at this instant. 
Presumably an electron will remain at rest in a 
Euclidean reference system with constant light 
velocity in the absence of an external field. 
Therefore it appears as if the electron under con- 
sideration has a choice as to whether it shall re- 
main at rest in S, S’ or another of the infinitely 
many equivalent reference systems with con- 
stant acceleration relative to S. Actually, how- 
ever, no indeterminacy exists. For the angular 
distribution relative to S of the charge of an 
electron of finite dimensions is different according 
as the electron is permanently at rest in S or S’. 
What determines this angular distribution the 
present theory does not indicate. But a given 
angular distribution specifies the reference sys- 
tem in which the electron remains at rest. 

Consider an electron permanently at rest in S’. 
Then, relative to S, the electron, starting from 
rest in a field-free space, moves away with con- 
stant acceleration and ever increasing velocity. 
Although no work is performed by external 
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forces, the kinetic energy of the electron con- 
tinually increases. No violation of the conserva- 
tion principle relative to S is involved, however. 
For, as the velocity of the electron grows, its 
linear dimensions relative to S increase, and its 
mass relative to S decreases. We have, then, a 
conversion of mass into energy. As the velocity 
of the electron approaches that of light, this 
process of conversion approaches completion. We 
have here a possible method (although probably 
not precisely that occurring in nature) of conver- 
sion of matter into energy. The converse trans- 
formation takes place during a retardation. 

It is hoped to deal with these matters quanti- 
tatively in a succeeding communication. In addi- 
tion it would seem desirable to investigate equiv- 
alent reference systems having other types of 
motion, particularly relative rotation, in the 
hope of finding a rational detailed description of 
atomic structure. 

Further consideration leads to the suspicion 


FEBRUARY 1, 1936 


PHYSICAL REVIEW 


SILBERSTEIN 


that it was not necessary to give a detailed proof 
of the fact that the geometry of the reference 
system S’ considered in this Part is Euclidean, 
For the particle-observer P’ to whom the refer- 
ence system 5S’ is adjoined is in exactly the same 
situation with respect to light-signals as is the 
particle-observer P to whom the reference system 
S is adjoined. Therefore, as our geometry is one 
based on light-signals, the geometry of a reference 
system adjoined to P’ must be identical with that 
of a reference system adjoined to P. If one is 
Euclidean, the other must be also. Nevertheless, 
most of the analysis presented would be required 
to find the space-time transformation between 
Sand 5S’. 

The author wishes to acknowledge his in- 
debtedness to his colleague Professor N. I. 
Adams, Jr. for a number of suggestions regarding 
the presentation of the subject matter of this 
paper, and to Mr. T. J. Carroll for verifying the 
algebra leading to Eqs. (46) and (47). 


VOLUME 49 


Two-Centers Solution of the Gravitational Field Equations, and the Need for a 
Reformed Theory of Matter 


LupWIK SILBERSTEIN, Toronto, Ontario, Canada 
(Received November 25, 1935) 


S has been shown by Levi-Civita! and as 
the reader may ascertain more directly in a 
perfectly straightforward way, the field equations 
outside of matter, Rx =0, are satisfied by the 
axially symmetrical line-element 


ds* = e*dx4? — e~*"[e?*(dx,?+dx-2")+x,%dx537], (1) 


where v and X, functions of x1, x2 only, satisfy, 
respectively, the Laplace equation 


1 @ dv 0° 
Vey=— eae (« + +) +——=0 (2) 
X1 OX, Ox, OX? 


and the condition (equivalent to two partial 
differential equations) 


Ov\? Ov \? Ov Ov 
an-x] (—) -(—) fintan — —dx2, (3) 
Ox, OX2 OX, 0X2 


1 Levi-Civita, Rend. Ac. Linc., Note VIII, Rome (1919), 


where d\=(0/0x;)dxi1+(AA/Ax2)dx2 is a total 
differential, namely, in virtue of (2). 

Since (2) is linear and homogeneous, the 
superposition of any integrals is again an integral 
of that equation. 

The object of this paper is to derive a solution 
of (2), (3) corresponding to two mass centers A, 
B, a field, that is, which has singularities at A 
and B only, and not (as in R. Bach’s and H. 
Weyl’s physically trivial solution®) along the 
straight segment joining these two points. 

I may mention that I have constructed such 
a solution (a stationary one) in December, 1933 
and have then communicated it to Einstein, 
pointing out, rather emphatically, that this is a 
case of a perfectly rigorous solution of his field 
equations and yet utterly inadmissible physi- 

2R. Bach and H. Weyl, Math. Zeits. XII, 134, Berlin 
1922; see especially page 141 et seq. 
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cally, so that one cannot henceforth treat 
“material particles’ as singularities of the field. 
This has, in fact, induced Einstein to attempt, 
in collaboration with N. Rosen, a new theory 
of matter.* More recently, Professor Epstein of 
Pasadena has asked me to communicate to him 
the details of my derivation of that embarrassing 
solution which I have gladly done. Now, both 
Einstein and Epstein have expressed the opinion 
that a publication in toto of the derivation of 
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this two-mass centers (or better, two singular 

points) solution is, in view of its fundamental ‘ 

implications, by all means desirable. And it will, 

therefore, be given here in all detail. (It is, 

moreover, a perfectly straightforward affair.) 
But before turning to two “‘mass centers” it 

will be well to consider briefly one such center, 

since the solution (1) then by no means reduces 

to the well-known radially symmetrical Schwarz- 

schild field. 


SINGLE Mass CENTER 


Take as a solution of (2), 


v=-—L/r, 


r=x ? +27, 


where L= W//c*, say, is a constant. Then (3) is readily integrated and the line-element (1) becomes 


ds? =e-*4! "dx? — e24!"[ el? 


or, with «;=7 sin 6, x2=r cos 0, x3=¢, 


ds? = e~®L! dx? — e2L!t[ e(—L? sin? 


Let us now turn to our main subject. 


Two Mas 


Take, as a solution of the Laplace Eq. (2), v= 
re?=xP+(xe+a)*, re2=x?+(x2—b)*, where a and 
apart of the two singular points). Then (3) gives 


z*) (dxy?+dx2*) +x,°dx;? ], 


17° (dr?4+-r2d62) +r? sin? ddg?). @) 


s CENTERS 


—L,/r,—Le2/re, where L; and Lz are constants and 
b are again constants (a+ being the x:-distance 


1 OAr ra) 1/L,;? L-,? 2x +a—b 
- | - ( +- ) Jee (5) 
2x17 0X2 OXe 4\r;' ro! ry ro 


1 OA 1 av\2 Ov 2 
and -— -|( ) -( ) |J-r-a/ 
xX) OX, Ox, OXe 


which can also be written 


Ly Lo 2 Ly Le aL, bLs 2 
sf )  2xo( , ) ri ( _ ) 
r;? ro? r;? ro® r,? re* J 


Or re] (a+x2)? 14,’ (xe—b)? 1 : xi \*% (xeta)(x2—b)x, 
era ara ne (era 
Ox, Ox, r;* r;? 2 ro! re? 2) rile ry*re* 


x? L;? L2? * x1 3 (xe+a)(xe—b)x, 
Whence A= — ~(= a )+20.L2 | ( ) = fast 1100), (6a) 
2 r,4 ret rife r,*ro 


where // is an arbitrary function of x2 alone. 
On the other hand, by (5), 


x," L;? L;? 2L Lex}? *(7;?— 1") 
A= -(=+ )+ — | ————dx2+G(x), (6b) 
a+b . 


2 r,3 ro! 


* Einstein and Rosen, Phys. Rev. 48, 73 (1935). 
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where G is an arbitrary function of x only. 
The first terms in (6a) and (6b) are identical. 


If we now put 


1 x)? L;* L2? P 
an [r+ (= +—) =y(x1,%2) (7) 
2L,L¢ 2\7r1* re! 


and choose a=} (which by a proper placing of 
the origin of x2, can always be obtained), we have 
Ou / Ox, = (%1/71372") (x12 — x2? +a”), 
Ou /OX2 = 2%12x2/7r1' re, (8) 
where 7;7=x,°+(xe+a)?, re? =x)°+(x2—a)?. 
Our result, so far, is 
N= — (x1°/2)(Li?/ri4§+-L2?/re*) +2L Ley (x1, x2), 


where yu is to be determined from the two Eqs. 
(8). 

Let us consider the iso-uv lines, »=const. 
The differential equation of any such line is, by 
(8), 

— dx_/dx,= (1/2x 1x2) (x12 — x2? +a") 
or, with u=x, du/dx,—u/xj+(x,+a?/x) =0, 
which is Euler’s well-known equation. 

Its general integral is 

u=Cx,+a?—x;?, 


where C is a constant. Thus the equation 
u=constant becomes 


(x1?+x2?—a?) /x,=C. 
Hence yu must be a function of v= (x,?>+ x2? —a?)/x, 
alone, and the form of this function is to be 


determined from (8). 
Now, 0u/dx2=p'(v)-2x2/x%1, while the second 
of Eqs. (8) is du/Ax_e= 2x,2x2/r;'re*, so that 


u’(v) = (x1/rire)*. (9) 
But this would imply that rr2/x, is (identically) 
a function of v alone. Is this the case? Yes. 
For we have 
117 ro? =x 14+ (x9? —a?)?+2x)7(x2? +a?) 

= x)4+ (xo? — a”)? +2x17(x2? — a?) +4a’x,? 
= (x? +x? —a,”)?+4a?x,?. 

Therefore 7,?r_2/x,;=v?+4a*, which is indeed a 
function of v alone. 


We now have, by (9), 
du/dv=(v?+4a?)-}, 
u=v/4a?(4a?+v?)!+A, 


where A is an arbitrary constant. Thus, recalling 


whence 


the meaning of »v, 
u—A =(1/4a?)(1 —4a?x,?/r,?r-")!, 


and if we require \=0, and therefore also n.=0, 
at infinity, we must put A = —1/4a’. 
Ultimately, therefore, we have 


v= —L,/r1;—Le/re, 


(10) 


x7/L; L?? 2L,L2 D*x?\ 3 
nae sh mel 4 
2A\r‘ ref D ry re? 
where the square root is to be taken with the 
positive sign and D=2a is the mutual distance 
of the two centers A and B. This, with (1), 
completes our solution of the field equations 
Ry=0. 

The only singularities of (10), and therefore 
of (1), are 

r3=0 and re=0, 

i.e., at A and B, while the segment AB of the 
axis of symmetry is free from singularities. In 
fact, for x;=0, the expressions (10) reduce to 


v= —L,/r,:—L2/re, 
A= —(x;? ‘2)(L7? ryit+Le, re'), 


and this has only singularities at the points A, B 
themselves. 

Thus our solution (1), (10) corresponds to a 
complete absence of matter (tensor of matter) 
everywhere, except only at the points A,B 
themselves, the ‘‘mass centers.’’ There are thus 
(unlike Bach and Weyl’s solution) no stresses 
between A and B to keep them apart like a 
stiff rod. And yet (the solution being stationary, 
invariable in time) the two points are fixed 
relatively to each other, instead of falling 
towards each other, which flagrantly contradicts 
man’s most ancient, primitive experience. 

Here, then, is an example of a perfectly 
rigorous solution of Einstein’s field equations 
which does not at all correspond to reality. If, 
therefore, these equations are to be retained, 
one cannot consider ‘‘material particles’ (mass 
points) as singularities of the field, as mentioned 
above. 
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LETTERS TO THE EDITOR 


Prompt publication of brief reports of important discoveries in physics may be secured by 
addressing them to this department. Closing dates for this department are, for the first issue of the 
month, the twentieth of the preceding month; for the second issue, the fifth of the month. The Board 
of Editors does not hold itself responsible for the opinions expressed by the correspondents. 


Communications should not in general exceed 600 words in length. 


Zeeman Effect in the Atmospheric Oxygen Bands; 
Production of a Strong Magnetic Field Over a 
Length of 80 cm 


To obtain Zeeman effects in these bands, the production 
of suitable magnetic fields in a considerable length is 
necessary. The usual magnets give fields only in spaces of 
a few centimeters length, while to get an observable 
0: absorption spectrum, a path of light at least a few 
meters in length is desired. 

The Institute’s new magnet produces fields up to 20,000 
gauss and more over a length of 80 cm. The construction 
is that of a mantle-transformer, i.e., the iron frame covers 
the coils from outside completely and reduces the scattering 
of the field to a minimum. The coils, about 80 in number, 
are water-cooled copper tubes and may carry 400 amperes, 
but a current of about 50 amperes produces a field of 
20,000 gauss in a gap of 11X80 cm. Details of the 
instrument will be given elsewhere by the constructor, 
Engineer I. Ferencz, assistant at the Electrotechnical 
Institute of our University. 

A copper tube, provided with glass windows and filled 
with tank oxygen under a pressure of 5 atmospheres, was 
placed between the (movable) magnet poles. The con- 
tinuous light of a “‘Tungsram” special lamp, 10 volts, 
7.5 amperes, passed down the 80-cm long absorption tube 
and was concentrated by a cylindrical lens on the slit 
of the 30,000 lines/inch, 6.5-meter concave grating spec- 
trograph in Paschen mounting. The path between the 
slit and the grating, the grating itself and the path between 
grating and plateholder in the region of about 7600A, 
were separated from the air of the grating room by im- 
pregnated carton tubes and boxes. Commercial carbon 





dioxide gas was let into this space constantly to keep away 
the air from the about 12-meter long path of the light, 
reducing in this way the very disturbing O2 absorption in 
the grating room itself as much as possible. 

The bands, being a '2—*Z transition, must show,—since 
the 'Z state is insensible to magnetic fields,—substantially 
Zeeman effects due to the *E ground state. At about twenty 
thousand gauss the value of Avporm, increases up to about 
1 cm™, while the distance between successive band lines 
amounts to not more than a few wave number units itself; 
in this way the Zeeman effect at such rather strong fields 
becomes perceptible by the appearance of new lines be- 
tween, or on the side of, the no-field band-lines, which are 
always present, due to insufficient expelling of air from 
the light’s path. 

The long wave-length side patterns of the P branch lines 
appear at first as shown in Fig. 1. It appears that the dis- 
placements are unequal for different lines. A detailed 
report, comparing the calculated and observed positions 
and the intensity and polarization of the patterns, will be 
given elsewhere in the near future by the writer's col- 
laborator in this work: Miss Iolanthe Zemplén. 

As the magnet seems to be unique at present in giving 
strong fields in an unusual length, the writer would be very 
indebted to anyone who might suggest problems,—espe- 
cially nonspectroscopic,—which might be solved with this 
instrument. 

R. SCHMID 


Physica! Institute of the Royal Hungarian University 
for Technical and Economical Sciences, 
Budapest, Hungary, 
January 18, 1936, 
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LETTERS TO 


The Isotopic Constitution of Rubidium, Zinc and Argon 


With the same mass spectrograph as was employed for 
the study of potassium! a similar investigation has been 
made for rubidium and zinc. No new isotopes were found. 
From measurements of the peak heights of the known 
isotopes their relative abundances were computed. By 
measuring the background in the regions where one might 
expect other isotopes it was possible to set upper limits for 
their abundance. 

For rubidium an abundance ratio Rb**/Rb*’ = 2.68 +0.02 
was found. Assuming a packing fraction of —8.2 for 
rubidium, the chemical weight computed from the above 
abundance ratio would be 85.45 in good agreement with 
the international value, 85.44. The following upper limits of 
abundances relative to Rb® can be set: Rb®, Rb*!, Rb*, 
Rb*® and Rb®, 1/100,000; Rb**, 1/60,000; Rb**, 1/12,000; 
Rb**, 1/13,000; Rb*, 1/22,000. 

In zinc the isotopes 64, 66, 68, 67 and 70 were found. 
The relative abundances agree well with Aston’s? revised 
values. With Aston’s value for the packing fraction of 
zinc, the chemical weight of zinc, when calculated, was 
found to be 65.31. This is in fair agreement with Aston’s 
value, 65.33, but is in distinct disagreement with the 
International value, 65.38. Recently, Stenvinkel and 
Svensson’ from their investigation of the band spectum of 
zinc hydride claim to have found the isotopes Zn® and 
Zn*’—both more abundant than Zn”. From the present 
work it is possible to say that Zn®, Zn®* and Zn® are 
present in zinc to less than 1 part in 40,000, 80,000 and 
60,000, respectively. These abundances compared to Zn” 
would be 1/200, 1/400 and 1/300, respectively. 

Argon was studied in another mass spectrograph very 
similar to that employed in the above investigations. The 
presence of the rare isotope, A®*, reported by Zeeman and 
de Gier* was confirmed. From a measure of the peak 
heights it was found that A*/A**=325 and A*/A®=5.1. 
In addition, the following upper !:mits can be set for the 
abundances of hypothetical isotopes relative to A*®: A*! 
and A®®, 1/10,000; A®# and A®, 1/20,000. 

I wish to express my appreciation to Professor John T. 
Tate for his valuable suggestions and criticisms of the work. 

ALFRED O. NIER 

Department of Physics, 

University of Minnesota, 
January 17, 1936. 
1 Nier, Phys. Rev. 48, 283 (1935). 
? Aston, Mass Spectra and Isotopes (1933), p. 119 


3 Stenvinkel and Svensson, Nature 135, 955 (1935). 
4Zeeman and de Gier, Proc. K. Akad. Amsterdam 37, 3 (1934) 


An Infrared Absorption Band of Heavy Water Vapor 


An infrared band of D,O vapor has been observed at 
1.52914 with a one-meter absorption tube heated to 100°C. 
Although the sample of D,O was 98 percent, as given by 
the manufacturer, the maximum absorption in the zero 
branch was only 5 percent. This is considerably less than 
the absorption of the corresponding band of H,0. 

The line separation is irregular in the positive and nega- 
tive branches, indicating a combination band, but the 
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TaBLe I. DO band. 


bu , m \ 

1.5090 6626.8 

12 
1.5118 6614.7 

7.1 
1.5134 6607.6 

10.0 
1.5157 6597.6 

11 
1.5184 6585.9 

9.9 
1.5207 6576.0 

17 
1.5248 6558.3 

18.6 
*1.5291 6539.7 

27.1 
1.5355 6512.6 

9.3 
1.5377 6503.3 

11.4 
1.5404 6491.9 

12.7 
1.5434 6479.2 

Band center 


resolving power was insufficient to obtain the fine structure 
However the center of the zero branch was definitely 
located and the following comparisons are drawn. 

Sleator and Barker! found the bands corresponding to 
the fundamental frequencies v; and v3, and Bonner? has 
calculated vz assuming the D.O molecule to be identical 
with the H,O molecule except for mass differences. 


MOLECULE v1 v v Y 
DO 1179 cm™! 2655)? 2784 cm™! 2790 
HO 1596 (3600)8 3742 3756.5)2 
De) 0.7 

- .7375 ).7387 7 7427 

>, (HzO) 375 0.73 0.7440 0.742 


From the average value of the ratio of the fundamental 
frequencies for D,O and H,O the band for DO found at 
1.52914 should correspond to an H,O band at 1.1327,. 
There is an H,O band center at 1.1340u. Mecke? gives this 
band at 1.1354 and represents it by the vibration quantum 
numbers (1, 1, 1). 

The calculated value of (1, 1, 1) vibration frequency for 
D.O by Bonner’s? equations: 


vo= Wi Vit W2V2+ W333, 
Wi =X1—- Xu —- Xw2—- 3 X13, 
W2=X2—X2—3X12— 3X23, 
Ws=X3—X33—3X13— 3X23, 


gives v»=6544 cm while the experimental value is 6540 
cm}, 

A similar calculation of v) by Mecke’s equations gives 
vo =6543 cm”, so that neither set of equations is to be 
preferred in the light of the present experimental data. 
(See Table I.) 

W. V. Norris 
University of Oregon 

H. J. UNGER 
University of Idaho, S. B. 

R. E. HOLMQuIST 
University of Oregon 

January 14, 1936. 
1W. W. Sleator and E. F. Barker, Phys. Rev. 46, 336 (1934 


2L. G. Bonner, Phys. Rev. 46, 458 (1934 
> R. Mecke, Zeits. f. Physik, 81, 313 (1933). 
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LETTERS TO 


A Lower Limit to the Normal State Eigenvalue of the 
Nuclear Three-Body Problem 


Certain relations exist among the lowest eigenvalues of 
the several Hamiltonian operators 


H=— }(Ai+4A2+As3) —Af(ri2) —Af(ris), (1) 
H’ = — }(Ai+A2+As3) —A2'f (riz) —As'f(ris), (2) 


H" = —}(O:+42+4:)—A"f(ri2) 
—A'f(ri3)—A'f(res), (3) 


H'"’ =— 3 (A: +Ae+As) —A fi rie) 
— Aa!" f (113) —As’"f (723), (4) 


when the potentials are related by the canditions 
Ai’ +A3' =3A" =A,'"+Aq'"+A3'" =2A. (5) 


(The function f(r) has the general form of a potential well. 
The normal state eigenfunctions and eigenvalues are 
designated by the symbols yo, Yo’, yo’, yo’, Eo, Eo’, Eo”, 
Ey", each primed symbol being associated with the 
correspondingly primed Hamiltonian.) 

Because of the symmetry properties of H’’, the eigen- 
function yo” is completely symmetrical in the coordinates 
of the three particles. A variational calculation with Jo” as 
the trial function yields at once the inequalities 


Eo”=Ev and Ey’ and Ey’. (6) 


The completely symmetrical Hamiltonian has the highest 
normal state eigenvalue. 
The eigenfunction Yo has the symmetry property 


Yo(1, 2, 3) =yYo(1, 3, 2). (7) 


Applying again the minimum property of the variational 
method, 
ESS +++ SWoH'podr =E, (8) 
and 
Ep’ SEpSEo"". (9) 


It is readily seen that the inequalities (6) and (9) are also 
valid for the case of Majorana forces between unlike 
particles. To obtain a lower limit to the eigenvalues Ey and 
Eo" set Ax’ =2A, A;'=0. Then if particle 3 and the center 
of gravity of the complete system have no kinetic energy 
the Hamiltonian H’ reduces to the expression 


—Ai2—2Af (ri) (10) 


and yo’, Eo’ are normal state eigenfunction and eigenvalue, 
respectively, of the equation 


[A+2A f(r) +Eo' yo’ =0. (11) 


Now the Schrédinger equation for the deuteron is 
simply 


{A+Af(r)+E(H2)} ¢=0, (12) 


which differs from (11) in that the potential well is only 
one-half as deep. In the physically interesting range of 
values for the radius of action of the neutron-proton force 
the eigenvalue of Eq. (12) is a small difference between 
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large potential and kinetic energy terms. It is thus evident 
without calculation that 


— Ey’ >>— E(H*)(=4me). (13) 


Actual numerical values for E»’ have been found by 
setting f(r)=e-ar* and determining A to fit the binding 
energy of the deuteron. The results are shown in the 
following table: 


(1 /a4) x 10" em 3.00 2.00 1.00 
— Eo'(me units) 24 36 84 


The lower limit Ep’ lies far below the upper limit found by 
variational calculations and has therefore little practical 
value. The same type of symmetry argument suffices ot 
demonstrate that inequalities similar to (6) and (9) exist 
also in the alpha-particle problem. 
EUGENE FEENBERG 
University of Wisconsin, 
January 7, 1936. 


The Production of High Velocity Particles in a Cyclotron 
by the Use of Multiphase Oscillators 


The voltage amplification of the present cyclotron is 
limited due to the building up of space charge within the 
duants. For a given accelerating potential, space charge 
may be cut down only at the expense of reduction in 
intensity of the ion beam. The space-charge density for a 
given final velocity with a given magnetic field is propor- 
tional to the number of times the ion must circulate before 
attaining that final velocity. The number of times the ion 
must circulate may be cut down by using higher potentials 
on the duants. This is undesirable, however, in many 
respects. 

In connection with work already begun here on the con- 
struction of an apparatus for the production of high 
velocity ions, it is proposed to experiment with oscillators 
producing 2-phase, 3-phase, --- m-phase potentials con- 
nected, respectively, to accelerating electrodes consisting 
of quadrants, sextants, --- 2m sectors. 

The accelerating field now used in the cyclotron may be 
considered to consist of two electric fields of equal intensity 
rotating in opposite directions. The ions are accelerated 
twice per revolution. With an n-phase oscillator connected 
to an m-sector cyclotron the ion will be accelerated 2n 
times per revolution. Thus, the space charge may be re- 
duced 1/nth of that in the present cyclotron for a given ion 
current and magnetic field. It is hoped that such an 
arrangement may lead to higher final velocities as well as 
to higher intensities. 

ROBERT J. Moon 
WituiAM D. HARKINS 


George Herbert Jones Chemical Laboratory, 
University of Chicago, 
January 18, 1936. 


LETTERS TO 


Quantum Mechanics and Physical Reality 


The present writer wishes to add a few words to the 


discussion of this subject in the current literature." 

Most of us have in the back of our minds an idea of 
‘physical reality” which is satisfactory in classical physics 
and which we tend to carry along into quantum mechanics. 
Discussion of a single free particle will be sufficient to 
bring out the essential point. We believe that the particle is 
somewhere (it has position) and that it is moving in some 
definite way (it has momentum) and so on for any other 
physical properties. In classical physics, these quantities 
would be, in principle, simultaneously measurable and the 
state of the particle could be described by giving their 
values. In spite of the fact that the values of these quanti- 
ties cannot be determined simultaneously (uncertainty 
principle), we still are inclined to believe in the state in 
which they all have definite values as an objective physical 
reality because we can measure any one of these quantities, 
chosen arbitrarily, and get a definite answer. 

Quantum mechanics does not concern itself with the 
“state”’ of the particle in the above sense but only with our 
knowledge of that ‘‘state.’’ The “state” is unknowable to 
the extent required by the uncertainty relations and any 
discussion of it belongs in the realm of speculation rather 
than physics. If we are to continue to employ this idea of 
the “state” of the particle, we must be careful not to base 
any physical arguments on it and we must not confuse it 
with our quantum-mechanical description of our knowledge 
about the particle. 

In quantum mechanics, our knowlege about the particle 
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is represented by a wave function, ¢, which is the solution 
of a wave equation. Our knowledge, based on measurement, 
determines the choice of the wave equation and the 
boundary conditions. This knowledge, embodied in ¢, may 
be used to predict the results of further measurements, 
Every observable, 6, has associated with it a quantum- 
mechanical operator, B, and the predicted values of } 
are given by /¢*Bedr. Usually this prediction takes the 
form of stating the probabilities for various possible values 
to be found by measurement.” 

Viewed in this light, the case discussed by Einstein, 
Podolsky and Rosen simply dissolves. Two systems have 
interacted and then separated. Nothing that we do to the 
first system after this affects the “state” of the second. 
But measurements on the first system affect our knowledge 
of the second and therefore affect the wave function which 
describes that knowledge. Different measurements on the 
first system give us different information about the second 
and therefore different wave functions and different 
predictions as to the results of measurement on the second 
system. 

The writer believes the ideas set forth here to be com- 
patible with Bohr’s point of view and hopes they will 
prove helpful in the dissemination of that point of view. 

HuGu C. WoLre 

College of the City of New York, 

December 30, 1935. 

1 Einstein, Podolsky and Rosen, Phys. Rev. 47, 777 (1935). N. Bohr, 
Phys. Rev. 48, 696 (1935). E. Schrédinger, Proc. Camb. Phil. Soc, 
31, 555 (1935). 

2 No attempt is made here to give an adequate discussion of this 


topic. The writer wishes to thank Dr. W. H. Furry for a preview of a 
forthcoming article which deals instructively with this question. 
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